wo2gyy
2012-11-24, 05:44
%k2=1,c=2,w1=5,w2=3,q'2,1=5/1372,q2,1=3/686
k2=1;
y2=0;
y4=0;
y6=0;
for k1=0:50;
for w=0:3;
for i=0:w;
for j=0: (6-2*w);
x1=nchoosek(6-2*w,j);
y1=((-1)^j)*x1*((1-j*5/1372)^k1);
y2=y1+y2;
end;
x2=nchoosek(w,i);
y3=((-1)^i)*x2*((1-i*3/686)^(k2-1))*y2;
y4=y3+y4;
end;
x3=nchoosek(3,w);
y5=x3*y4;
y6=y6+y5;
end;
z=0.5*y6;
k1=0:50;
figure(2);hold on;semilogy(k1,z);plot(k1,z);grid on
end;
k2=1;
y2=0;
y4=0;
y6=0;
for k1=0:50;
for w=0:3;
for i=0:w;
for j=0: (6-2*w);
x1=nchoosek(6-2*w,j);
y1=((-1)^j)*x1*((1-j*5/1372)^k1);
y2=y1+y2;
end;
x2=nchoosek(w,i);
y3=((-1)^i)*x2*((1-i*3/686)^(k2-1))*y2;
y4=y3+y4;
end;
x3=nchoosek(3,w);
y5=x3*y4;
y6=y6+y5;
end;
z=0.5*y6;
k1=0:50;
figure(2);hold on;semilogy(k1,z);plot(k1,z);grid on
end;