Dear All, On doing the following integral I get different results analytically and numerically. # Define the integrand P3 := proc (x, m) options operator, arrow; -(1/2)*(-1+4*x^2*m^2+8*x^3*m+2*x*m-3*x^2+m^2)/((2*x^2+1+2*sqrt((x+m)^2/((x+m+1)*(x+m-1)))*x^2)*(x+m+1)*(x+m-1)*sqrt(x)) end proc; IntegralP3 := proc (m) options operator, arrow; Int(P3(x, m), x = 0 .. infinity) end proc; evalf(IntegralP3(1.0)); -2.262741700 int(P3(x,1),x=0..infinity); 0 evalf(-8*sqrt(2)/5); -2.262741699 Does anybody know why I get different results? Drew

Please Wait...