您所在位置：主题：求助，为什么我的程序部能运行？
yaoguoling
[专家分：0] 发布于 2010-09-05 17:04:00
E=1:10:100;
w=1;
m=2;
Eg=7.8;
e=1.6;
r=(w*(m*Eg)^1/2)./(e*E);
r1=r.^2./(1+r.^2);
r2=1./(1+r.^2);
[k,e]=ellipke(r1);
[k,e]=ellipke(r2);
x=[2*Eg.*E(r2)]/(pi*h*w*r1.^2);
z=pi*(2*int(x+1)-2*x+n)^1/2;
H(z)=int(exp(y^2-z^2),y,0,z);
Q(r,x)=(pi/(2*K(r2)).^1/2*(symsum(-(n*pi*(K(r2)-E(r2)))./E(r1),0,inf)...
*H(z);
W(E)=2*w*(w*m)^3/2*Q(r,x)*exp(-(pi*int(x+1).*(K(r1)-E(r2))./(E(r2))...
*(9*pi*.(r1.^3/4);
w=1;
m=2;
Eg=7.8;
e=1.6;
r=(w*(m*Eg)^1/2)./(e*E);
r1=r.^2./(1+r.^2);
r2=1./(1+r.^2);
[k,e]=ellipke(r1);
[k,e]=ellipke(r2);
x=[2*Eg.*E(r2)]/(pi*h*w*r1.^2);
z=pi*(2*int(x+1)-2*x+n)^1/2;
H(z)=int(exp(y^2-z^2),y,0,z);
Q(r,x)=(pi/(2*K(r2)).^1/2*(symsum(-(n*pi*(K(r2)-E(r2)))./E(r1),0,inf)...
*H(z);
W(E)=2*w*(w*m)^3/2*Q(r,x)*exp(-(pi*int(x+1).*(K(r1)-E(r2))./(E(r2))...
*(9*pi*.(r1.^3/4);
