Search All of the Math Forum:

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

Notice: We are no longer accepting new posts, but the forums will continue to be readable.

Topic: An independent integration test suite
Replies: 128   Last Post: Dec 8, 2013 3:21 PM

 Messages: [ Previous | Next ]
 clicliclic@freenet.de Posts: 1,245 Registered: 4/26/08
Re: The A. F. Timofeev symbolic integration test suite
Posted: Nov 16, 2013 7:47 AM

Albert Rich schrieb:
>
> On Friday, November 15, 2013 3:43:59 AM UTC-10, clicl...@freenet.de wrote:
>

> >> Ok, for integrands of the form (c+d x)^n/(a+b x) when n is symbolic,
> >> the next version of Rubi will return
> >>
> >> (c+d*x)^n/(b*n*(b*(c+d*x)/(d*(a+b*x)))^n)*
> >> 2F1(-n,-n,1-n,-(b*c-a*d)/(d*(a+b*x)))
> >>
> > > The simpler equivalent rule derived using Euler's transformation is
> >> not used since it is harder to simplify its derivative back to the
> >> original integrand.

>
> > But this is just an arbitrary property of the differentiator, right?
> > Another differentiator may give the result you would now obtain by
> > applying Euler's transformation first (and undoing it on non-elementary
> > hypergeometrics that remain in the derivative).
> >
> > So your reason is no good reason; you are just bending to the dictate of
> > WRI. I suggest that Rubi redefines 2F1 differentiation instead. The
> > optimality of Rubi's antiderivatives determines the rules to which WRI
> > must bend!
> >
> > I have spoken.
> >

>
> Ok, the Messiah has spoken. I modified the optimal antiderivatives for Timofeev Chapter 8 examples 6a.n, 6b.n and 14 in the test-suite to reflect use of Euler's transformation and posted the revised pdf file at
>
> http://www.apmaths.uwo.ca/~arich/TimofeevChapter8TestResults.pdf
>
> However, full credit is still given to those integrators that return
> valid, but suboptimal, antiderivatives for these problems.
>
> Note that the revised pdf file also includes Chapter 8 test-suite
> results for various integrators, including the forthcoming version 4.3
> of Rubi...
>

