diff options
Diffstat (limited to 'Engineering_Physics_by_V_Yadav/4-Polarization.ipynb')
| -rw-r--r-- | Engineering_Physics_by_V_Yadav/4-Polarization.ipynb | 850 |
1 files changed, 850 insertions, 0 deletions
diff --git a/Engineering_Physics_by_V_Yadav/4-Polarization.ipynb b/Engineering_Physics_by_V_Yadav/4-Polarization.ipynb new file mode 100644 index 0000000..69cf5ed --- /dev/null +++ b/Engineering_Physics_by_V_Yadav/4-Polarization.ipynb @@ -0,0 +1,850 @@ +{ +"cells": [ + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "# Chapter 4: Polarization" + ] + }, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 4.10: Thickness_of_quarter_wave_plate.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"// Scilab Code Ex4.10:: Page-4.23 (2009)\n", +"clc; clear;\n", +"mu_o = 1.658; // Refractive index of ordinary wave\n", +"mu_e = 1.486; // Refractive index of extraordinary wave\n", +"lambda = 5893e-008; // Wavelength of light used, m\n", +"// As (mu_o - mu_e)*t = lambda/4, solving for t\n", +"t = lambda/(4*(mu_o - mu_e)); // Thickness of quarter-wave plate, cm\n", +"\n", +"printf('\nThe thickness of quarter-wave plate = %3.1e cm', t);\n", +"\n", +"// Result \n", +"// The thickness of quarter-wave plate = 8.6e-005 cm \n", +"\n", +" " + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 4.11: Least_thickness_of_plate_for_which_emergent_beam_is_plane_polarised.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"// Scilab Code Ex4.11:: Page-4.23 (2009)\n", +"clc; clear;\n", +"mu_o = 1.5442; // Refractive index of ordinary wave\n", +"mu_e = 1.5533; // Refractive index of extraordinary wave\n", +"lambda = 5000e-008; // Wavelength of light used, m\n", +"// As (mu_o - mu_e)*t = lambda/4, solving for t\n", +"t = lambda/(4*(mu_e - mu_o)); // Least thickness of plate for which emergent beam is plane polarised, cm\n", +"\n", +"printf('\nThe least thickness of plate for which emergent beam is plane polarised = %4.2e cm', t);\n", +"\n", +"// Result \n", +"// The least thickness of plate for which emergent beam is plane polarised = 1.37e-003 cm \n", +"\n", +" " + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 4.12: Difference_in_refractive_indices_of_rays.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"// Scilab Code Ex4.12:: Page-4.23 (2009)\n", +"clc; clear;\n", +"lambda = 5893e-008; // Wavelength of light used, m\n", +"t = 0.005; // Thickness of the crystal, cm\n", +"// As for quarter wave plate, mu_diff*t = (mu_o - mu_e)*t = lambda/4, solving for mu_diff\n", +"mu_diff = lambda/(4*t); // The difference in refractive indices of rays, cm\n", +"printf('\nThe least thickness of plate for which emergent beam is plane polarised = %4.2e cm', mu_diff);\n", +"\n", +"// Result \n", +"// The least thickness of plate for which emergent beam is plane polarised = 2.95e-003 cm \n", +"\n", +" " + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 4.13: The_thickness_of_a_half_wave_plate.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"// Scilab Code Ex4.13:: Page-4.24 (2009)\n", +"clc; clear;\n", +"mu_o = 1.54; // Refractive index of ordinary wave\n", +"mu_e = 1.45; // Refractive index of extraordinary wave\n", +"lambda = 5500e-008; // Wavelength of light used, m\n", +"// As for a half wave plate, (mu_o - mu_e)*t = lambda/4, solving for t\n", +"t = lambda/(2*(mu_o - mu_e)); // The thickness of a half wave plate for wavelength, cm\n", +"\n", +"printf('\nThe thickness of a half wave plate for wavelength = %4.2e cm', t);\n", +"\n", +"// Result \n", +"// The thickness of a half wave plate for wavelength = 3.06e-004 cm \n", +"\n", +" " + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 4.14: The_thickness_of_a_quarter_wave_plate.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"// Scilab Code Ex4.14:: Page-4.24 (2009)\n", +"clc; clear;\n", +"mu_o = 1.55; // Refractive index of ordinary wave\n", +"mu_e = 1.52; // Refractive index of extraordinary wave\n", +"lambda = 5500e-008; // Wavelength of light used, m\n", +"// As for a half wave plate, (mu_o - mu_e)*t = lambda/4, solving for t\n", +"t = lambda/(4*(mu_o - mu_e)); // The thickness of a quarter wave plate for wavelength, cm\n", +"\n", +"printf('\nThe thickness of a quarter wave plate for wavelength = %4.2e cm', t);\n", +"\n", +"// Result \n", +"// The thickness of a quarter wave plate for wavelength = 4.58e-004 cm \n", +"\n", +" " + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 4.15: The_thickness_of_a_half_wave_plate_quartz.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"// Scilab Code Ex4.15:: Page-4.24 (2009)\n", +"clc; clear;\n", +"mu_o = 1.51; // Refractive index of ordinary wave\n", +"mu_e = 1.55; // Refractive index of extraordinary wave\n", +"lambda = 6000e-008; // Wavelength of light used, m\n", +"// As for a half wave plate, (mu_o - mu_e)*t = lambda/4, solving for t\n", +"t = lambda/(2*(mu_e - mu_o)); // The thickness of a quarter wave plate for wavelength, cm\n", +"\n", +"printf('\nThe thickness of a half wave plate quartz = %4.2e cm', t);\n", +"\n", +"// Result \n", +"// The thickness of a half wave plate quartz = 7.50e-004 cm \n", +"\n", +" " + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 4.16: Difference_between_refractive_indices.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"// Scilab Code Ex4.16:: Page-4.24 (2009)\n", +"clc; clear;\n", +"lambda = 5890e-008; // Wavelength of light used, m\n", +"t = 7.5e-004; // Thickness of the crystal, cm\n", +"// As for quarter wave plate, mu_diff*t = (mu_e - mu_o)*t = lambda/4, solving for mu_diff\n", +"mu_diff = lambda/(4*t); // The difference in refractive indices of rays, cm\n", +"printf('\nThe difference between refractive indices = %6.4f cm', mu_diff);\n", +"\n", +"// Result \n", +"// The difference between refractive indices = 0.0196 cm \n", +"\n", +" " + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 4.17: Specific_rotation_of_superposition.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"// Scilab Code Ex4.17:: Page-4.34 (2009)\n", +"clc; clear;\n", +"theta = 15.2; // Angle through which plane of polarization is rotated, degrees\n", +"c = 0.2; // Concentration of sugar, g/cc\n", +"l = 25; // Length of sugar, cm\n", +"S = 10*theta/(l*c); // Specific rotation of superposition, degrees\n", +"\n", +"printf('\nThe specific rotation of superposition = %4.1f cm', S);\n", +"\n", +"// Result \n", +"// The specific rotation of superposition = 30.4 cm \n", +"\n", +" " + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 4.18: Strength_of_sugar_solution.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"// Scilab Code Ex4.18: : Page-4.34 (2009)\n", +"clc; clear;\n", +"theta = 15.2; // Angle through which plane of polarization is rotated, degrees\n", +"S = 65; // Specific rotation of sugar solution, degrees\n", +"l = 15; // Length of sugar, cm\n", +"// As S = 10*theta/(l*c), solving for c\n", +"c = 10*theta/(l*S); // Concentration of sugar, g/cc\n", +"\n", +"printf('\nThe strength of sugar solution = %4.2f g/cc', c);\n", +"\n", +"// Result \n", +"// The strength of sugar solution = 0.16 g/cc \n", +"\n", +" " + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 4.19: Quantity_of_sugar_contained_in_the_tube_in_the_form_of_solution.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"// Scilab Code Ex4.19:: Page-4.34 (2009)\n", +"clc; clear;\n", +"theta = 15; // Angle through which plane of polarization is rotated, degrees\n", +"S = 69; // Specific rotation of sugar solution, degrees\n", +"l = 10; // Length of sugar, cm\n", +"V = 50; // Volume of the tube, cc\n", +"// As S = 10*theta/(l*c), solving for c\n", +"c = 10*theta/(l*S); // Concentration of sugar, g/cc\n", +"M = c*V; // Mass of sugar in solution, g\n", +"\n", +"printf('\nThe quantity of sugar contained in the tube in the form of solution = %5.2f g', M);\n", +"\n", +"// Result \n", +"// The quantity of sugar contained in the tube in the form of solution = 10.87 g \n", +"\n", +" " + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 4.1: Refractive_index_of_the_material.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"// Scilab Code Ex4.1:: Page-4.5 (2009)\n", +"clc; clear;\n", +"ip = 60; // Polarizing angle, degrees\n", +"mu = tand(ip); // Refractive index of the material from Brewster's law \n", +"printf('\nThe refractive index of the material = %5.3f', mu);\n", +"\n", +"// Result \n", +"// The refractive index of the material = 1.732 " + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 4.20: Specific_rotation_of_sugar_solution_from_the_given_data.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"// Scilab Code Ex4.20:: Page-4.35 (2009)\n", +"clc; clear;\n", +"theta = 8; // Angle through which plane of polarization is rotated, degrees\n", +"M = 10; // Amount of sugar, g\n", +"l = 14; // Length of the tube, cm\n", +"V = 44; // Volume of sugar solution, cc\n", +"c = M/V; // Concentration of sugar, g/cc\n", +"S = 10*theta/(l*c); // Specific rotation of sugar solution from the given data, degrees\n", +"\n", +"printf('\nThe specific rotation of sugar solution from the given data = %4.1f degrees', S);\n", +"\n", +"// Result \n", +"// The specific rotation of sugar solution from the given data = 25.1 degrees \n", +"\n", +" " + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 4.21: Angle_of_rotation_of_the_plane_of_polarization.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"// Scilab Code Ex4.21:: Page-4.35 (2009)\n", +"clc; clear;\n", +"m = 15; // Amount of sugar, g\n", +"S = 66; // Specific rotation of sugar solution from the given data, degrees\n", +"l = 20; // Length of the tube, cm\n", +"V = 100; // Volume of sugar solution, cc\n", +"c = m/V; // Concentration of sugar, g/cc\n", +"// As S = 10*theta/(l*c), solving for theta\n", +"theta = S*l*c/10; // Angle of rotation of the plane of polarization, degrees\n", +"\n", +"printf('\nThe angle of rotation of the plane of polarization = %4.1f degrees', theta);\n", +"\n", +"// Result \n", +"// The angle of rotation of the plane of polarization = 19.8 degrees \n", +"\n", +" " + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 4.22: Angle_of_rotation_of_the_optically_active_solution.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"// Scilab Code Ex4.22: : Page-4.35 (2009)\n", +"clc; clear;\n", +"l = 5; // Length of the tube, dm\n", +"m = 50; // Amount of sugar, g\n", +"S = 50; // Specific rotation of sugar solution, degrees\n", +"V = 150; // Volume of sugar solution, cc\n", +"c = m/V; // Concentration of sugar, g/cc\n", +"// As S = theta/(l*c), solving for theta\n", +"theta = S*l*c; // Angle of rotation of the optically active solution\n", +"\n", +"printf('\nThe angle of rotation of the optically active solution = %4.1f degrees', theta);\n", +"\n", +"// Result \n", +"// The angle of rotation of the optically active solution = 83.3 degrees \n", +"\n", +" " + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 4.23: Angle_of_rotation_in_a_tube_of_new_length.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"// Scilab Code Ex4.23:: Page-4.35 (2009)\n", +"clc; clear;\n", +"l = 3; // Length of the tube, dm\n", +"theta = 17.0; // Angle of rotation of the plane of polarization, degrees\n", +"c = 1.0; // For simplicity assume concentration of solution to be unity, g/cc\n", +"l_prime = 2.5; // New length of the tube, dm\n", +"c_prime = 1.25*c; // Concentration of solution with 25 cm length of tube, g/cc\n", +"theta_prime = theta*l_prime*c_prime/(l*c); // Angle of rotation in a tube of new length\n", +"\n", +"\n", +"printf('\nThe angle of rotation in a tube of new length of %3.1f cm = %4.1f degrees', l_prime, theta_prime);\n", +"\n", +"// Result \n", +"// The angle of rotation in a tube of new length of 2.5 cm = 17.7 degrees \n", +"\n", +" " + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 4.24: Mass_of_sugar_in_the_solution_contained_in_the_tube.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"// Scilab Code Ex4.24:: Page-4.36 (2009)\n", +"clc; clear;\n", +"l = 17; // Length of the tube, cm\n", +"V = 37; // Volume of sugar solution, cc\n", +"theta = 15; // Angle of rotation of the plane of polarization, degrees\n", +"S = 68; // Specific rotation of sugar solution, degrees\n", +"// As S = 10*theta/(l*c), solving for c\n", +"c = 10*theta/(l*S); // Concentration of sugar solution, g/cc\n", +"m = c*V; // Mass of sugar in the solution contained in the tube, g\n", +"\n", +"printf('\nThe mass of sugar in the solution contained in the tube = %3.1f g', m);\n", +"\n", +"// Result \n", +"// The mass of sugar in the solution contained in the tube = 4.8 g \n", +"\n", +" " + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 4.25: Percentage_purity_of_the_sugar_sample.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"// Scilab Code Ex4.25:: Page-4.36 (2009)\n", +"clc; clear;\n", +"m = 80; // Mass of sugar in the solution, g\n", +"theta = 9.9; // Angle of rotation of the plane of polarization, degrees\n", +"l = 20; // Length of the tube, cm\n", +"S_pure = 66; // Specific rotation of pure sugar solution, degrees per dm per (g/cc)\n", +"c = 0.08; // Concentration of sugar solution, g/cc\n", +"S = 10*theta/(l*c); // calculated specific rotation of sugar solution, degrees per dm per (g/cc)\n", +"percent_purity = S/S_pure*100; // Percentage purity of sugar sample, percent\n", +"\n", +"printf('\nThe percentage purity of the sugar sample = %5.2f percent', percent_purity);\n", +"\n", +"// Result \n", +"// The percentage purity of the sugar sample = 93.75 percent \n", +"\n", +" " + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 4.26: Angle_of_rotation_produced_by_the_polarimeter_plate.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"// Scilab Code Ex4.26:: Page-4.42 (2009)\n", +"clc; clear;\n", +"lambda = 6600e-010; // Wavelength of circularly polarized light, cm\n", +"mu_R = 1.53914; // Refractive index of right-handed circularly polarized light\n", +"mu_L = 1.53920; // Refractive index of left-handed circularly polarized light\n", +"t = 0.0005; // Thickness of polarimeter plate, m\n", +"theta = %pi/lambda*(mu_L-mu_R)*t; // Angle of rotation produced by the polarimeter plate, radian\n", +"\n", +"printf('\nThe angle of rotation produced by the polarimeter plate = %4.2f degrees', theta*180/%pi);\n", +"\n", +"// Result \n", +"// The angle of rotation produced by the polarimeter plate = 8.18 degrees \n", +"\n", +" " + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 4.2: Polarization_by_reflection.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"// Scilab Code Ex4.2:: Page-4.6 (2009)\n", +"clc; clear;\n", +"ip = 57; // Polarizing angle, degrees\n", +"mu = tand(ip); // Refractive index of the material from Brewster's law \n", +"printf('\nThe refractive index of the material = %4.2f', mu);\n", +"\n", +"// Result \n", +"// The refractive index of the material = 1.54 " + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 4.3: Angle_of_refraction_of_the_ray.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"// Scilab Code Ex4.3:: Page-4.6 (2009)\n", +"clc; clear;\n", +"mu = 1.53; // Refractive index of the material from Brewster's law \n", +"// As mu = tand(ip), solving for ip\n", +"ip = atand(mu); // Polarizing angle, degrees\n", +"// But mu = sind(ip)/sind(r), solving for r\n", +"r = asind(sind(ip)/mu); // Angle of refraction, degrees\n", +"\n", +"printf('\nThe angle of refraction of the ray = %4.1f degrees', r);\n", +"\n", +"// Result \n", +"// The angle of refraction of the ray = 33.2 degrees " + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 4.4: Angle_of_minimum_deviation_for_green_light.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"// Scilab Code Ex4.4:: Page-4.6 (2009)\n", +"clc; clear;\n", +"ip = 60; // Polarizing angle, degrees\n", +"A = 60; // Angle of equilateral prism, degrees\n", +"mu = tand(ip); // Refractive index of the material from Brewster's law \n", +"// For angle of minimum deviation in prism, delta_m, refractive index\n", +"// mu = sind((A+delta_m)/2)/sind(A/2), solving for delta_m\n", +"delta_m = 2*asind(mu*sind(A/2))-A; // Angle of minimum deviation, degrees\n", +"\n", +"printf('\nThe angle of minimum deviation for green light = %2d degrees', ceil(delta_m));\n", +"\n", +"// Result \n", +"// The angle of minimum deviation for green light = 60 degrees " + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 4.5: Polarizing_angles_of_the_materials_for_given_refractive_indices.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"// Scilab Code Ex4.5:: Page-4.7 (2009)\n", +"clc; clear;\n", +"mu = [1.33 1.65 1.55]; // Refractive indices of the material\n", +"// As mu = tand(ip), solving for ip\n", +"ip = atand(mu); // Brewster's law gives polarizing angle, degrees\n", +"for i =1:1:3 \n", +"printf('\nmu = %4.2f, ip = %4.1f degrees', mu(i), ip(i));\n", +"end\n", +"\n", +"// Result \n", +"// mu = 1.33, ip = 53.1 degrees\n", +"// mu = 1.65, ip = 58.8 degrees\n", +"// mu = 1.55, ip = 57.2 degrees " + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 4.6: Angle_of_rotation_of_analyser.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"// Scilab Code Ex4.6:: Page-4.8 (2009)\n", +"clc; clear;\n", +"E0 = 1; // For simplicity assume maximum intensity through polarizer and analyser to be unity, unit\n", +"E = 1/6*E0; // One-sixth of the maximum intensity, unit\n", +"// From Malus law, E = E0*cosd(theta)^2, solving for theta\n", +"theta = acosd(sqrt(E)); // Angle through which analyser should be rotated, degrees\n", +"printf('\nThe angle of rotation of analyser = %4.1f degrees', theta);\n", +"\n", +"// Result \n", +"// The angle of rotation of analyser = 65.9 " + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 4.7: Angles_of_rotation_of_analyser_for_given_transmitted_light_intensities.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"// Scilab Code Ex4.7:: Page-4.8 (2009)\n", +"clc; clear;\n", +"E0 = 1; // For simplicity assume maximum intensity through polarizer and analyser to be unity, unit\n", +"light_fraction = [0.25 0.45 0.65 0.75 0.0];\n", +"for i = 1:1:5\n", +"E = light_fraction(i)*E0; // Light fraction of the maximum intensity, unit\n", +"// From Malus law, E = E0*cosd(theta)^2, solving for theta\n", +"theta = acosd(sqrt(E)); // Angle through which analyser should be rotated, degrees\n", +"printf('\nE = %4.2fE0, theta = %4.1f degrees', light_fraction(i), theta);\n", +"end\n", +"\n", +"// Result \n", +"// E = 0.25E0, theta = 60.0 degrees\n", +"// E = 0.45E0, theta = 47.9 degrees\n", +"// E = 0.65E0, theta = 36.3 degrees\n", +"// E = 0.75E0, theta = 30.0 degrees\n", +"// E = 0.00E0, theta = 90.0 degrees " + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 4.8: Angle_of_minimum_deviation_for_green_light.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"// Scilab Code Ex4.8:: Page-4.9 (2009)\n", +"clc; clear;\n", +"ip = 60; // Polarizing angle, degrees\n", +"mu = tand(ip); // Brewster's law giving refractive index\n", +"A = 60; // Angle of prism, degrees\n", +"d = (mu - 1)*A; // Angle of minimum deviation for green light, degrees\n", +"\n", +"printf('\nThe angle of minimum deviation for green light = %5.2f degrees', d);\n", +"\n", +"// Result \n", +"// The angle of minimum deviation for green light = 43.92 degrees \n", +" " + ] + } +, +{ + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Example 4.9: Ratio_of_ordinary_to_extraordinary_ray_intensities.sce" + ] + }, + { +"cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], +"source": [ +"// Scilab Code Ex4.9:: Page-4.9 (2009)\n", +"clc; clear;\n", +"theta = 30; // Angle which the plane of vibration makes with the incident beam, degrees\n", +"// As intensity of ordinary and extraordinary ray are\n", +"// E_E = A^2*cosd(theta)^2 and E_O = A^2*sind(theta)^2, solving for E_E/E_O\n", +"EE_ratio_EO = cotd(30)^2; // Ratio of ordinary and extraordinary ray intensities \n", +"\n", +"printf('\nThe ratio of ordinary to extraordinary ray intensities = %d', EE_ratio_EO);\n", +"\n", +"// Result \n", +"// The ratio of ordinary to extraordinary ray intensities = 3 \n", +" " + ] + } +], +"metadata": { + "kernelspec": { + "display_name": "Scilab", + "language": "scilab", + "name": "scilab" + }, + "language_info": { + "file_extension": ".sce", + "help_links": [ + { + "text": "MetaKernel Magics", + "url": "https://github.com/calysto/metakernel/blob/master/metakernel/magics/README.md" + } + ], + "mimetype": "text/x-octave", + "name": "scilab", + "version": "0.7.1" + } + }, + "nbformat": 4, + "nbformat_minor": 0 +} |