Account Options

  1. Sign in
The old Google Groups will be going away soon, but your browser is incompatible with the new version.
Google Groups Home
« Groups Home
Message from discussion Integration algorithm
The group you are posting to is a Usenet group. Messages posted to this group will make your email address visible to anyone on the Internet.
Your reply message has not been sent.
Your post was successful
 
From:
To:
Cc:
Followup To:
Add Cc | Add Followup-to | Edit Subject
Subject:
Validation:
For verification purposes please type the characters you see in the picture below or the numbers you hear by clicking the accessibility icon. Listen and type the numbers you hear
 
Manuel Bronstein  
View profile  
 More options Nov 24 2000, 3:00 am
Newsgroups: comp.soft-sys.math.maple
From: Manuel Bronstein <manuel.bronst...@sophia.inria.fr>
Date: 2000/11/24
Subject: Re: Integration algorithm

MYWY Becker wrote:

> Does anyone know whether the integration algorithm in Maple is
> Risch-complete? That is, if Maple's answer is not an elementary
> function, can you be sure that there doesn't exist an elementary answer?
> In other words, has Maple implemented the *entire* Risch-like algorithm,
> in particular the *entire* algebraic case?

It is not Risch-complete, even in the purely transcendental
case (recursion problems in the logarithmic case). It is not
complete either in the purely algebraic case (see below),
and certainly not in the mixed algebraic-transcendental case.

Neither is Axiom (despite various rumors). The difference
is that Axiom issues an error message when an unimplemented
branch of the algorithm is hit, so an unevaluated integral
in Axiom is a proof that the integral is not elementary.

Neither is Mathematica, despite all the hype and ads
(source code is unavailable but "black box" experiments
show Mma's integrator to be in the, ahem, "low" category).

For the sceptics, here is an elementary integral of an
algebraic function, whose integral is a simple logarithm,
missed by the Risch integrators of both Maple and Mma:

Mathematica 4.0:
g = x / Sqrt[x^4 + 10 x^2 - 96 x - 71]
Integrate[g,x]
f = -Log[(x^6+15 x^4-80 x^3+27 x^2-528 x+781) Sqrt[x^4+10 x^2-96 x-71]
     - x^8 - 20 x^6 + 128 x^5 - 54 x^4 + 1408 x^3 - 3124 x^2 - 10001]/8
Simplify[D[f,x] - g]

Maple 5.5 and Maple 6:
g := x / sqrt(x^4 + 10*x^2 - 96*x - 71);
int(g,x);
int(convert(g,RootOf),x);
f := -log((x^6+15*x^4-80*x^3+27*x^2-528*x+781) *
sqrt(x^4+10*x^2-96*x-71)
     - x^8 - 20*x^6 + 128*x^5 - 54*x^4 + 1408*x^3 - 3124*x^2 - 10001)/8;
normal(diff(f,x)-g);

For fairness, here is an easy one that Axiom cannot decide:
integrate(sqrt atan x,x)

-- Manuel Bronstein
-- Manuel.Bronst...@sophia.inria.fr
-- http://www.inria.fr/cafe/Manuel.Bronstein/


 
You must Sign in before you can post messages.
To post a message you must first join this group.
Please update your nickname on the subscription settings page before posting.
You do not have the permission required to post.