f=2000;omega=2*pi*f;
yg=10^(-3)/220+j*2*pi*2000*23.7*10^(-9);
ys=1/146-j*2*pi*2000*8.85*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)));
 
   k=0
beta=1/(4*l)*(atan((2*t*sin(fi))/(1-t*t))+1*pi);
gam=alfa+j*beta;
r=real(zc*gam);
l=imag(zc*gam)/omega;
g=real(zc/gam);
c=imag(zc/gam)/omega;
