From a08ed173692ce16b79002733003049ec32a02485 Mon Sep 17 00:00:00 2001 From: Eugeniy Mikhailov Date: Thu, 28 Feb 2013 14:58:35 -0500 Subject: moved octave using code to separate folder --- beam_tracing/octave/demo_prism_disk_coupling.m | 29 ++++++++++++++++++++++++++ 1 file changed, 29 insertions(+) create mode 100644 beam_tracing/octave/demo_prism_disk_coupling.m (limited to 'beam_tracing/octave/demo_prism_disk_coupling.m') diff --git a/beam_tracing/octave/demo_prism_disk_coupling.m b/beam_tracing/octave/demo_prism_disk_coupling.m new file mode 100644 index 0000000..e809968 --- /dev/null +++ b/beam_tracing/octave/demo_prism_disk_coupling.m @@ -0,0 +1,29 @@ +%% 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); + -- cgit v1.2.3