initial commit / add all books

26 files changed, 175 insertions, 0 deletions

diff --git a/1271/CH5/EX5.1/1.txt b/1271/CH5/EX5.1/1.txt
new file mode 100755
index 000000000..d762266be
--- /dev/null
+++ b/1271/CH5/EX5.1/1.txt
@@ -0,0 +1,3 @@
+ Critical angle = 68.068581 degree, 

+ Numerical aperture = 0.567715. 

+ Maximum incidence angle = 34.591019 degree. 

diff --git a/1271/CH5/EX5.1/example5_1.sce b/1271/CH5/EX5.1/example5_1.sce
new file mode 100755
index 000000000..41458e675
--- /dev/null
+++ b/1271/CH5/EX5.1/example5_1.sce
@@ -0,0 +1,11 @@
+clc 
+// Given that
+mu1 = 1.52 // refractive index for core
+mu2 = 1.41 // refractive index for cladding
+// Sample Problem 1 on page no. 5.15
+printf("\n # PROBLEM 1 # \n")
+theta_c = asin(mu2 / mu1) * (180 / %pi)
+NA = sqrt(mu1^2 - mu2^2)
+theta_0 = asin(NA) * (180 / %pi)
+printf("\n Standard formula used \n theta_c = asin(mu2 / mu1) * (180 / pi). \n NA = sqrt(mu1^2 - mu2^2). \n theta_0 = asin(NA) * (180 / pi). \n")
+printf("\n Critical angle = %f degree. \n Numerical aperture = %f,\n Maximum incidence angle = %f degree.",theta_c,NA,theta_0)
diff --git a/1271/CH5/EX5.10/10.txt b/1271/CH5/EX5.10/10.txt
new file mode 100755
index 000000000..fe03bfa14
--- /dev/null
+++ b/1271/CH5/EX5.10/10.txt
@@ -0,0 +1,2 @@
+Number of modes in first case = 44.104795. 

+ Number of modes in second case = 4410.479467.

diff --git a/1271/CH5/EX5.10/example5_10.sce b/1271/CH5/EX5.10/example5_10.sce
new file mode 100755
index 000000000..03558e9a9
--- /dev/null
+++ b/1271/CH5/EX5.10/example5_10.sce
@@ -0,0 +1,11 @@
+clc 
+// Given that
+mu1 = 3.6 // refractive index for core
+mu2 = 3.55 // refractive index for cladding
+// Sample Problem 10 on page no. 5.19
+printf("\n # PROBLEM 10 # \n")
+NA = sqrt(mu1^2 - mu2^2)//calculation for numerical aperture
+Mm1 = 0.5 * (%pi * 5 * NA)^2//calculation for no. of modes in first case
+Mm2 = 0.5 * (%pi * 50 * NA)^2//calculation for no. of modes in second case
+printf("\n Standard formula used \n Mm=1/2(pi*d*NA/lambda)^2. NA=sqrt(mu1^2-mu2^2). \n")
+printf("\n Number of modes in first case = %d. \n Number of modes in second case = %d.",Mm1,Mm2)
diff --git a/1271/CH5/EX5.11/11.txt b/1271/CH5/EX5.11/11.txt
new file mode 100755
index 000000000..7aa82f219
--- /dev/null
+++ b/1271/CH5/EX5.11/11.txt
@@ -0,0 +1 @@
+Maximum diameter of core = 1.020804e-05 meter
\ No newline at end of file
diff --git a/1271/CH5/EX5.11/example5_11.sce b/1271/CH5/EX5.11/example5_11.sce
new file mode 100755
index 000000000..a3168cd32
--- /dev/null
+++ b/1271/CH5/EX5.11/example5_11.sce
@@ -0,0 +1,11 @@
+clc 
+// Given that
+lambda = 1.25e-6 //wavelength of light in meter
+mu1 = 1.46 // refractive index for core
+mu2 = 1.457 // refractive index for cladding
+// Sample Problem 11 on page no. 5.20
+printf("\n # PROBLEM 11 # \n")
+NA = sqrt(mu1^2 - mu2^2)//calculation for numerical aperture
+k = (2.4 * lambda) / ( %pi * NA)
+printf("\n Standard formula used \n d<8*lambda/(pi*NA)\n")
+printf("\n Maximum diameter of core = %f micro meter",k*1e6)
diff --git a/1271/CH5/EX5.12/12.txt b/1271/CH5/EX5.12/12.txt
new file mode 100755
index 000000000..8ef3405bf
--- /dev/null
+++ b/1271/CH5/EX5.12/12.txt
@@ -0,0 +1 @@
+ Absorption coefficient = 69.897000 dB/km. 
\ No newline at end of file
diff --git a/1271/CH5/EX5.12/example5_12.sce b/1271/CH5/EX5.12/example5_12.sce
new file mode 100755
index 000000000..c2e476418
--- /dev/null
+++ b/1271/CH5/EX5.12/example5_12.sce
@@ -0,0 +1,10 @@
+clc 
+// Given that
+L = 0.1 // length of fiber in km
+p = 5e-6 // power of signal in watt
+p_ = 1e-6 // power of signal inside the fiber in watt
+// Sample Problem 12 on page no. 5.20
+printf("\n # PROBLEM 12 # \n")
+alpha = (10 * log10(p / p_)) / L//calculation for absorption coefficient
+printf("\n Standard formula used \n alpha=10/L*log(Pi/Po).\n")
+printf("\n Absorption coefficient = %f dB/km. ",alpha)
diff --git a/1271/CH5/EX5.13/13.txt b/1271/CH5/EX5.13/13.txt
new file mode 100755
index 000000000..8ac45fa54
--- /dev/null
+++ b/1271/CH5/EX5.13/13.txt
@@ -0,0 +1 @@
+ Output power = 7.918550e-05 watt.  
\ No newline at end of file
diff --git a/1271/CH5/EX5.13/example5_13.sce b/1271/CH5/EX5.13/example5_13.sce
new file mode 100755
index 000000000..d058d9137
--- /dev/null
+++ b/1271/CH5/EX5.13/example5_13.sce
@@ -0,0 +1,11 @@
+clc 
+// Given that
+L = 3 // length of optical fiber in km
+l = 6 // losses in dB
+p = 5e-3 // input power in watt 
+// Sample Problem 13 on page no. 5.20
+printf("\n # PROBLEM 13 # \n")
+alpha = (l * 3) / L
+p_ = p / (exp((2.303 * alpha * L) / 10))
+printf("\n Standard formula used \n alpha = (l * 3) / L. \n p_ = p / (exp((2.303 * alpha * L) / 10)). \n")
+printf("\n Output power = %f mW. ",p_*1e3)
diff --git a/1271/CH5/EX5.2/2.txt b/1271/CH5/EX5.2/2.txt
new file mode 100755
index 000000000..1f62a21d2
--- /dev/null
+++ b/1271/CH5/EX5.2/2.txt
@@ -0,0 +1,2 @@
+Numerical aperture = 0.556776.

