s11=0.57*exp(j*(-80)/180*pi);
s12=0.016*exp(j*(+169)/180*pi);
s21=1.79*exp(j*(+74)/180*pi);
s22=0.85*exp(j*(-26)/180*pi);
s_param=[s11,s12;s21,s22];
delta=abs(det(s_param))
k=(1-abs(s11)^2-abs(s22)^2+delta^2)./(2*abs(s12.*s21))

B1=1+abs(s11)^2-abs(s22)^2-delta^2;
C1=s11-conj(s22)*delta;
B2=1+abs(s22)^2-abs(s11)^2-delta^2;
C2=s22-conj(s11)*delta;

Gamma_ms=(B1-sqrt(B1^2-4*abs(C1)^2))./(2*C1)
Gamma_ml=(B2-sqrt(B2^2-4*abs(C2)^2))./(2*C2)

zs=(1+Gamma_ms)/(1-Gamma_ms);
zl=(1+Gamma_ml)/(1-Gamma_ml);
z1=1/(1/real(zs)-1);
z2=1/(1/real(zl)-1);
Z1=z1*50
Z2=z2*50
L=50*imag(zs)/(500*10^6*2*pi)*10^9
C=50/(imag(-zl)*(500*10^6*2*pi))*10^12