1 files changed, 23 insertions, 0 deletions
diff --git a/1871/CH3/EX3.13/Ch03Ex13.sce b/1871/CH3/EX3.13/Ch03Ex13.sce
new file mode 100755
index 000000000..14ea23426
--- /dev/null
+++ b/1871/CH3/EX3.13/Ch03Ex13.sce
@@ -0,0 +1,23 @@
+// Scilab code Ex3.13 : Pg:115 (2008)
+clc;clear;
+P = 5;    // Power of combination of a convex lens and a plano-convex lens, dioptre
+mu1 = 1.50;    // Refractive index of crown glass
+mu2 = 1.60;    // Refractive index of flint glass
+omega1 = 0.01;    // Dispersive power of the crown glass convex lens
+omega2 = 0.02;    // Dispersive power of flint glass plano-convex lens
+F = 100/P;    // Focal length of the combination of lenses, cm
+f_ratio = -omega2/omega1;    // Ratio of f2 to f1
+// From thin lens formula, 1/F = 1/f1 + 1/f2 and as f2 = f_ratio*f1, solving for f1
+f1 = -F/f_ratio;    // Focal length of flint glass lens, cm
+f2 = f_ratio*f1;    // Focal length of crown glass lens, cm
+mu = 1.60;    // Refractive index of flint glass
+R2 = %inf;    // Radius of curvature of the first surface of convex lens, cm
+// Now from lens maker's formula
+R1 = (mu - 1)/(1/f2+(mu-1)/R2);    // Radius of curvature of second surface of convex lens, cm
+mu = 1.5;    // Refractive index of crown glass
+R2 = R1;    // Radius of curvature of the first surface of convex lens, cm
+R1_prime = (mu - 1)/(1/f1+(mu-1)/R2);    // Radius of curvature of second surface of concave lens, cm
+printf("\nThe radii of curvature of the convex lens are = %-3.1f cm and %2.0f cm", R1_prime, R1);
+
+// Result 
+// The radii of curvature of the convex lens are = 8.6 cm and -12 cm