%% Calculates incident angle for proper coupling into the disc via prism

% angle of the prism faces in degrees
prism_angle_in_degrees = 45; 

prism_angle = prism_angle_in_degrees*pi/180; 

%% prism index of refraction
% Rutile (TiO2) see http://refractiveindex.info/?group=CRYSTALS&material=TiO2
n_rutile_o = 2.4885;  % p - polarization
n_rutile_e = 2.75324; % s - polarization

n_prism=n_rutile_o ; % for horizontal or p polarization

%% disk material index of refraction
% Magnesium Fluoride (MgF2)  see http://refractiveindex.info/?group=CRYSTALS&material=MgF2
n_MgF2_o = 1.3751;
n_MgF2_e = 1.38679;

% Measured 
n_MgF2_p = 1.465; % p-polarization


n_disk=n_MgF2_p;
coupling_description='Rutile prism, MgF_{2} disk, p-polarization';

% calculate coupling angle and draw coupling solution with annotations
prism_disk_coupling(prism_angle_in_degrees, n_disk, n_prism, coupling_description);