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