Subscript[, 0] = 0.14;
Subscript[\[Sigma], v] = 5.0*10^22;
e = 1.602*10^-19;
m = 9.1092897*10^-31;
c = 3*10^8;
Subscript[\[Epsilon], 0] = 8.85*10^-12;
\[HBar] = 1.055*10^-34;
T = 100*10^3;
R = 10*10^-9;
P = 150*10^3;
Subscript[\[Omega], 0] = 3.6*10^14;
Subscript[n, r] = 3.2;
\[Mu] = 4*\[Pi]*10^-7;
n = 1;
Subscript[E, n] = (2*n + 3/2)*\[HBar]*Subscript[\[Omega], 0];
Subscript[E, g][P, T] = 1519 - (0.5405 T^2)/(T + 204) + 10.7 P;
Subscript[m, e] =
m*(1 + 7510/(Subscript[E, g][P, T] + 341) + 15020/
Subscript[E, g][P, T])^-1;
Subscript[m, h] = m*(0.09 - 0.20*10^-3*P - 3.55*10^-5*T);
\[Mu] = Subscript[m, e] + Subscript[m, h];
w = Subscript[\[Omega], 0]/(1 - 2*P*0.14);
\[Beta] = Sqrt[(\[Mu]*w)/\[HBar]];
Subscript[N, n] = (1/\[Beta] Sqrt[\[Pi]] 2^(2*n) ((2*n + 1)!))^(-1/2);
z = 2;
\[Phi] = Subscript[N, n]*
Exp[-N[(1/2*\[Beta]^2*z^2), $MachinePrecision]]*
HermiteH[2*n + 1, \[Beta]*z];

你看看这帖：
https://tieba.baidu.com/p/6784934028