//chapter12,Example12_2,pg 357

lam=1*10^-6//wavelength

n1=1.53

n2=1.5

NA=sqrt((n1^2)-(n2^2))

a=(2.405*lam)/(2*%pi*NA)

printf("core radius\n")

printf("a=%.8f m",a)