clc();
clear;
//To calculate the numerical aperture
n1=1.49;                    //refractive index of core
n2=1.46;                    //refractive index of cladding
//NA=sqrt((n1^2)-(n2^2))
NA=sqrt((n1^2)-(n2^2))      //Numerical aperture
printf("The numerical aperture is %f",NA);