//Chapter 3, Prblem 16, figure 3.32
clc
//from figure 3.32
zo=50                       //in ohm
zl=50                       //in ohm

//calculating the S parameter
z1=(zo*zl)/(zo+zl)
s11=(z1-zo)/(z1+zo)
z2=(zo*zl)/(zo+zl)
s22=(z2-zo)/(z2+zo)
s21=(2*z1)/(50+z1)
s12=(2*z2)/(50+z2)

s=[s11 s12;s21 s22]

printf("s11 (magnitude) = %.3f \n\t(angle) = 180 degree\n\n",-s11)
printf("s12 (magnitude) = %.3f \n\t(angle) = 0 degree\n\n",s12)
printf("s21 (magnitude) = %.3f \n\t(angle) = 0 degree\n\n",s21)
printf("s22 (magnitude) = %.3f \n\t(angle) = 180 degree\n\n",-s22)
disp(s,"S (magnitude) = ")