#LAPLACE
#USEFLOAT=FALSE
#SIMPLIFY=TRUE
#ERRBND=15
#EXPDSP=5
#PREC=8
#F1=f(x)=
#F2=g(x)=
#F3=h(x)=f(x)/g(x)
#F4=matrix(
#F5=vector(
#F6=diff(f)
#F7=diff(g)
#F8=diff(h)
#F9=taylor(f,0)
#F10=taylor(f,1)
#F11=%(run type s:startup-sequence)
#F12=
#F13=
#F14=
#F15=
#F16=
#F17=
#F18=
#F19=
#F20=
const(epsilon,n)
alpha2(alpha1)=asin(sin(epsilon)*sqrt(n^2-sin(alpha1)^2) + cos(epsilon)*sin(alpha1))
delta(alpha1)=alpha1 + alpha2(alpha1) - epsilon