clc;
clear all;
n1 = 1.55; // Refractive index of core
n2 = 1.50; // Refractice index of clad
NA = sqrt(n1^2-n2^2); // Numerical Aperture
disp('',NA,'The numerical aperture of fiber is')