1 files changed, 26 insertions, 0 deletions
diff --git a/3638/CH17/EX17.5/Ex17_5.sce b/3638/CH17/EX17.5/Ex17_5.sce
new file mode 100644
index 000000000..682e5a5c8
--- /dev/null
+++ b/3638/CH17/EX17.5/Ex17_5.sce
@@ -0,0 +1,26 @@
+//Introduction to Fiber Optics by A. Ghatak and K. Thyagarajan, Cambridge, New Delhi, 1999
+//Example 17.5
+//OS=Windows XP sp3
+//Scilab version 5.5.2
+clc;
+clear;
+//given
+//Since the reflectivity of fiber is 90%,
+R=0.9;//Reflection coefficient of fiber
+L=25e-3;//Length of fiber in m
+lambdaB=800e-9;//Bragg wavelength in m
+neff=1.4517;//Corresponding value of effective index in LP01 mode
+I=0.5;//Transverse overlap integral of modal distribution
+
+//Now, (tanh(k*L))^2=R
+//Rearranging terms, we get:
+k=atanh(sqrt(R))/L;//Corresponding coupling coefficient in m^(-1)
+mprintf("\n k=%.3f mm^(-1)",k/1e3);//Dividing by 10^3 to convert into mm^(-1)
+
+//Rearranging terms of expression k=%pi*Deltan*I/lambdaB
+Deltan=k*lambdaB/(%pi*I);//Change in refractive index
+mprintf("\n Deltan=%.2e",Deltan);//Unitless quantity
+//The answers vary due to round off error
+
+DeltaLambda=lambdaB^2/(%pi*neff*L)*sqrt((k*L)^2+(%pi)^2);//Corresponding bandwidth in m
+mprintf("\n DeltaLambda=%.2f nm",DeltaLambda/1e-9);//Division by 10^(-9) to convert into nm