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.: