Bạn đang xem bản rút gọn của tài liệu. Xem và tải ngay bản đầy đủ của tài liệu tại đây (1.42 MB, 73 trang )
e0=1;m0=9.1093892-31;h=1.05459*10^(-34);
m1=.067*m0;e=2.07*e0;n0=1e23;
k=1.3807e-23;
a11=2*pi^2*e^2*k.*T*V^2./((c*sqrt(Xinf)*r0*vs^2*m1^2.*omeg2.^3).*h.^2);
y=0;
for n=1:3
for t=1:mt
for n1=1:3
if (Imnn5(t,n,n1)~=0)&(Imnn5(t,n,n1)~=inf);
al=a11.*(D(0,1,omeg1,omeg2,T,d,E02,n,n1)...
-D(0,-2,omeg1,omeg2,T,d,E02,n,n1)...
-1/2.*(H(0,1,omeg1,omeg2,T,d,E01,E02,gama,n,n1)...
-H(0,-1,omeg1,omeg2,T,d,E01,E02,gama,n,n1))...
+3/32.*(G(0,1,omeg1,omeg2,T,d,E01,E02,gama,n,n1)...
-G(0,-1,omeg1,omeg2,T,d,E01,E02,gama,n,n1))...
+1/4.*(H(-1,1,omeg1,omeg2,T,d,E01,E02,gama,n,n1)...
-H(-1,-1,omeg1,omeg2,T,d,E01,E02,gama,n,n1)...
+H(1,1,omeg1,omeg2,T,d,E01,E02,gama,n,n1)...
-H(1,-1,omeg1,omeg2,T,d,E01,E02,gama,n,n1))...
-1/16.*(G(-1,1,omeg1,omeg2,T,d,E01,E02,gama,n,n1)...
-G(-1,-1,omeg1,omeg2,T,d,E01,E02,gama,n,n1)...
+G(1,1,omeg1,omeg2,T,d,E01,E02,gama,n,n1)...
-G(1,-1,omeg1,omeg2,T,d,E01,E02,gama,n,n1))...
68
+1/64.*(G(-2,1,omeg1,omeg2,T,d,E01,E02,gama,n,n1)...
-G(-2,-1,omeg1,omeg2,T,d,E01,E02,gama,n,n1)...
+G(2,1,omeg1,omeg2,T,d,E01,E02,gama,n,n1)...