rg=11077;l=.56;
cg=24.0e-9;
f=20e+3;
om=2*3.14*f;
re=1/rg;
im=2*pi*f*cg;
yg=re+im*i;
zg=1/yg;
rs=167;
cs=6e-9;
f=20e+3;
om=2*3.14*f;
re=1/rs;
im=2*pi*f*cs;
ys=re-im*i;
zs=1/ys;
zc=sqrt(zs*zg);
tgh=sqrt(zs/zg);
t1=abs(tgh);
fi=angle(tgh);
alfa=(1/(4*l))*log((1+t1^2+2*t1*cos(fi))/(1+t1^2-2*t1*cos(fi)));
for k=1:4 
   beta(k)=(1/(2*l))*(atan(2*t1*sin(fi)/(1-t1^2))+(k-1)*pi);

gam(k)=alfa+beta(k)*i;
t2(k)=gam(k)*zc;
t3(k)=gam(k)/zc;
ll(k)=(1/om)*(imag(t2(k)));
cap(k)=(1/om)*(imag(t3(k)));
end

