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,235 of 10,432    |
|    clicliclic@freenet.de to antispam@math.uni.wroc.pl    |
|    Re: alarum: Risch integrator fails to di    |
|    14 Jan 17 23:40:51    |
      antispam@math.uni.wroc.pl schrieb:       >       > clicliclic@freenet.de wrote:       > >       > > A state-of-the-art implementation of a Risch integrator       > >       > > )version       > >       > > Value = "FriCAS 1.3.0 compiled at Thu Sep 1 17:01:37 UTC 2016"       > >       > > responds to       > >       > > integrate((1 + x)/((1 + x + x^2)*(a + b*x^3)^(1/3)), x)       > >       > > with       > >       > > >> Error detected within library code:       > > catdef: division by zero       > >       > > In his "Symbolic Integration Tutorial", Bronstein appears to believe in       > > the existence of a workable integration algorithm for algebraic       > > integrands involving a simple radical extension. Didn't he know what he       > > was talking about?       >       > Writing paper he ceratainly knew what is is writing. But writing       > code he apparently assumed that he will be lucky. Or to put it       > differetly: code he wrote is enough to "demostrate" that algorithm       > works. But there remain "uninteresting" programming probles to       > resolve. To the point: his code uses random integers. If       > substituting those integers into one of denominators gives zero,       > then you will see this error. Correct algorithm would retry       > with different values. But his code just gives error. Recent       > changes to FriCAS lead to change in random choice and for       > some integrals (like this one) the problem became very visible.       >       > One can try to hide the problem using more randomness. Or       > one can do "uninteresting" programming and properly implement       > retry.       >       > BTW: This is FriCAS bug and proper place to report this is       > FriCAS mailing list or Sourceforge bugtracke. And this       > problem was was already reported (after relase of 1.3.0).       >              Oh, I thought I was the first to try this particular variation of a       familar cube-root integrand. But to implement a loop over different       pseudo-random numbers should be easy, and its omission may be a mere       oversight. I will submit this integrand again to the next version.              Another thing has me wondering. Your antiderivatives of rational       functions can be discontinuous because the ATAN part is not always       rewritten to involve polynomial arguments only, as proposed by Rioboo.       For example, FriCAS 1.3.0 computes (expressed in Derive notation)              INT((2560*x^3 - 400*x^2 - 576*x - 84)/(320*x^4 + 80*x^3 - 12*x^2        + 24*x + 9), x) = 2*LN(320*x^4 + 80*x^3 - 12*x^2 + 24*x + 9)        + 2*SQRT(11)*ATAN(SQRT(11)*(10*x + 3)/(80*x^2 + 10*x - 9))              This is surprising especially since Bronstein provides pseudocode of       Rioboo's algorithm in his book (Ch. 2.8).              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