%Deviation of a ray passing through a prism as a function of the angle of
%incidence

%Apex Angle in degrees
alpha = 60;

%Index of Refraction
%for purple
n = 1.525;

%Incident angle in degrees
i_angle = 1:0.5:89;

%Deviation angle
delta = i_angle - alpha + asind(sind(alpha.*sqrt(n^2 - sind(i_angle).^2)) - cosd(alpha).*sind(i_angle));

figure;
plot(i_angle, delta);
xlabel('Incident angle in degrees');
ylabel('Deviation in degrees');

hold on;

%Prism 2, for purple
n = 1.674;
delta = i_angle - alpha + asind(sind(alpha.*sqrt(n^2 - sind(i_angle).^2)) - cosd(alpha).*sind(i_angle));
plot(i_angle, delta, 'r');
legend('Prism 1 for Purple Wavelength, n = 1.525', 'Prism 2 for Purple Wavelength, n = 1.674');