f=80000;omega=2*pi*f;
yg=10^(-3)/.287+j*omega*2.46*10^(-9);
ys=1/81+j*omega*9.75*10^(-9);
zc=1/sqrt(yg*ys);
tt=sqrt(yg/ys);l1=.55;
t=abs(tt);fi=angle(tt);alfa=1/(4*l1)*log((1+t*t+2*t*cos(fi))/(1+t*t-2*t*cos(fi)));bet=ones(4,1);
k=1
bet=(1/(2*l1))*(atan((2*t*sin(fi))/(1-t*t))+(k-1)*pi);

r=ones(3,1);
l=ones(3,1);
c=ones(3,1);
g=ones(3,1);
for k=1:3
bet(k,1)=(1/(2*l1))*(atan((2*t*sin(fi))/(1-t*t))+(k-1)*pi);
gam(k,1)=alfa+j*bet(k,1);
r(k,1)=real(zc*gam(k,1));
l(k,1)=imag(zc*gam(k,1))/omega;
g(k,1)=real((gam(k,1))/zc);
c(k,1)=imag((gam(k,1))/zc)*(1/omega);
end
