根据一个模型，求解其中两个参数的变化关系，并通过图像表示出来，可是我总是画得不对，方程就是Plot里面的，但是显示不出来
Clear["Global`*"]
\[Kappa]1 = 49079/(6.28*3.68); \[Kappa]c =
5.58/(4*6.28*3.68); \[CapitalDelta]c = 1; \[CapitalDelta]1 = 1; \
\[Gamma]0 = 1*10^(-5); J0 = 9.1/36800;
g1 = 0.1*727/(6.28*3.68);
\[Delta] = pr + 1;
Plot[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) == Epu^2, {Epu, 0, 50}]

隐函数的绘图可采用"ContourPlot"，“Plot”不能绘制隐函数图像。
Clear["Global`*"]
\[Kappa]1 = 49079/(6.28*3.68); \[Kappa]c =
5.58/(4*6.28*3.68); \[CapitalDelta]c = 1; \[CapitalDelta]1 = 1; \
\[Gamma]0 = 1*10^(-5); J0 = 9.1/36800;
g1 = 0.1*727/(6.28*3.68);
\[Delta] = pr + 1;
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) == Epu^2, {Epu, 0, 50}, {na, -10,
10}]