initial commit / add all books

1 files changed, 45 insertions, 0 deletions

diff --git a/3638/CH17/EX17.1/Ex17_1.sce b/3638/CH17/EX17.1/Ex17_1.sce
new file mode 100644
index 000000000..7fa24ba5d
--- /dev/null
+++ b/3638/CH17/EX17.1/Ex17_1.sce
@@ -0,0 +1,45 @@
+//Introduction to Fiber Optics by A. Ghatak and K. Thyagarajan, Cambridge, New Delhi, 1999

+//Example 17.1

+//OS=Windows XP sp3

+//Scilab version 5.5.2

+clc;

+clear;

+//given

+n1=1.4532;//refractive index of core

+n2=1.45;//refractive index of cladding

+a=5e-6;//fiber core radius in m

+d=12e-6;//Distance between the fiber axes in m

+dbar=d/a;//Ratio of distance between fiber axes to the core radius

+delta=((n1)^2-(n2)^2)/((n1)^2);//Dimensionless quantity

+

+//Case (i):

+lambda0=1.3e-6;//Free space wavelength in m

+k0=2*%pi/lambda0;//free space wave number in rad/m

+V=k0*a*sqrt((n1^2)-(n2^2));//dimensionless waveguide parameter

+//The approximate expression for k consists of constants A, B and C

+A=5.2789-3.663*V+0.3841*(V^2);//Expression for constant A in terms of 'V'

+B=-0.7769+1.2252*V-0.0152*(V^2);//Expression for constant B in terms of 'V'

+C=-0.0175-0.0064*V-0.0009*(V^2);//Expression for constant C in terms of 'V'

+k=(%pi/(2*a))*sqrt(delta)*exp(-(A+B*dbar+C*(dbar)^2));//Expression for Coupling Coefficient in m^(-1)

+mprintf("\n For lambda=1.3 um:");

+mprintf("\n k=%f mm^(-1)",k/1e3);//Dividing by 10^3 to conevert into mm^(-1)

+//The answers vary due to round off error

+Lc=%pi/(2*k);//Corresponding coupling length in m

+mprintf("\n Lc =%.2f mm",Lc/1e-3);//Dividing by 10^(-3) to convert into mm

+P2=(sin(k*Lc/2))^2;//The coupled power at given wavelength

+mprintf("\n P2=%.2f",P2);

+

+//Case (ii):

+lambda0=1.35e-6;//Free space wavelength in m

+k0=2*%pi/lambda0;//free space wave number in rad/m

+V=k0*a*sqrt((n1^2)-(n2^2));//dimensionless waveguide parameter

+//The approximate expression for k consists of constants A, B and C

+A=5.2789-3.663*V+0.3841*(V^2);//Expression for constant A in terms of 'V'

+B=-0.7769+1.2252*V-0.0152*(V^2);//Expression for constant B in terms of 'V'

+C=-0.0175-0.0064*V-0.0009*(V^2);//Expression for constant C in terms of 'V'

+k=(%pi/(2*a))*sqrt(delta)*exp(-(A+B*dbar+C*(dbar)^2));//Expression for Coupling Coefficient in m^(-1)

+mprintf("\n For lambda=1.35 um:");

+mprintf("\n k=%f mm^(-1)",k/1e3);//Dividing by 10^3 to conevert into mm^(-1)

+//The answers vary due to round off error

+P2=(sin(k*Lc/2))^2;//The coupled power at given wavelength

+mprintf("\n P2=%.2f",P2);//The answers vary due to round off error