Advanced Search
All Posts : MaplePrimes Commons General Technical Discussions

If I evaluate :

assume(0 <= a); full := integrate((sin(t)-t*cos(t))^2*cos(a*t)/t^5, t = 0 .. infinity)

I seem to be getting the wrong result. Not only does the expression look strange (imaginary parts for a real integrand), but it seems also to give wrong results e.g.:

>                   evalf(subs(a = 3000, full))
> print();
                       -3999.7500 + 0. I
 

while the upper bound (as obtained by setting a=0) is 1/4.

Is my observation correct (i.e. a problem with Maple 11) or did I overlook something?

Thanks,

 

Joost
infinity

Please Wait...