real function moller_asy(theta) real theta thr=theta*3.14159/180.0 azz = -sin(thr)**2*(8.-sin(thr)**2)/(4.-sin(thr)**2)**2 moller_asy=azz return end