From 7f60ea012dd2524dae921a2a35adbf7ef21f2bb6 Mon Sep 17 00:00:00 2001
From: prashantsinalkar
Date: Tue, 10 Oct 2017 12:27:19 +0530
Subject: initial commit / add all books

---
 3638/CH8/EX8.1/Ex8_1.jpg | Bin 0 -> 83279 bytes
 3638/CH8/EX8.1/Ex8_1.sce |  31 +++++++++++++++++++++++++++++++
 2 files changed, 31 insertions(+)
 create mode 100644 3638/CH8/EX8.1/Ex8_1.jpg
 create mode 100644 3638/CH8/EX8.1/Ex8_1.sce

(limited to '3638/CH8/EX8.1')

diff --git a/3638/CH8/EX8.1/Ex8_1.jpg b/3638/CH8/EX8.1/Ex8_1.jpg
new file mode 100644
index 000000000..e111b4f38
Binary files /dev/null and b/3638/CH8/EX8.1/Ex8_1.jpg differ
diff --git a/3638/CH8/EX8.1/Ex8_1.sce b/3638/CH8/EX8.1/Ex8_1.sce
new file mode 100644
index 000000000..571d79a35
--- /dev/null
+++ b/3638/CH8/EX8.1/Ex8_1.sce
@@ -0,0 +1,31 @@
+//Introduction to Fiber Optics by A. Ghatak and K. Thyagarajan, Cambridge, New Delhi, 1999
+//Example 8.1
+//OS=Windows XP sp3
+//Scilab version 5.5.2
+clc;
+clear;
+//given Case(1) 
+n2=1.45;//refractive index of cladding
+a=3e-6;//radius of core in m
+delta=0.0064//fractional change in refractive index
+lambda0=1.546e-6;//wavelength  in m
+n1=n2/(1-delta);//refractive index of core
+V=2*(%pi)*a*sqrt((n1^2)-(n2^2))/lambda0;//corresponding dimensionless V number
+mprintf("\n For fiber 1:");
+mprintf("\n V=%.1f at lambda0=%.3f um ",V,lambda0/1e-6);//Division by 10^(-6) to convert into um
+b=0.41616;//value of dimensionless propagation constant corresponding to V=2 as per given table
+B=sqrt((n2^2)+b*((n1^2)-(n2^2)));//corresponding value of Beta/k0
+mprintf("\n Beta/k0=%f",B);//The answers vary due to round off error
+
+//given Case(2)
+n2=1.45;//refractive index of cladding
+a=2e-6;//radius of core in m
+delta=0.010//fractional change in refractive index
+lambda0=1.288e-6;//wavelength  in m
+n1=n2/(1-delta);//refractive index of core
+V=2*(%pi)*a*sqrt((n1^2)-(n2^2))/lambda0;//corresponding dimensionless V number
+mprintf("\n For fiber 2:");
+mprintf("\n V=%.1f at lambda0=%.3f um ",V,lambda0/1e-6);//Division by 10^(-6) to convert into um
+b=0.41616;//value of dimensionless propagation constant corresponding to V=2 as per given table
+B=sqrt((n2^2)+b*((n1^2)-(n2^2)));//corresponding value of Beta/k0
+mprintf("\n Beta/k0=%f",B);//The answers vary due to round off error
-- 
cgit