+ Maximum incidence angle = 33.833159 degree.

diff --git a/1271/CH5/EX5.2/example5_2.sce b/1271/CH5/EX5.2/example5_2.sce
new file mode 100755
index 000000000..c3eb8e9c8
--- /dev/null
+++ b/1271/CH5/EX5.2/example5_2.sce
@@ -0,0 +1,10 @@
+clc 
+// Given that
+mu1 = 1.6 // refractive index for core
+mu2 = 1.5 // refractive index for cladding
+// Sample Problem 2 on page no. 5.16
+printf("\n # PROBLEM 2 # \n")
+NA = sqrt(mu1^2 - mu2^2)//calculation for numerical aperture
+theta_0 = asin(NA) * (180 / %pi)//calculation for maximum incidence angle
+printf("Standard formula used \n NA=aqrt(mu1^2-mu2^2),\n sin(theta_)=NA. \n")
+printf("\n Numerical aperture = %f.\n Maximum incidence angle = %f degree.",NA,theta_0)
diff --git a/1271/CH5/EX5.3/3.txt b/1271/CH5/EX5.3/3.txt
new file mode 100755
index 000000000..cab1a7785
--- /dev/null
+++ b/1271/CH5/EX5.3/3.txt
@@ -0,0 +1,4 @@
+ Critical angle = 80.633211 degree. 

+ Numerical aperture = 0.244131. 

+ Acceptance angle = 14.130491 degree.

+ Fractional refractive index = 0.013333.

