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,825 of 10,432   
   clicliclic@freenet.de to Nasser M. Abbasi   
   Re: what special about 18 in this integr   
   23 Feb 18 08:45:14   
   
   "Nasser M. Abbasi" schrieb:   
   >   
   > On 2/22/2018 1:24 PM, clicliclic@freenet.de wrote:   
   > >   
   > > "Nasser M. Abbasi" schrieb:   
   > >>   
   > >> On 2/22/2018 7:57 AM, clicliclic@freenet.de wrote:   
   > >>>   
   > >>> "Nasser M. Abbasi" schrieb:   
   > >>>>   
   > >>>> How high a number "n" can your cas integrate this?   
   > >>>>   
   > >>>>           1/((1 + x^n)*(1 + x^2))   
   > >>>>   
   > >>>> Fricas 1.3.2 and Rubi 4.14.7 go up to n=18. Then after that, they   
   > >>>> return unevaluated.   
   > >>>>   
   > >>>> Mathematica gives an answer for higher than 18, but the result   
   > >>>> after 18 is in terms of Root objects. I tried up to n=100.   
   > >>>>   
   > >>>> Maple also can go higher than n=18, and it also gives results   
   > >>>> in terms of RootOf. Tried up to n=200. So to remove these Roots,   
   > >>>> one has to evaluate the answer numerically.  I assume this is why   
   > >>>> Rubi and Fricas stop at n=18.   
   > >>>>   
   > >>>> Any insight why 18 is the limit here? Is it due to some   
   > >>>> factorization done, which after n=18 produces polynomials that   
   > >>>> can't be solved exactly for higher order?   
   > >>>>   
   > >>   
   > >> [...]   
   > >>   
   > >> Sorry, let me clarify things again. FriCAS hangs at n=29.   
   > >> I have it running for one hr now at n=29. Will leave it   
   > >> running and will check on it after I come back from school.   
   > >>   
   > >> Here is what I found:   
   > >> ======================   
   > >> Rubi:  Up to n=18 OK, then unevaluated after n=18.   
   > >> Maple: Up to n=18 OK, then uses ROOT objects for n>18. Tried to 200   
   > >> Mathematica: Up to n=18 OK, then uses ROOT objects for n>18. Tried   
   > >> to 100 FriCAS: Up to n=28 OK, then "hangs" or still trying.....   
   > >>   
   > >> So the winner so far is FriCAS on this test, it does it   
   > >> for n=28, with no root objects. root objects can only be   
   > >> evaluated numerically to get rid of them.   
   > >>   
   > >   
   >   
   > > I am glad to learn that the super-massive FriCAS bug has vanished   
   > > at a second glance! In fact, I see no problem here at all. While it   
   > > remains uncertain how long the computation for n=29 takes, this is   
   > > likely to happen for some power anyway ...   
   > >   
   >   
   > Yes, sorry, I must have mixed it up with 1/((1 + x^n)*(1 + x^2))   
   > for `n` being symbolic (no numerical value). In this call,   
   > all CAS systems return unevaluated.   
   >   
   > I like to also report that Fricas finally gave answer for n=29.   
   > many pages long. But it did not hang.   
   >   
   > good job Fricas.   
   >   
      
   It may be appropriate here to recall another integrand for which   
   something basic went wrong inside FriCAS:   
      
   integrate(1/((p*sqrt(1 - x^2) - 1)*sqrt(p*sqrt(1 - x^2) - x^2)), x)   
      
   which used to exit with:   
      
   >> Error detected within library code:   
      Denominator not equal to 1   
      
   If this integral is found to be non-elementary, as I expect it to, it   
   should simply be returned unevaluated.   
      
   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