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,755 of 10,432    |
|    clicliclic@freenet.de to Waldek Hebisch    |
|    Re: Announce: FriCAS 1.2.4 has been rele    |
|    15 Feb 15 19:11:03    |
      Waldek Hebisch schrieb:       >       > clicliclic@freenet.de wrote:       > >       > > clicliclic@freenet.de schrieb:       > > >       > > > integrate((2 - (k + 1)*x)/((1 - (k + 1)*x)       > > > *(x*(1 - x)*(1 - k*x))^(1/3)), x)       > > >       > > > >> Error detected within library code:       > > > impossible       > > >       > > > It is good to be cognizant of one's limitations ;).       > > >       > >       > > ... 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?       >       > Beware that Kauer's heuristics is incomplete. AFAICS the       > your first example is one when this heurisctic fail. Namely,       > rationalizing substitution gives integral of the form       >       > \sum w*log(x - w*(x*(1 - x)*(1 - k*x))^(1/3))       >       > where w runs over cube roots of 1/k. Derivative of the       > log has pole at 0 and at point over 1/(k+1). Trager       > method depends on rearanging integral in such a form that       > logands are regular at 0. So Trager method produces       > 3 divisors, each one corresponding to point with       > x = 1/(k+1) and y beeing one ot the roots of (x*(1 - x)*(1 - k*x)).       > The same happens in Kaures method. Next Kaures tries       > to find function having zero as specified by the divisor       > or its power and arbitrary behaviour at infinity. For this       > to succeed the divisor have to be torsion divisor. But       > this is unlikely to be the case. At least my partial       > implementation of Kauers method gives up on this example.       >       > The example is amenable to Trager method. If I replace k       > by k^3, then FriCAS can compute the integral. Such replacement       > is equivalent to extending field of constants by cube root       > of k and FriCAS could do this automatically. It is not       > done ATM because this can lead to excessive execution time       > (in this case it is 287.49 sec on fast machine). I did       > not analyze the second example, but it seem to be similar.       >       > Let me add that that trying Kauers method I substituted       > 5 for k. Without such substitution Groebner basis       > computations were quite slow.       >              Wow, Kauers is done in! [How about publishing a corresponding note       explaining why torsion divisors are needed here and Kauers' Groebner       bases cannot furnish them? (I could supply the theory of Goursat       pseudo-elliptics for cube roots of quadratics and cubics.)]              Yes, my antiderivative of the first example contains (-k)^(1/3):              INT((2 - (k + 1)*x)/((1 - (k + 1)*x)*(x*(1 - x)*(1 - k*x))^(1/3)), x)        = - 1/(2*(-k)^(1/3))*LN(x*(1 - x)*(1 - k*x) - k*x^3)        + 3/(2*(-k)^(1/3))*LN((x*(1 - x)*(1 - k*x))^(1/3) + (-k)^(1/3)*x)        - SQRT(3)/(-k)^(1/3)        *ATAN(1/SQRT(3)*(1 - 2*(-k)^(1/3)*x/(x*(1 - x)*(1 - k*x))^(1/3)))              Expecting the FriCAS subsubversion for an improved Trager procedure to       be x < 9,              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