diff --git a/1271/CH5/EX5.3/example5_3.sce b/1271/CH5/EX5.3/example5_3.sce
new file mode 100755
index 000000000..dfebde6bf
--- /dev/null
+++ b/1271/CH5/EX5.3/example5_3.sce
@@ -0,0 +1,13 @@
+clc 
+// Given that
+mu_0 = 1 // refractive index of air
+mu1 = 1.5 // refractive index for core
+mu2 = 1.48 // refractive index for cladding
+// Sample Problem 3 on page no. 5.16
+printf("\n # PROBLEM 3 # \n")
+theta_c = asin(mu2 / mu1) * (180 / %pi)
+delta_mu = (mu1 - mu2) / mu1
+NA = sqrt(mu1^2 - mu2^2)
+theta_0 = asin(NA) * (180 / %pi)
+printf("\n Standard formula used \n theta_c = asin(mu2 / mu1) * (180 / pi). \n delta_mu = (mu1 - mu2) / mu1. \n NA = sqrt(mu1^2 - mu2^2). \n theta_0 = asin(NA) * (180 / pi). \n ")
+printf("\n Critical angle = %f degree. \n Numerical aperture = %f. \n Acceptance angle = %f degree.\n Fractional refractive index = %f.",theta_c,NA,theta_0,delta_mu)
diff --git a/1271/CH5/EX5.4/4.txt b/1271/CH5/EX5.4/4.txt
new file mode 100755
index 000000000..7d7ecd6ed
--- /dev/null
+++ b/1271/CH5/EX5.4/4.txt
@@ -0,0 +1,2 @@
+ Numerical aperture = 0.560357. 

+ Maximum incidence angle = 34.080492 degree.

diff --git a/1271/CH5/EX5.4/example5_4.sce b/1271/CH5/EX5.4/example5_4.sce
new file mode 100755
index 000000000..017688d16
--- /dev/null
+++ b/1271/CH5/EX5.4/example5_4.sce
@@ -0,0 +1,10 @@
+clc 
+// Given that
+mu1 = 1.62 // refractive index for core
+mu2 = 1.52 // refractive index for cladding
+// Sample Problem 4 on page no. 5.17
+printf("\n # PROBLEM 4 # \n")
+NA = sqrt(mu1^2 - mu2^2)
+theta_0 = asin(NA) * (180 / %pi)
+printf("\n Standard formula used \n NA = sqrt(mu1^2 - mu2^2). \n theta_0 = asin(NA) * (180 / pi). \n")
+printf("\n Numerical aperture = %f. \n Maximum incidence angle = %f degree.",NA,theta_0)
diff --git a/1271/CH5/EX5.5/5.txt b/1271/CH5/EX5.5/5.txt
new file mode 100755
index 000000000..3a690be31
--- /dev/null
+++ b/1271/CH5/EX5.5/5.txt
@@ -0,0 +1,2 @@
+ Refractive index for core = 1.424373.

+ Refractive index for cladding = 1.407281.

diff --git a/1271/CH5/EX5.5/example5_5.sce b/1271/CH5/EX5.5/example5_5.sce
new file mode 100755
index 000000000..b49aaa36b
--- /dev/null
+++ b/1271/CH5/EX5.5/example5_5.sce
@@ -0,0 +1,10 @@
+clc 
+// Given that
+NA = 0.22 // numerical aperture
+delta_mu = 0.012 // fractional refractive index
+// Sample Problem 5 on page no. 5.17
+printf("\n # PROBLEM 5 # \n")
+mu1 = sqrt(NA^2 / (1 - (1 - delta_mu)^2))
+mu2 = (1 - delta_mu) * mu1
+printf("\n Standard formula used \n mu1 = sqrt(NA^2 / (1 - (1 - delta_mu)^2)). \n mu2 = (1 - delta_mu) * mu1. \n")
+printf("\n Refractive index for core = %f.\n Refractive index for cladding = %f.",mu1,mu2)
diff --git a/1271/CH5/EX5.6/6.txt b/1271/CH5/EX5.6/6.txt
new file mode 100755
index 000000000..ba1897e93
--- /dev/null
+++ b/1271/CH5/EX5.6/6.txt
@@ -0,0 +1,4 @@
+ Numerical aperture = 0.639375. 

+ Acceptance angle = 39.745208 degree. 

+ Number of reflections at maximum incidence = 6468.493987. 

+ Number of reflections in second case = 3204.446178. 

diff --git a/1271/CH5/EX5.6/example5_6.sce b/1271/CH5/EX5.6/example5_6.sce
new file mode 100755
index 000000000..8fdf61106
--- /dev/null
+++ b/1271/CH5/EX5.6/example5_6.sce
@@ -0,0 +1,15 @@
+clc 
+// Given that
+d = 0.0064 // diameter of fiber in cm
+mu1 = 1.53 // refractive index for core
+mu2 = 1.39 // refractive index for clad
+L = 90 // length of fiber in cm
+mu_0 = 1 // refractive index of air
+// Sample Problem 6 on page no. 5.17
+printf("\n # PROBLEM 6 # \n")
+NA = sqrt(mu1^2 - mu2^2)
+theta_0 = asin(NA) * (180 / %pi)
+N1 = L / (d * sqrt((mu1 / (mu_0 * sin(theta_0 * (%pi / 180))))^2 - 1))
+N2 = L / (d * sqrt((mu1 / (mu_0 * sin(theta_0 * (%pi / 360))))^2 - 1))
+printf("\n Standard formula used \n NA = sqrt(mu1^2 - mu2^2). \n theta_0 = asin(NA) * (180 / pi). \n N = L / (d * sqrt((mu / (mu_0 * sin(theta_0 * (pi / 180))))^2 - 1)). \n ")
+printf("\n Numerical aperture = %f.\n Acceptance angle = %f degree. \n Number of reflections at maximum incidence = %f. \n Number of reflections in second case = %f. ",NA,theta_0,N1,N2)
diff --git a/1271/CH5/EX5.7/7.txt b/1271/CH5/EX5.7/7.txt
new file mode 100755
index 000000000..f855992b0
--- /dev/null
+++ b/1271/CH5/EX5.7/7.txt
@@ -0,0 +1,2 @@
+ The normalised frequency = 40.655905,

