summaryrefslogtreecommitdiff
path: root/1952/CH3/EX3.1/Ex3_1.sce
blob: 0d12185f4b66641e5b14a0a20f32bda38ce9c033 (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
// chapter 3 , Example 3.1 , pg 84
n1=1.5//core refractive index
n2=1.47//cladding refractive index
n0=1//refractive index of air
NA=sqrt(n1^2-n2^2)//numerical aperture
alpha_m =asin(NA/n0)//angle of acceptance   (in radian)
phi_m=asin((n0*sin(alpha_m))/n1)//   no*sin(alpha_m)=n1*sin(phi_m) (in radian)
phi_c=asin(n2/n1)  //critical angle   (in radian)
printf("NA=%.1f \n",NA)
printf("alpha_m=%.2f  degree\n",(alpha_m*180)/%pi)
printf("phi_m=%.2f degree\n",(phi_m*180)/%pi)
printf("phi_c=%.2f  degree",(phi_c*180)/%pi)


//data given is n2=1.97  which is not possible since refractive index of cladding should always be less than refractive index of core
//in calculation n2=1.47