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: The Charlwood Fifty
Replies: 52   Last Post: Jun 24, 2013 10:24 PM

Advanced Search

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

Posts: 983
Registered: 4/26/08
Re: The Charlwood Fifty
Posted: Jun 8, 2013 1:17 PM
  Click to see the message monospaced in plain text Plain Text   Click to reply to this topic Reply


Albert Rich schrieb:
>
> I just posted a revised pdf file of the Charlwood Fifty integration
> test-suite at
>
> http://www.apmaths.uwo.ca/~arich/CharlwoodIntegrationProblems.pdf
>
> [...] Also posted is a Mathematica package file of the test-suite in
> machine readable form at
>
> http://www.apmaths.uwo.ca/~arich/CharlwoodProblems.m
>
> [...] Hopefully all the antiderivatives in Charlwood Fifty test-suite
> are now optimal...
>


I haven't checked the evaluations in the file systematically, but a new
look has revealed further possibilities for improvement.

The present solutions of problems #21 and #22 from Charlwood's appendix,
INT(x^3*ASIN(x)/SQRT(1-x^4), x) and INT(x^3*ASEC(x)/SQRT(x^4-1), x), can
be written as:

1/4*(x*SQRT(1-x^4)/SQRT(1-x^2)
+ LN(1-x^2) - LN(-x + x^3 + SQRT(1-x^2)*SQRT(1-x^4)))
- 1/2*SQRT(1-x^4)*ASIN(x)

1/2*(SQRT(x^4-1)*ASEC(x) - SQRT(x^4-1)/(x*SQRT(1 - 1/x^2))
- LN(x - x^3) + LN(1 - x^2 - x*SQRT(x^4-1)*SQRT(1 - 1/x^2)))

In my eyes, ATANH constitutes a more natural option here than LN:

1/4*(x*SQRT(1-x^4)/SQRT(1-x^2) + ATANH(x*SQRT(1-x^2)/SQRT(1-x^4)))
- 1/2*SQRT(1-x^4)*ASIN(x)

1/2*(SQRT(x^4-1)*ASEC(x) - SQRT(x^4-1)/(x*SQRT(1 - 1/x^2))
+ ATANH(x*SQRT(1 - 1/x^2)/SQRT(x^4-1)))

Or, using piecewise-constant prefactors:

SQRT(1-x^4)/(4*SQRT(1-x^2)*SQRT(1+x^2))
*(x*SQRT(1+x^2) + ATANH(x/SQRT(1+x^2))) - 1/2*SQRT(1-x^4)*ASIN(x)

1/2*(SQRT(x^4-1)*ASEC(x) - x*SQRT(x^2+1)/SQRT(x^4-1)*SQRT(1 - 1/x^2)
*(SQRT(x^2+1) - ATANH(1/SQRT(x^2+1))))

A much simpler evaluation of integral #48 from the appendix is:

INT(ATAN(SQRT(x) - SQRT(x+1)), x) =
= (x+1)*ATAN(SQRT(x) - SQRT(x+1)) - SQRT(x)/2

Finally, in the antiderivative #50, the rationalization of radicand
denominators in the numerical prefactors may be considered:

INT(ATAN(x*SQRT(1-x^2)), x) =
= x*ATAN(x*SQRT(1-x^2))
- SQRT(2*SQRT(5) + 2)/2*ATAN(SQRT(2*SQRT(5) + 2)/2*SQRT(1-x^2))
+ SQRT(2*SQRT(5) - 2)/2*ATANH(SQRT(2*SQRT(5) - 2)/2*SQRT(1-x^2))

Martin.


