home bbs files messages ]

Forums before death by AOL, social media and spammers... "We can't have nice things"

   sci.math.symbolic      Symbolic algebra discussion      10,432 messages   

[   << oldest   |   < older   |   list   |   newer >   |   newest >>   ]

   Message 9,144 of 10,432   
   clicliclic@freenet.de to Albert Rich   
   Re: Rubi 4.9.2 do not integrate Timofeev   
   21 Aug 16 18:12:49   
   
   Albert Rich schrieb:   
   >   
   > On Friday, August 19, 2016 at 8:42:33 PM UTC-10, clicl...@freenet.de wrote:   
   >   
   > > Hah: There are still cube-root Goursat integrals FriCAS cannot do!   
   > > And wow: a monic quadratic in ?^3 already too complicated! My   
   > > Goursat integrability check, however, immediately certifies the   
   > > whole integral as elementary:   
   > >   
   > >   goursat5a((a + b*x)/(3 + x^2), x, 1, 0, -1, 0)   
   > >   
   > >   [false, false, false, true, true]   
   > >   
   > > the final true signifying integrability. And its evaluation does not   
   > > naturally decompose into pure a and b parts:   
   > >   
   > > INT((a + b*x)/((3 + x^2)*(1 - x^2)^(1/3)), x)   
   > >  = 2^(1/3)*(a + 3*b)/24*(LN((1 - x)^3 + 2*(1 - x^2))   
   > >  - 3*LN((1 - x) + 2^(1/3)*(1 - x^2)^(1/3))   
   > >  + 2*SQRT(3)*ATAN(1/SQRT(3)*(1 - 2^(2/3)*((1 - x)/(1 - x^2)^(1/3)))))   
   > >  - 2^(1/3)*(a - 3*b)/24*(LN((1 + x)^3 + 2*(1 - x^2))   
   > >  - 3*LN((1 + x) + 2^(1/3)*(1 - x^2)^(1/3))   
   > >  + 2*SQRT(3)*ATAN(1/SQRT(3)*(1 - 2^(2/3)*((1 + x)/(1 - x^2)^(1/3)))))   
   > >   
   > > Alternative formulations are possible, but none has one part pure a   
   > > and the other pure b.   
   > > [...]   
   >   
   > Martin's antiderivative for (a + b*x)/((3 + x^2)*(1 - x^2)^(1/3)) may   
   > be optimal.  But if either a or b is 0, there exists better ones.   
   > Using the substitution u=x^2, the following antiderivative of x/((3 +   
   > x^2)*(1 - x^2)^(1/3)) is relatively trivial:   
   >   
   >     Sqrt[3]/(2*2^(2/3)) *   
   >       ArcTan[(1 + 2^(1/3)*(1 - x^2)^(1/3))/Sqrt[3]] -   
   >     1//(4*2^(2/3)) * Log[3 + x^2] +   
   >     3/(4*2^(2/3)) * Log[2^(2/3) - (1 - x^2)^(1/3)]   
   >   
   > The following antiderivative of 1/((3 + x^2)*(1 - x^2)^(1/3)) involves   
   > two arctangents and two logs, instead of Martin's two arctangents and   
   > four logs:   
   >   
   >     1/(2*2^(2/3)*Sqrt[3]) *   
   >       ArcTan[(1 - 2^(2/3)*(1 - x)^(2/3)/(1 + x)^(1/3))/Sqrt[3]] -   
   >     1/(2*2^(2/3)*Sqrt[3]) *   
   >       ArcTan[(1 - 2^(2/3)*(1 + x)^(2/3)/(1 - x)^(1/3))/Sqrt[3]] -   
   >     1/(4*2^(2/3)) * Log[(1 - x)^(2/3) + 2^(1/3)*(1 + x)^(1/3)] +   
   >     1/(4*2^(2/3)) * Log[2^(1/3)*(1 - x)^(1/3) + (1 + x)^(2/3)]   
   >   
      
   Your antiderivative for a = 0 implies a linear relation between its   
   single LN+LN+ATAN triplet and my pair of such triplets; this relation   
   can be used to eliminate either triplet from my antiderivative. In the   
   expressions thus obtained, one triplet is multiplied simply by a, the   
   other again by a linear combination of a and b. My antiderivative and   
   those two are comparable in size and therefore, I suppose, optimality.   
      
   For reasons of transparency, in my antiderivatives I usually ignore the   
   possibility of the triplet logarithms with polynomial arguments being   
   constant or compensating mutually, either fully or in part.   
      
   > I am trying to derive a general rule for integrating expressions of   
   > the form 1/((a + b*x^2)*(c + d*x^2)^(1/3)).  Can Martin's Goursat   
   > integrability checker determine the relationship, if any, between the   
   > parameters a, b, c and d that must hold in order for an elementary   
   > antiderivative to exist?  Of course, help anyone can provide with the   
   > rule itself would be much appreciated...   
   >   
      
   By rescaling the integration variable, removing a positive factor from   
   the radical, and removing another factor from the integral, the   
   parameters in your expression can be reduced to b = 1, c = +-1, d = +-1.   
   The integrand 1/((p + x^2)*(1 + x^2)^(1/3)) is Goursat pseudo-elliptic   
   if and only if (p + 3)*(p - 9) = 0, while 1/((p + x^2)*(1 - x^2)^(1/3))   
   and 1/((p + x^2)*(x^2 - 1)^(1/3)) are Goursat pseudo-elliptic if and   
   only if (p - 3)*(p + 9) = 0.   
      
   However, I do not consider it worthwhile to add such knowledge to Rubi.   
   In a rule-based integrator, the whole lot of Goursat pseudo-elliptics   
   can and should be caught in a single check of integrands predetermined   
   to be algebraic and to possess a suitable single radicand; integrands   
   found to be pseudo-elliptic would then be decomposed as needed and the   
   components rationalized by substitution, all others should be expressed   
   in terms of Legendre's canonical integrals.   
      
   Compared to making Rubi available on free algebra systems, all of this   
   is unimportant.   
      
   Martin.   
      
   --- SoupGate-Win32 v1.05   
    * Origin: you cannot sedate... all the things you hate (1:229/2)   

[   << oldest   |   < older   |   list   |   newer >   |   newest >>   ]


(c) 1994,  bbs@darkrealms.ca