getf('multilayer.sci');
angles = linspace(0,90,90);
[re,te]=coefRTpuissances(1,1,angles*%pi/180,1.5,1,[]);
[rm,tm]=coefRTpuissances(1,1,angles*%pi/180,1.5,0,[]);
            
plotA=scf();
plot2d(angles,re,style=4);
plot2d(angles,rm,style=5);
legend("TE","TM");
xtitle("n1=1, n2=1.5","theta1","R");
            

n=sqrt(1.5);
[re,te]=coefRTpuissances(1,1,angles*%pi/180,1.5,1,[n,0.25]);
[rm,tm]=coefRTpuissances(1,1,angles*%pi/180,1.5,0,[n,0.25]);  
            
plotB=scf();
plot2d(angles,re,style=4);
plot2d(angles,rm,style=5);
legend("TE","TM");
xtitle("","theta1","R");
            
