defun my_gamma_in(S, ZL, Z02) { decl gammaL = (ZL-conj(Z02))/(ZL+Z02); decl gamma_in=S(1,1)+S(1,2)*S(2,1)*gamma_L/(1-S(2,2)*gamma_L); return gamma_in; }