close all
n1=input('Brytningsindex på infallssida?  ');
n2=input('Brytningsindex på brytningssida?  ');
i=[0:0.2:90]*pi/180+eps;
b=asin(n1/n2 .*sin(i));
rtm=-tan(i-b) ./tan(i+b);
rte=-sin(i-b) ./sin(i+b);
Rtm=abs(rtm .^2);
Rte=abs(rte .^2);
svar=input('Resultatet ligger i rte,rtm,Rte och Rtm. Vill Du ha en plot? y/n ','s');
if (svar=='n') | (svar=='N') | (isempty(svar)) break
end
figure(1)
set(1,'position',[10 10 500 200])
plot(i*180/pi,abs(rte),'b',i*180/pi,abs(rtm),'r')
grid
xlabel('Infallsvinkel (grader)')
ylabel('Amplitudreflektans (belopp) ')
text(10,0.8,'E-fält')
figure(2)
set(2,'position',[10 240 500 200])
plot(i*180/pi,Rte,'b',i*180/pi,Rtm,'r')
grid
xlabel('infallsvinkel (grader)')
ylabel('Intensitetsreflektans ')
text(10,0.8,'Intensitet')