Integral(exp(-R*alpha)*LegendreP[l, (-R**2 + t**2 + 1)/(2*t)], R)