// Example 2.17.d: Numerical Aperture clc; clear; close; n1=1.5;//Waveguide Refractive Index d= 0.0005;// Cange in core-cladding refractive index NA=n1*sqrt(2*d); disp(NA,"numerical aperture is")