Indefinite integral is incorrect
indefinite_integral(sqrt(1+cos(x)**2), x).full_simplify()
gives 1/6*sin(x)^3
, which is incorrect.
indefinite_integral(sqrt(1+cos(x)**2), x).full_simplify()
gives 1/6*sin(x)^3
, which is incorrect.
This is a bug. Furthermore, it is not a maxima
bug, as it is often the case.
Here, Sage truly screws things up itself:
Maxima doesn't give a false answer:
sage: maxima.integrate(sqrt(1+cos(x)^2),x).sage()
integrate(sqrt(cos(x)^2 + 1), x)
When one tries to "ease" the problem,
maxima doesn't recognize the "obvious", but does not give a false answer:
sage: maxima.integrate(sqrt(1-m*sin(x)^2),x).sage()
integrate(sqrt(-m*sin(x)^2 + 1), x)
Sage does:
sage: integrate(sqrt(1-m*sin(x)^2),x)
1/4*m*x - 1/8*m*sin(2*x)
BTW: what is expected:
sage: elliptic_e(x,1/2).diff(x)
sqrt(-1/2*sin(x)^2 + 1)
One can easily check that sympy
, giac
and fricas
all fail to integrate,
but that none of them gives misleading "answers".
This one does not seem to be related to existing indefinite integral bugs, and is an original, genuine, Sage-specific one. Reported as Trac #26563.
Asked: 2018-10-26 07:08:37 -0600
Seen: 174 times
Last updated: Oct 27 '18
Two ways of integrating x↦xⁿsin(x) give contradictory results. Bug?
Why is Sage calling 1 a variable?
How to load optional gap packages
solve(x) does not fully isolate x. Can Sympy use assume()tions?
Bug in cos function in version 7.2
Error computing curvature of graph submanifold
sage cannot get the simplest form on the limit expression with binomial and exponential operations
Where did you get the
indefinite_integral
function ?in 8.5.b0:
@tmonteli, I use sage 8.1
Using integral(sqrt(1+cos(x)**2), x).full_simplify() as suggested still results in an answer which is incorrect, doesn't it? This is a nonelementary integral. See here or here