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 8,746 of 10,432   
   clicliclic@freenet.de to Waldek Hebisch   
   Re: Announce: FriCAS 1.2.4 has been rele   
   25 Jan 15 14:16:01   
   
   Waldek Hebisch schrieb:   
   >   
   > clicliclic@freenet.de wrote:   
   > >   
   > > ... and another quick FriCAS 1.2.3 on-line experiment:   
   > >   
   > > integrate((a + b*x + c*x^2)/((1 - x + x^2)*(1 - x^3)^(1/3)), x)   
   > >   
   > > >> Error detected within library code:   
   > >    integrate: implementation incomplete   
   > >    (residue poly has multiple non-linear factors)   
   > >   
   > > Is there no workable treatment of "multiple non-linear factors" other   
   > > than Kauer's Groebner-basis heuristics?   
   > >   
   >   
   > This one can be handled by simple method: since the integrand   
   > is linear in a, b, c each term can be integrated separately.   
   > FriCAS can integrate each part, however this is slow   
   > (I used 'setSimplifyDenomsFlag(true}' before integration to   
   > get resonable time) and result is ugly (sum over roots of   
   > sextic of large function).   
   >   
      
   Hmm, an integrator that can (incompletely) do individual terms but not   
   their sum. But this trick doesn't help with Example 21 from Timofeev   
   Chapter 4:   
      
   INT(((x - 1)^2*(x + 1))^(1/3)/x^2, x) = ?   
      
   Another way would be to split the rational factor into its complex   
   partial fractions:   
      
   (a + b*x + c*x^2)/(1 - x + x^2) = c   
    + (3*(b + c) + SQRT(3)*#i*(2*a + b - c))/(3*(2*x + SQRT(3)*#i - 1))   
    + (3*(b + c) - SQRT(3)*#i*(2*a + b - c))/(3*(2*x - SQRT(3)*#i - 1))   
      
   But this once again produces ">> Error detected within library code:   
   impossible".   
      
   However, the full result then doesn't involve any unresolvable roots:   
      
   INT(1/(1-x^3)^(1/3),x)=1/2*LN((1-x^3)^(1/3)+x)-1/SQRT(3)*ATAN(((~   
   1-x^3)^(1/3)-2*x)/(SQRT(3)*(1-x^3)^(1/3)))   
      
   INT(1/((2*x+SQRT(3)*#i-1)*(1-x^3)^(1/3)),x)=2^(2/3)*(1+SQRT(3)*#~   
   i)/32*(LN(4+((2*x+1-SQRT(3)*#i)/(2*(1-x^3)^(1/3)))^3)-3*LN(2^(2/~   
   3)+(2*x+1-SQRT(3)*#i)/(2*(1-x^3)^(1/3)))+2*SQRT(3)*ATAN(1/SQRT(3~   
   )*(1-(2*x+1-SQRT(3)*#i)/(2^(2/3)*(1-x^3)^(1/3)))))   
      
   For the third integral replace #i by -#i.   
      
   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