www.pudn.com > algorithms.zip > horn.m, change:2010-05-14,size:502b


permability=input('input the permability'); 
epsilon=input('input the epsilon'); 
frequency=input('input the frequency'); 
landa=300000000/frequency; 
a=input('input a'); 
eta=sqrt(epsilon/permability); 
z0=inline('((eta*(landa/pi*b1))*((2*b1/landa)*((cos(pi*a/landa))^2)*log(csc(pi*b1/(2*b))))+((pi*a1)/(2*landa))+(0.25*sin((2*pi*a1)/landa))+((b1/b)*((cos((pi*a1)/landa)^2)/sin((pi/landa)*(a-a1)))*(((pi*(a-a1))/(2*landa))-(0.25*sin((2*pi/landa)*(a-a1))))))','b1','a1','b','a','eta','landa')