since Maxima can integrate (x^n* exp(%i*(a+b*x)), x) it seems to me it should be able to handle x^n*cos(a+b*x) which is the real part of that. And the answer is the realpart of the answer to the exp(_) expression. (even declare([x,a,b,n],real), assume(x>0)) is no help. RJF