hold on end legend('Frp') else subplot(1,3,1) qc=asin(n2./n1); qa=0:0.0001:qc; qb=asin(n1.*sin(qa)./n2); rs=-sin(qa-qb)./sin(qa+qb); rp=tan(qa-qb)./tan(qa+qb); ts=2.*cos(qa).*sin(qb)./sin(qa+qb); tp=2.*cos(qa).*sin(qb)./sin(qa+qb)./cos(qa-qb);