#LAPLACE 1
#->
const(epsilon,n)
#<-
#->
alpha_2(alpha_1)=asin(sin(epsilon)*sqrt(n^2-sin(alpha_1)^2) + cos(epsilon)*sin(alpha_1))
#<-
#->
delta(alpha_1)=alpha_1 + alpha_2(alpha_1) - epsilon
#<-
#->

#<-