使用的程序如下：
res11 = {}; res22 = {};
Do[{h = 0.01 i ; Subscript[\[CurlyEpsilon], b] = 2.3;
Subscript[\[CurlyEpsilon], \[Infinity]] = 4.6;
a = 7; \[HBar] = 6.58*10^-16;
Subscript[\[CurlyEpsilon], 0] = 8.85*10^-21; \[Mu] = 3.8*1.6*10^-19;
Subscript[\[Mu], r] = Sin[Pi/2]*Cos[0]*\[Mu];
Subscript[\[Mu], \[Theta]] = Sin[Pi/2]*Sin[0]*\[Mu];
Subscript[\[Mu], \[CapitalPhi]] = Cos[Pi/2]*\[Mu];
Subscript[\[Omega], p] = 9/\[HBar];
Subscript[\[Gamma], p] = 0.1/\[HBar];
Subscript[\[Omega], 1] = Subscript[\[Omega], p]/Sqrt[
Subscript[\[CurlyEpsilon], \[Infinity]] +
Subscript[\[CurlyEpsilon], b]*2];
Subscript[\[Omega], n] = Subscript[\[Omega], p]/Sqrt[
Subscript[\[CurlyEpsilon], \[Infinity]] +
Subscript[\[CurlyEpsilon], b]*(n + 1)/n];
x = a^(2 n + 1)/(a + h)^(2 n + 4);
y = 1 + 1/(2 n);
u = Subscript[\[Omega], p]/(
4*Pi*Subscript[\[CurlyEpsilon], 0]*\[HBar]);
c = (Subscript[\[Omega], n]/Subscript[\[Omega], p])^3;
f = x*y*u*c;
Subscript[g, n] =
Subscript[\[Mu], r]^2*(n + 1)^2*f +
Subscript[\[Mu], \[Theta]]^2*(n*(n + 1))/2*f +
Subscript[\[Mu], \[CapitalPhi]]^2*(n*(n + 1))/2*f;
g = Subscript[g, n]*1/(1.6*10^-19);
m = Sum[g, {n, 2, 10}];
Subscript[g, M] = Sqrt[m];
Subscript[\[Omega], M1] = Subscript[\[Omega], n]*g;
Subscript[g, d] =
Sqrt[(3*Subscript[\[Mu], r]^2*a^3*Subscript[\[Omega], 1]^3)/(
2*Pi*\[HBar]*Subscript[\[CurlyEpsilon], 0]*(a + h)^6*
Subscript[\[Omega], p]^2) + (
3*Subscript[\[Mu], \[Theta]]*Subscript[\[Mu], \[Theta]]*a^3*
Subscript[\[Omega], 1]^3)/(
8*Pi*\[HBar]*Subscript[\[CurlyEpsilon], 0]*(a + h)^6*
Subscript[\[Omega], p]^2) + (
3*Subscript[\[Mu], \[CapitalPhi]]*
Subscript[\[Mu], \[CapitalPhi]]*a^3*Subscript[\[Omega], 1]^3)/(
8*Pi*\[HBar]*Subscript[\[CurlyEpsilon], 0]*(a + h)^6*
Subscript[\[Omega], p]^2)]*Sqrt[1/(1.6*10^-19)];
Subscript[\[Omega], M] = Sum[Subscript[\[Omega], M1], {n, 2, 10}]/m;
Subscript[\[Gamma], QE] = Subscript[\[Gam

ma], p] - 4*Su

代码不齐全啊