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,886 of 10,432    |
|    clicliclic@freenet.de to clicliclic@freenet.de    |
|    Re: elementarily integrable or not?    |
|    20 Apr 18 21:18:27    |
      clicliclic@freenet.de schrieb:       >       > antispam@math.uni.wroc.pl schrieb:       > >       > > clicliclic@freenet.de wrote:       > > >       > > > clicliclic@freenet.de schrieb:       > > > >       > > > > Here's a much more deadly failure:       > > > >       > > > > setSimplifyDenomsFlag(true)       > > > >       > > > > integrate((7*x + 6*sqrt(17*sqrt(2) + 23) + 9*sqrt(2) + 34)       > > > > /((x^2 + 2*x*(3*sqrt(sqrt(2) + 1) + 1) + 6*sqrt(58*sqrt(2) + 82)       > > > > - 18*sqrt(2) - 26)*sqrt(x^3 - 30*x - 56)), x)       > > > >       > > > > is returned unevaluated on the web interface. But my oracle says that       > > > > the integand is pseudo-elliptic according to Goursat:       > > > >       > > > > goursat3((7*x + 6*SQRT(17*SQRT(2) + 23) + 9*SQRT(2) + 34)       > > > > /(x^2 + 2*x*(3*SQRT(SQRT(2) + 1) + 1) + 6*SQRT(58*SQRT(2) + 82)       > > > > - 18*SQRT(2) - 26), x, -56, -30, 0, 1)       > > > >       > > > > [false, false, false, true]       > > > >       > > > > The last answer counts.       > > > >       > > >       > > > Since the oracle was characteristically cryptic, I should supply some       > > > details. Putting x = (6 - 3*SQRT(2))*t - 4 one finds:       > > >       > > > SQRT(x^3 - 30*x - 56) =       > > > 3*SQRT(12 - 6*SQRT(2))*SQRT(t*(1 - t)*(1 - (SQRT(2) - 1)^2*t))       > > >       > > > So Legendre's modulus k of the elliptic radical takes the nice special       > > > value SQRT(2) - 1. The quadratic in the denominator of the integrand       > > > divides the division polynomial of order four for the elliptic curve       > > > y^2 = x^3 - 30*x - 56. The antiderivative consists of two terms:       > > >       > > > 6*(2*SQRT(SQRT(2) + 1) - SQRT(2))/7       > > > * INT((7*x + 6*SQRT(17*SQRT(2) + 23) + 9*SQRT(2) + 34)       > > > / ((x^2 + 2*x*(3*SQRT(SQRT(2) + 1) + 1) + 6*SQRT(58*SQRT(2) + 82)       > > > - 18*SQRT(2) - 26)*SQRT(x^3 - 30*x - 56)), x)       > > > = 2/(SQRT(3*SQRT(2) + 6) - SQRT(3)*2^(3/4))       > > > * ATANH((SQRT(3*SQRT(2) + 6) - SQRT(3)*2^(3/4))*(x - 3*SQRT(2) - 2)       > > > / SQRT(x^3 - 30*x - 56)) - SQRT(6 - 3*SQRT(2))/3       > > > * ATAN(SQRT(6*SQRT(2) + 12)*(x + 4)/SQRT(x^3 - 30*x - 56))       > > >       > > > Having this integral declared non-elementary constitutes a bug that       > > > seriously undermines trust in the FriCAS integrator.       > > >       > >       > > You hit known limitation of current implementation. Your integral       > > contains dependent roots. If I eliminate dependent roots as below:       > >       > > setSimplifyDenomsFlag(true)       > > ff := (7*x + 6*sqrt(17*sqrt(2) + 23) + 9*sqrt(2) + 34) _       > > /((x^2 + 2*x*(3*sqrt(sqrt(2) + 1) + 1) + 6*sqrt(58*sqrt(2) + 82) _       > > - 18*sqrt(2) - 26)*sqrt(x^3 - 30*x - 56))       > >       > > s2 := sqrt(2)*sqrt(sqrt(2) + 1)^5       > > s3 := (sqrt(2) + 3)*sqrt(sqrt(2) + 1)       > > ff2 := eval(ff, [sqrt(58*sqrt(2) + 82) = s2, sqrt(17*sqrt(2) + 23) = s3])       > > integrate(ff2, x)       > >       > > I get:       > >       > > [solution snipped]       > >       > > which is large but seem to correspond to your desired value. OTOH       > > there are 3 other combinations of signs of roots and FriCAS returns       > > unevaluated for those combinations.       > >       > > Note: FriCAS should find out that roots are dependent, but this part       > > currently remains unimplemented.       >       > Yes, the only numerical roots appearing here are sqrt(2) and       > sqrt(sqrt(2) + 1), corresponding to sqrt(w) and sqrt(w - 1) in the PS       > to my post of 04 Apr 2018, but I didn't connect this hazy background       > fact with my rather clear knowledge that FriCAS may give wrong results       > in the case of dependent roots. My interest was simply in checking the       > integrability of certain algebraics; this question has been resolved       > without using FriCAS once its results were found unreliable, by the       > way.       >       > Does your statement about the signs of roots mean that the correct       > method of reducing dependent roots must be found by trying different       > sign combinations? Your s2, s3 look very natural to me.       >              The sign problem may reflect that there is just a single numerical root       of degree four in this integrand. Its rational part:              (7*x + 6*sqrt(17*sqrt(2) + 23) + 9*sqrt(2) + 34)       /(x^2 + 2*x*(3*sqrt(sqrt(2) + 1) + 1) + 6*sqrt(58*sqrt(2) + 82)        - 18*sqrt(2) - 26)              multiplied by 6*(- sqrt(sqrt(2) + 1) + sqrt(2) + 2)/7 equals:              (3*(t^2 + 8*t + 22)*x - t^3 + 12*t^2 + 90*t + 392)       /((x - t)*((2*t + 8)*x + t^2 + 4*t + 18))              for t = - 6*sqrt(sqrt(2) + 1) + 3*sqrt(2) + 2, one among four roots of       the quartic t^4 - 8*t^3 - 84*t^2 - 464*t - 956 = 0.              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