www.pudn.com > pinsanquxian.rar > ex1.m, change:2004-07-02,size:310b


function y=ex1(c)  
h=0.00276;cd=6100;cs=3300;j=1 
for w=1000000:1000000:1000000; 
    for c=1000; 
    b=w/c; 
     kdt=sqrt((w/cd)^2-b^2); 
     kst=sqrt((w/cs)^2-b^2); 
     y=tan(kst*h/2)/tan(kdt*h/2)+4*b^2*kst*kdt/(kst^2-b^2)^2; 
     if y<=0.001 
        z(j)=c;  
        j=j+1; 
     end 
     end 
end