Date Subject Author
5/23/13
Read The Charlwood Fifty
Albert D. Rich
5/23/13
Read Re: The Charlwood Fifty
Nasser Abbasi
5/23/13
Read Re: The Charlwood Fifty
clicliclic@freenet.de
6/16/13
Read Re: The Charlwood Fifty
clicliclic@freenet.de
6/16/13
Read Re: The Charlwood Fifty
Nasser Abbasi
6/17/13
Read Re: The Charlwood Fifty
clicliclic@freenet.de
6/17/13
Read Re: The Charlwood Fifty
Nasser Abbasi
6/17/13
Read Re: The Charlwood Fifty
clicliclic@freenet.de
6/17/13
Read Re: The Charlwood Fifty
Nasser Abbasi
6/17/13
Read Re: The Charlwood Fifty
Waldek Hebisch
6/18/13
Read Re: The Charlwood Fifty
Albert D. Rich
6/24/13
Read Re: The Charlwood Fifty
clicliclic@freenet.de
6/24/13
Read Re: The Charlwood Fifty
Albert D. Rich
5/23/13
Read Re: The Charlwood Fifty
Albert D. Rich
5/23/13
Read Re: The Charlwood Fifty
Waldek Hebisch
5/24/13
Read Re: The Charlwood Fifty
Andreas Dieckmann
5/25/13
Read Re: The Charlwood Fifty
clicliclic@freenet.de
6/1/13
Read Re: The Charlwood Fifty
clicliclic@freenet.de
6/3/13
Read Re: The Charlwood Fifty
Albert D. Rich
6/6/13
Read Re: The Charlwood Fifty
clicliclic@freenet.de
6/6/13
Read Re: The Charlwood Fifty
Nasser Abbasi
6/6/13
Read Re: The Charlwood Fifty
clicliclic@freenet.de
6/8/13
Read Re: The Charlwood Fifty
clicliclic@freenet.de
6/8/13
Read Re: The Charlwood Fifty
Albert D. Rich
6/8/13
Read Re: The Charlwood Fifty
clicliclic@freenet.de
6/8/13
Read Re: The Charlwood Fifty
clicliclic@freenet.de
6/9/13
Read Re: The Charlwood Fifty
Albert D. Rich
6/9/13
Read Re: The Charlwood Fifty
Nasser Abbasi
6/9/13
Read Re: The Charlwood Fifty
Albert D. Rich
6/9/13
Read Re: The Charlwood Fifty
Nasser Abbasi
6/10/13
Read Re: The Charlwood Fifty
clicliclic@freenet.de
6/10/13
Read Re: The Charlwood Fifty
Nasser Abbasi
6/10/13
Read Re: The Charlwood Fifty
Nasser Abbasi
6/12/13
Read Re: The Charlwood Fifty
Nasser Abbasi
6/13/13
Read Re: The Charlwood Fifty
clicliclic@freenet.de
6/13/13
Read Re: The Charlwood Fifty
Nasser Abbasi
6/13/13
Read Re: The Charlwood Fifty
Waldek Hebisch
6/13/13
Read Re: The Charlwood Fifty
Nasser Abbasi
6/13/13
Read Re: The Charlwood Fifty
Waldek Hebisch
6/14/13
Read Re: The Charlwood Fifty
clicliclic@freenet.de
6/14/13
Read Re: The Charlwood Fifty
Waldek Hebisch
6/14/13
Read Re: The Charlwood Fifty / Sage and sympy
Richard Fateman
6/14/13
Read Re: The Charlwood Fifty
clicliclic@freenet.de
6/9/13
Read Re: The Charlwood Fifty and Macsyma
Richard Fateman
6/9/13
Read Re: The Charlwood Fifty and Macsyma
Albert D. Rich
6/9/13
Read Re: The Charlwood Fifty, another Macsyma result
Richard Fateman
6/9/13
Read Re: The Charlwood Fifty, another Macsyma result
Albert D. Rich
6/10/13
Read Re: The Charlwood Fifty, another Macsyma result
Richard Fateman
6/10/13
Read Re: The Charlwood Fifty, another Macsyma result
clicliclic@freenet.de
6/10/13
Read Re: The Charlwood Fifty, another Macsyma result
Waldek Hebisch
6/11/13
Read Re: The Charlwood Fifty, another Macsyma result
Nasser Abbasi
6/11/13
Read Re: The Charlwood Fifty, another Macsyma result
clicliclic@freenet.de
6/11/13
Read Re: The Charlwood Fifty, another Macsyma result
Nasser Abbasi

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.