The correct integral is exp(-x/L)/L*r/x**2*sin(theta)