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];


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=Gamma_ms*50
ZL=Gamma_ml*50
zl=1/(1/real(ZL)-1/50)
zs=1/(1/real(ZS)-1/50)