m = 20;
\[HBar] = 1.05457;
\[Omega]1 = 50
G/(2*Pi) == -12;
Subscript[\[CapitalGamma], m]/(2*Pi) == 41;
\[Kappa]/(2*Pi) == 15;
Subscript[\[CapitalOmega], m]/(2*Pi) == 51.8;
\[CapitalDelta] = -Subscript[\[CapitalOmega], m];
\[CapitalDelta]1 = \[CapitalDelta] - G*x1;
Subscript[\[Eta], c] = 1/2;
\[Epsilon]1 = Sqrt[P1/(\[HBar]*\[Omega]1)];
a1 = (Sqrt[Subscript[\[Eta], c]*\[Kappa]]/
((-I)*\[CapitalDelta]1*
TextCell["+\[Kappa]/2"]))*\[Epsilon]1;
f[P1_] = P1 /. Solve[x1 ==
-((\[HBar]*G*Abs[a1]^2)/
(m*Subscript[\[CapitalOmega], m]^2)), P1];
Plot[f[x], {x, 0, 200}, PlotRange ->
{{0, 200}, {0, 8}}]
参照文献，写了下来，运行没报错也没出图

https://tieba.baidu.com/p/2964416898 第19条