//To calculate the numerical aperture
n1 = 1.55;    //refractive index of core
n2 = 1.50;    //refractive index of cladding
NA = sqrt(n1^2 - n2^2);
printf("numerical aperture is %5.3f",NA);