f=10000;omega=2*pi*f;
yg=10^(-3)/9.9+j*2*pi*10000*22.0*10^(-9);
ys=1/151-j*2*pi*10000*9.38*10^(-9);
zc=1/sqrt(yg*ys);
tt=sqrt(yg/ys);l=.55;
t=abs(tt);fi=angle(tt);alfa=1/(4*l)*log((1+t*t+2*t*cos(fi)))/(1+t*t-2*t*cos(fi));
for k=1:3
beta(k)=(1/(4*l))*(atan(2*t*sin(fi)/(1-t*t))+(k-1)*pi);
gam(k)=alfa+j*beta(k);
r(k)=real(zc*gam(k));
l(k)=imag(zc*gam(k))/omega;
g=real(zc/gam(k));
c=imag(zc/gam(k))/omega;
end
