First time here? Check out the FAQ!

Ask Your Question
3

Indefinite integral is incorrect

asked 6 years ago

proy87 gravatar image

updated 2 years ago

tmonteil gravatar image

indefinite_integral(sqrt(1+cos(x)**2), x).full_simplify() gives 1/6*sin(x)^3, which is incorrect.

Preview: (hide)

Comments

1

Where did you get the indefinite_integral function ?

tmonteil gravatar imagetmonteil ( 6 years ago )

in 8.5.b0:

sage: integral(sqrt(1+cos(x)**2), x).full_simplify()
1/6*sin(x)^3
FrédéricC gravatar imageFrédéricC ( 6 years ago )

@tmonteli, I use sage 8.1

proy87 gravatar imageproy87 ( 6 years ago )

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

dazedANDconfused gravatar imagedazedANDconfused ( 6 years ago )

1 Answer

Sort by » oldest newest most voted
3

answered 6 years ago

Emmanuel Charpentier gravatar image

updated 6 years ago

slelievre gravatar image

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.

Preview: (hide)
link

Your Answer

Please start posting anonymously - your entry will be published after you log in or create a new account.

Add Answer

Question Tools

Stats

Asked: 6 years ago

Seen: 539 times

Last updated: Oct 27 '18