二楼附上程序。
Clear["Global`*"]
\[Kappa]1 = 49079/(6.28*3.68); \[Kappa]c = 5.58/(4*6.28*3.68);
\[CapitalDelta]1 = 1; \[Gamma]0 = 1*10^(-7); J0 = 9.1/36800;
g1 = 0*727/(6.28*3.68); aa = 100^2;
p1 = ContourPlot[
na*((\[CapitalDelta]c - 2*J0^2*na -
g1^2*\[CapitalDelta]1/(\[CapitalDelta]1^2 + (\[Kappa]1/
2)^2))^2 + (g1^2*\[Kappa]1/(\[CapitalDelta]1^2 + (\
\[Kappa]1/2)^2) + \[Kappa]c)^2/4) == aa, {\[CapitalDelta]c, -1,
2}, {na, 0, 11000000}, ContourStyle -> {Blue},
PlotLegends -> {"g=0"}]
Clear["Global`*"]
\[Kappa]1 = 49079/(6.28*3.68); \[Kappa]c = 5.58/(4*6.28*3.68);
\[CapitalDelta]1 = 1; \[Gamma]0 = 1*10^(-7); J0 = 9.1/36800;
g1 = 0*727/(6.28*3.68); aa = 100^2;
p1 = ContourPlot[
na*((\[CapitalDelta]c - 2*J0^2*na -
g1^2*\[CapitalDelta]1/(\[CapitalDelta]1^2 + (\[Kappa]1/
2)^2))^2 + (g1^2*\[Kappa]1/(\[CapitalDelta]1^2 + (\
\[Kappa]1/2)^2) + \[Kappa]c)^2/4) == aa, {\[CapitalDelta]c, -1,
2}, {na, 0, 11000000}, ContourStyle -> {Blue},
PlotLegends -> {"g=0"}]