acceptance angle of meridonal rays fiber formula code example

Example: acceptance angle

n0 sin Q1= underroot of (n1^2 -n2^2)= NA     , where n0 refractive index of air=1,
Q1 is acceptance angle
n1 is the refractive index of core
n2 is the refrective index of cladding