% WELL DONE

lamda=532 * 10^-9;
d=38 * 10^-6;
theta_r=45 * pi/180;
theta_s=-45 * pi/180;
dn=[];


%efficiency_s=100%
%for m=1:5
%   theta_s(m)=theta_r+acos(m/(m+1/2));
%   dn(m)=1/d*(m+1/2)*lamda*(cos(theta_r)*cos(theta_s(m)))^.5;
%end;


efficiency_p=100%
for m=1:1
   theta_s(m)=theta_r+acos((m-1/2)/m);
   dn(m)=1/d*m*lamda*(cos(theta_r)*cos(theta_s(m)))^.1;
end;

%theta_s *pi/180
dn