//Example 16.7//

ncladding=1.460;//index of refraction for cladding
ncore=1.470;// index of refraction for glass-fiber core
thethac=asind(ncladding/ncore)
mprintf("The value of ncore taken while calculating is ncore=1.479 but in the question the value of ncore is given n=1.470")
mprintf("\nthethac = %f degree ",thethac)