wrong answer for integral

asked 2015-04-13 18:46:39 +0100

beeson gravatar image

updated 2023-01-09 23:59:37 +0100

tmonteil gravatar image

Sage 6.5:

sage: (A,t) = var('A,t')
sage: assume(A > 0)
sage: assume(A < 1/3)
sage: num = A^2 * cos(t)^2 + 2*abs(A* cos(2*t))
sage: denom = 2 + A^2*cos(t)^2 + 2*abs(A*cos(2*t))
sage: beta = -num/denom
sage: integral(beta,t,0,2*pi)
-2*pi

This can't be right as beta has absolute value less than 1 so its integral over 0 to 2 pi has to have absolute value less than 2 pi.

edit retag flag offensive close merge delete

Comments

Integrals involving abs are difficult for the Maxima subsystem of Sage. Mostly they won't work.

rws gravatar imagerws ( 2015-04-14 07:52:10 +0100 )edit