+ number of guided in the core = 826.451303

diff --git a/1271/CH5/EX5.7/example5_7.sce b/1271/CH5/EX5.7/example5_7.sce
new file mode 100755
index 000000000..697820c95
--- /dev/null
+++ b/1271/CH5/EX5.7/example5_7.sce
@@ -0,0 +1,11 @@
+clc 
+// Given that
+d = 0.05 // diameter of fiber in mm
+NA = 0.22 // numerical aperture
+lambda = 8.5e-4 // wavelength of light in mm
+// Sample Problem 7 on page no. 5.18
+printf("\n # PROBLEM 7 # \n")
+Vn = (%pi * d * NA) / lambda
+Mm = 0.5 * (Vn)^2
+printf("\n Standard formula used \n Vn = (pi * d * NA) / lambda. \n Mm = 0.5 * (Vn)^2. \n")
+printf("\n The normalized frequency = %f,\n number of guided in the core = %f",Vn,Mm)
diff --git a/1271/CH5/EX5.8/8.txt b/1271/CH5/EX5.8/8.txt
new file mode 100755
index 000000000..824f34c04
--- /dev/null
+++ b/1271/CH5/EX5.8/8.txt
@@ -0,0 +1,2 @@
+ Diameter of core < 7.896296e-06 meter,

+ number of modes = 115

diff --git a/1271/CH5/EX5.8/example5_8.sce b/1271/CH5/EX5.8/example5_8.sce
new file mode 100755
index 000000000..c551388b3
--- /dev/null
+++ b/1271/CH5/EX5.8/example5_8.sce
@@ -0,0 +1,12 @@
+clc 
+// Given that
+lambda = 1.25e-6 //wavelength of light in meter
+mu1 = 1.465 // refractive index for core
+mu2 = 1.460 // refractive index for cladding
+// Sample Problem 8 on page no. 5.18
+printf("\n # PROBLEM 8 # \n")
+NA = sqrt(mu1^2 - mu2^2)
+k = (2.4 * lambda) / ( %pi * NA)
+Mm = 0.5 * ((%pi * 50e-6 * NA) / lambda)^2
+printf("\n Standard formula used \n NA = sqrt(mu1^2 - mu2^2). \n k = (2.4 * lambda) / ( pi * NA). \n Mm = 0.5 * ((pi * 50e-6 * NA) / lambda)^2. \n ")
+printf("\n Diameter of core < %e meter,\n number of modes = %d",k,Mm)
diff --git a/1271/CH5/EX5.9/9.txt b/1271/CH5/EX5.9/9.txt
new file mode 100755
index 000000000..a7431e6cb
--- /dev/null
+++ b/1271/CH5/EX5.9/9.txt
@@ -0,0 +1,2 @@
+ Numerical aperture = 0.120768.

+ Number of modes = 159. 

diff --git a/1271/CH5/EX5.9/example5_9.sce b/1271/CH5/EX5.9/example5_9.sce
new file mode 100755
index 000000000..2bddb57ff
--- /dev/null
+++ b/1271/CH5/EX5.9/example5_9.sce
@@ -0,0 +1,12 @@
+clc 
+// Given that
+lambda = 0.85e-6 //wavelength of light in meter
+mu1 = 1.461 // refractive index for core
+mu2 = 1.456 // refractive index for clad
+d = 4e-5 // diameter of core in meter
+// Sample Problem 9 on page no. 5.19
+printf("\n # PROBLEM 9 # \n")
+NA = sqrt(mu1^2 - mu2^2)
+Mm = 0.5 * ((%pi * d * NA) / lambda)^2
+printf("\n Standard formula used \n NA = sqrt(mu1^2 - mu2^2). \n Mm = 0.5 * ((pi * d * NA) / lambda)^2. \n ")
+printf("\n Numerical aperture = %f.\n Number of modes = %d. ",NA,Mm)