I was looking at a question on another forum, and tried in Maple also, and Maple also have a problem with this integral. Here is a simple version. The problem is that int() gives different numerical answer from evalf(Int). Maple can't solve this analytically, so values have to be used for the integrand before calling int()
eq2:=(a/(a + c*z))^L*exp(-z)/sqrt(z);
int(eq2,z=0..infinity); # 177.245
evalf(Int(eq2,z=0..infinity)); # 1.7551
I think the 1.7551 is the correct value. My question is: Why did maple give wrong answer from int()? Is it analytical reason, or purely numerical?
Maple 17.02, windows 7.