//example 4.4, page 92
clc
n=1.6
r1=.080//in cm
r2=-0.080
P=(n-1)*((1/r1)-(1/r2))
printf("\n The power  is %f D",P)