Drexel dragonThe Math ForumDonate to the Math Forum



Search All of the Math Forum:

Views expressed in these public forums are not endorsed by Drexel University or The Math Forum.


Math Forum » Discussions » sci.math.* » sci.math.symbolic.independent

Topic: A bug in Reduce package 'algint'?
Replies: 9   Last Post: Feb 1, 2013 6:02 PM

Advanced Search

Back to Topic List Back to Topic List Jump to Tree View Jump to Tree View   Messages: [ Previous | Next ] Topics: [ Previous | Next ]
clicliclic@freenet.de

Posts: 978
Registered: 4/26/08
Re: A bug in Reduce package 'algint'?
Posted: Jan 22, 2013 4:24 PM
  Click to see the message monospaced in plain text Plain Text   Click to reply to this topic Reply


Axel Vogt schrieb:
> On 22.01.2013 18:42, clicliclic@freenet.de wrote:
> > acer schrieb:
> >>
> >> Just because Axel's brought up Maple (and not wishing to side-track Martin), in Maple 16.02,
> >>
> >> expr := sqrt(sqrt(a^2 + x^2) + x)/x:
> >> p := u = sqrt(a^2 + x^2) + x:
> >> new := student[changevar](p, Int(expr,x), u):
> >> sol := eval(value(new),p):
> >> lprint(sol);
> >>
> >> -((2*a^2*((a^2+x^2)^(1/2)+x)^2+a^4+((a^2+x^2)^(1/2)+x)^4)

> > /((a^2+x^2)^(1/2)+x)^2)^(1/2)
> > *((a^2+x^2)^(1/2)+x)*4^(1/2)*(-((a^2+x^2)^(1/2)+x)^(1/2)
> > +a^(1/2)*arctan(((a^2+x^2)^(1/2)+x)^(1/2)/a^(1/2))
> > +a^(1/2)*arctanh(((a^2+x^2)^(1/2)+x)^(1/2)/a^(1/2)))
> > /(a^2+((a^2+x^2)^(1/2)+x)^2)

> >>
> >> other := expand(radnormal(sol),power):
> >> lprint(other):
> >>
> >> 2*((a^2+x^2)^(1/2)+x)^(1/2)-2*a^(1/2)*arctan(((a^2+x^2)^(1/2)+x)^(1/2)/a^(1/2))-2*a^(1/2)*arctanh(((a^2+x^2)^(1/2)+x)^(1/2)/a^(1/2))
> >>
> >> What should we expect from the system, automatically?
> >>

> >
> > My outside viewpoint: if the integrand is given in terms of elementary
> > functions, an elementary antiderivative whenever one exists. Even if it
> > can be expressed more compactly in terms of higher functions. But then
> > an infinite 3F2 series is not really simple.
> >
> > The Wolfram Integrator produces an elementary answer, but includes an
> > unnecessary extra factor. Perhaps this is arrived at by automatic
> > simplification of 3F2(-1/4, -1/4, 1/4; 1/2, 3/4; -a^2/x^2) ?
> >
> > How does Maple arrive at 3F2 in the first place? Solution of an ODE?
> > Series expansion of the integrand?
> >
> > Side-tracking won't harm my ramblings,
> >

>
> I have not checked, how Maple 'finds' the solution, but it does
> not have implemented the 'catalog' of rules you have worked out
> and at hand now. Unfortunately.


I suspect you may be referring to a paper discussed last autumn,
something about umpteen theorems on a few pages. As Maple appears to be
comparatively weak on pattern matching, it would probably be better here
to implement this kind of reduction by setting up and solving the
systems of linear equations anew for each problem (by the Extended
Euclidean Algorithm if you like), and generalizing to arbitrary order at
the same time.

However, those reduction rules only suffice to handle algebraic
integrands not involving nested roots - for the latter you need
Bronstein's generalization of Hermite reduction, of which Davenport's
version was a forerunner, I believe. Waldek should have aquired a much
deeper knowledge of all this - perhaps you wait for his review of
"Variations on modern Hermite reduction" before starting to work?

I am a bit curious what Reduce (with 'algint' loaded) does return for
Vladimir's infamous

INT(SQRT(SQRT(x^4 + 1) + x^2)/((x + 1)^2*SQRT(x^4 + 1)), x)

<http://mathforum.org/kb/message.jspa?messageID=6271992>

Does anybody want to give it a try?

Martin.



Point your RSS reader here for a feed of the latest messages in this topic.

[Privacy Policy] [Terms of Use]

© Drexel University 1994-2014. All Rights Reserved.
The Math Forum is a research and educational enterprise of the Drexel University School of Education.