diff options
Diffstat (limited to '1535/CH5')
-rwxr-xr-x | 1535/CH5/EX5.1/Ch05Ex1.sci | 15 | ||||
-rwxr-xr-x | 1535/CH5/EX5.2/Ch05Ex2.sci | 9 | ||||
-rwxr-xr-x | 1535/CH5/EX5.3/Ch05Ex3.sci | 9 | ||||
-rwxr-xr-x | 1535/CH5/EX5.4/Ch05Ex4.sci | 12 | ||||
-rwxr-xr-x | 1535/CH5/EX5.5/Ch05Ex5.sci | 11 | ||||
-rwxr-xr-x | 1535/CH5/EX5.6/Ch05Ex6.sci | 10 |
6 files changed, 66 insertions, 0 deletions
diff --git a/1535/CH5/EX5.1/Ch05Ex1.sci b/1535/CH5/EX5.1/Ch05Ex1.sci new file mode 100755 index 000000000..defb9ad55 --- /dev/null +++ b/1535/CH5/EX5.1/Ch05Ex1.sci @@ -0,0 +1,15 @@ +// Scilab Code Ex5.1 : Polarization by reflection: Page-113 (2010)
+mu_g = 1.72; // Refractive index of glass
+mu_w = 4/3; // Refractive index of water
+// For polarization to occur on flint glass, tan(i) = mu_g/mu_w
+// Solving for i
+i = atand(mu_g/mu_w);
+printf("\nThe angle of incidence for complete polarization to occur on flint glass = %4.1f degrees", i);
+// For polarization to occur on water, tan(i) = mu_w/mu_g
+// Solving for i
+i = atand(mu_w/mu_g);
+printf("\nThe angle of incidence for complete polarization to occur on water = %5.2f degrees", i);
+
+// Result
+// The angle of incidence for complete polarization to occur on flint glass = 52.2 degrees
+// The angle of incidence for complete polarization to occur on water = 37.78 degrees
\ No newline at end of file diff --git a/1535/CH5/EX5.2/Ch05Ex2.sci b/1535/CH5/EX5.2/Ch05Ex2.sci new file mode 100755 index 000000000..545400db1 --- /dev/null +++ b/1535/CH5/EX5.2/Ch05Ex2.sci @@ -0,0 +1,9 @@ +// Scilab Code Ex5.2 : Percentage transmission of polarized light: Page-113 (2010)
+I0 = 1; // For simplicity, we assume the intensity of light falling on the second Nicol prism to be unity, watt per metre square
+theta = 30; // Angle through which the crossed Nicol is rotated, degrees
+I = I0*cosd(90-theta)^2; // Intensity of the emerging light from second Nicol, watt per metre square
+T = I/(2*I0)*100; // Percentage transmission of incident light
+printf("\nThe percentage transmission of incident light after emerging through the Nicol prism = %4.1f percent", T);
+
+// Result
+// The percentage transmission of incident light after emerging through the Nicol prism = 12.5 percent
\ No newline at end of file diff --git a/1535/CH5/EX5.3/Ch05Ex3.sci b/1535/CH5/EX5.3/Ch05Ex3.sci new file mode 100755 index 000000000..1cf8b3bcd --- /dev/null +++ b/1535/CH5/EX5.3/Ch05Ex3.sci @@ -0,0 +1,9 @@ +// Scilab Code Ex5.3 : Thickness of Quarter Wave Plate : Page-113 (2010)
+lambda = 6000e-008; // Wavelength of incident light, cm
+mu_e = 1.55; // Refractive index of extraordinary ray
+mu_o = 1.54; // Refractive index of ordinary ray
+t = lambda/(4*(mu_e - mu_o)); // Thickness of Quarter Wave plate of positive crystal, cm
+printf("\nThe thickness of Quarter Wave plate = %6.4f cm", t);
+
+// Result
+// The thickness of Quarter Wave plate = 0.0015 cm
\ No newline at end of file diff --git a/1535/CH5/EX5.4/Ch05Ex4.sci b/1535/CH5/EX5.4/Ch05Ex4.sci new file mode 100755 index 000000000..bf1d9b71f --- /dev/null +++ b/1535/CH5/EX5.4/Ch05Ex4.sci @@ -0,0 +1,12 @@ +// Scilab Code Ex5.4 : Behaviour of half wave plate for increased wavelength : Page-114 (2010)
+lambda = 1; // For simplicity, wavelength of incident light is assumed to be , cm
+mu_e = 1.55; // Refractive index of extraordinary ray
+mu_o = 1.54; // Refractive index of ordinary ray
+t = lambda/(2*(mu_e - mu_o)); // Thickness of Half Wave plate for given lambda, cm
+t_prime = 2*lambda/(2*(mu_e - mu_o)); // Thickness of Half Wave plate for twice lambda, cm
+printf("\nThe thickness of half wave plate is %2.1f times that of the quarter wave plate.", t/t_prime);
+printf("\nThe half wave plate behaves as a quarter wave plate for twice the wavelength of incident light.");
+
+// Result
+// The thickness of half wave plate is 0.5 times that of the quarter wave plate.
+// The half wave plate behaves as a quarter wave plate for twice the wavelength of incident light.
\ No newline at end of file diff --git a/1535/CH5/EX5.5/Ch05Ex5.sci b/1535/CH5/EX5.5/Ch05Ex5.sci new file mode 100755 index 000000000..1090d058b --- /dev/null +++ b/1535/CH5/EX5.5/Ch05Ex5.sci @@ -0,0 +1,11 @@ +// Scilab Code Ex5.5 : Phase retardation for quartz : Page-114 (2010)
+lambda = 500e-09; // Wavelength of incident light, m
+mu_e = 1.5508; // Refractive index of extraordinary ray
+mu_o = 1.5418; // Refractive index of ordinary ray
+t = 0.032e-03; // Thickness of quartz plate, m
+dx = (mu_e - mu_o)*t; // Path difference between E-ray and O-ray, m
+dphi = (2*%pi)/lambda*dx; // Phase retardation for quartz for given wavelength, rad
+printf("\nThe phase retardation for quartz for given wavelength = %5.3f pi rad", dphi/%pi);
+
+// Result
+// The phase retardation for quartz for given wavelength = 1.152 pi rad
\ No newline at end of file diff --git a/1535/CH5/EX5.6/Ch05Ex6.sci b/1535/CH5/EX5.6/Ch05Ex6.sci new file mode 100755 index 000000000..3553317ef --- /dev/null +++ b/1535/CH5/EX5.6/Ch05Ex6.sci @@ -0,0 +1,10 @@ +// Scilab Code Ex5.6 : Brewster angle at the boundary between two materials : Page-114 (2010)
+C = 52; // Critical angle for total internal reflection at a boundary between two materials, degrees
+// From Brewster's law, tand(i_B) = 1_mu_2
+// Also sind(C) = 1_mu_2, so that
+// tand(i_B) = sind(C), solving for i_B
+i_B = atand(sind(C)); // Brewster angle at the boundary, degrees
+printf("\nThe Brewster angle at the boundary between two materials = %2d degrees", i_B);
+
+// Result
+// The Brewster angle at the boundary between two materials = 38 degrees
\ No newline at end of file |