The evaluation of Example 14, INT((a + b*#e^(n*x))^(r/s), x), remains
less simple than the original version, but not significantly so. The
original evaluation involved 2F1 on the branch cut for b/a*#e^(n*x) > 0,
the new one does for -1 < b/a*#e^(n*x) < 0. How about selecting between
them according to SIGN(b/a)?

Would the evaluations of Examples 3.n, 5a.n, 5b.n and 17 perhaps also
profit from Euler's transformation?

I have rummaged my vaults and dug up a mildewed sheaf of papers with the
gospel on 2F1 differentiation and integration. The integration part
reads as follows (typed in without checking):

INT(F21(a,b,c,x), x)
= (c-1)/((a-1)*(b-1))*F21(a-1,b-1,c-1,x)
[a/=1, b/=1, c/=1]

INT(x^(b-2)*F21(a,b,c,x), x)
= 1/(b-1)*x^(b-1)*F21(a,b-1,c,x)
[b/=1]

INT(x^(c-1)*F21(a,b,c,x), x)
= 1/c*x^c*F21(a,b,c+1,x)

INT((1-x)^(b-2)*F21(a,b,c,x), x)
= (c-1)/((a-c+1)*(b-1))*(1-x)^(b-1)*F21(a,b-1,c-1,x)
[b/=1, c-a/=1, c/=1]

INT(x^(c-1)*(1-x)^(b-c-1)*F21(a,b,c,x), x)
= 1/c*x^c*(1-x)^(b-c)*F21(a+1,b,c+1,x)

This set should be augmented by applying Euler's transformation on both
sides of each formula. Inasmuch as the 2F1 integration rules are unknown
to Rubi, I suggest to implement the complete set. Older versions of
Mathematica didn't know all of them, maybe the current version still
doesn't?

A 2F1 differentiator has to apply the rules backwards. For simple
antiderivatives, an expression has to be conditioned such that the most
appropriate rule matches, not just the first one from the above list.

Looking forward to Rubi4.3forte,

Martin.

Date Subject Author
2/24/13 clicliclic@freenet.de
3/19/13 clicliclic@freenet.de
3/21/13 Waldek Hebisch
3/22/13 clicliclic@freenet.de
3/26/13 Waldek Hebisch
3/26/13 clicliclic@freenet.de
4/20/13 clicliclic@freenet.de
4/20/13 Nasser Abbasi
4/20/13 Rouben Rostamian
4/20/13 clicliclic@freenet.de
4/20/13 Rouben Rostamian
4/20/13 Axel Vogt
4/20/13 clicliclic@freenet.de
4/20/13 Axel Vogt
4/21/13 Axel Vogt
4/21/13 clicliclic@freenet.de
4/21/13 Waldek Hebisch
4/22/13 clicliclic@freenet.de
4/22/13 Axel Vogt
4/22/13 clicliclic@freenet.de
4/23/13 Waldek Hebisch
4/24/13 clicliclic@freenet.de
4/25/13 Waldek Hebisch
4/26/13 clicliclic@freenet.de
4/27/13 Waldek Hebisch
4/24/13 Richard Fateman
4/24/13 clicliclic@freenet.de
4/25/13 Richard Fateman
4/26/13 clicliclic@freenet.de
4/26/13 Axel Vogt
4/27/13 clicliclic@freenet.de
4/25/13 Waldek Hebisch
4/25/13 Peter Pein
4/25/13 Nasser Abbasi
4/26/13 Peter Pein
4/26/13 clicliclic@freenet.de
4/26/13 Peter Pein
4/26/13 clicliclic@freenet.de
4/26/13 Richard Fateman
4/27/13 clicliclic@freenet.de
4/27/13 Richard Fateman
6/30/13 clicliclic@freenet.de
6/30/13 Axel Vogt
7/1/13 clicliclic@freenet.de
7/1/13 Axel Vogt
7/1/13 Waldek Hebisch
7/2/13 clicliclic@freenet.de
7/2/13 clicliclic@freenet.de
7/2/13 clicliclic@freenet.de
7/2/13 Nasser Abbasi
7/2/13 Nasser Abbasi
7/4/13 clicliclic@freenet.de
7/4/13 Nasser Abbasi
7/4/13 Nasser Abbasi
7/5/13 clicliclic@freenet.de
7/5/13 Nasser Abbasi
7/9/13 clicliclic@freenet.de
7/10/13 Nasser Abbasi
7/10/13 Richard Fateman
7/10/13 Nasser Abbasi
7/10/13 clicliclic@freenet.de
8/6/13 clicliclic@freenet.de
9/15/13 Albert D. Rich
9/15/13 clicliclic@freenet.de
9/15/13 clicliclic@freenet.de
9/21/13 Albert D. Rich
9/21/13 clicliclic@freenet.de
9/22/13 daly@axiom-developer.org
9/24/13 daly@axiom-developer.org
9/30/13 daly@axiom-developer.org
9/22/13 Albert D. Rich
9/25/13 Albert D. Rich
9/25/13 Albert D. Rich
9/25/13 clicliclic@freenet.de
9/25/13 Albert D. Rich
9/26/13 Albert D. Rich
9/26/13 clicliclic@freenet.de
9/26/13 Albert D. Rich
9/29/13 clicliclic@freenet.de
10/1/13 Albert D. Rich
10/1/13 clicliclic@freenet.de
10/1/13 Albert D. Rich
10/5/13 clicliclic@freenet.de
10/5/13 Albert D. Rich
10/6/13 clicliclic@freenet.de
10/10/13 Albert D. Rich
10/10/13 Nasser Abbasi
10/11/13 clicliclic@freenet.de
11/6/13 Albert D. Rich
11/6/13 Nasser Abbasi
11/7/13 did
11/7/13 clicliclic@freenet.de
11/7/13 clicliclic@freenet.de
11/7/13 Albert D. Rich
11/12/13 clicliclic@freenet.de
11/12/13 Albert D. Rich
11/13/13 clicliclic@freenet.de
11/13/13 Albert D. Rich
11/14/13 clicliclic@freenet.de
11/14/13 Albert D. Rich
11/15/13 clicliclic@freenet.de
11/15/13 Albert D. Rich
11/16/13 clicliclic@freenet.de
11/16/13 clicliclic@freenet.de
11/21/13 Albert D. Rich
11/21/13 clicliclic@freenet.de
11/21/13 Nasser Abbasi
11/21/13 Albert D. Rich
11/21/13 Albert D. Rich
11/22/13 clicliclic@freenet.de
11/14/13 Albert D. Rich
11/15/13 clicliclic@freenet.de
11/15/13 Nasser Abbasi
11/16/13 clicliclic@freenet.de
11/16/13 Nasser Abbasi
11/7/13 did
11/7/13 clicliclic@freenet.de
4/20/13 Richard Fateman
4/21/13 clicliclic@freenet.de
4/20/13 Axel Vogt
4/20/13 clicliclic@freenet.de
4/20/13 Waldek Hebisch
4/21/13 G. A. Edgar
12/8/13 clicliclic@freenet.de
10/5/13 Albert D. Rich
10/6/13 clicliclic@freenet.de