r1 = 0.9;

t1 = sqrt(1-r1^2);



r2 = 0.8;

t2 = sqrt(1-r2^2);



phi = pi/3;



S = 1/(1-r1*r2*exp(-j*phi))*[r1-r2*exp(-j*phi) -t1*t2*exp(-j*phi/2);

                             -t1*t2*exp(-j*phi/2) r2-r1*exp(-j*phi)];

                         