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,246 of 10,432    |
|    clicliclic@freenet.de to Nasser M. Abbasi    |
|    Re: is your integrator output Davenport     |
|    24 Jan 17 12:04:05    |
      "Nasser M. Abbasi" schrieb:       >       > On 1/24/2017 1:33 AM, clicliclic@freenet.de wrote:       > >       > > "Never forget to check that the output is a continuous function!"       > >       > > Thus James H. Davenport reminded us in a conference presentation on the       > > "Complexity of Integration, Special Values, and Recent Developments" at       > > last year's International Congress on Mathematical Software.       > >       > > Now, in a posthumous publication in the Institutiones Calculi Integralis       > > Tomus 4, pages 36 - 48, Leonhard Euler showed that:       > >       > > SQRT(a + b*x^2 + c*x^4)/(a - c*x^4)       > >       > > has an elementary antiderivative. In fact, the integrand is Goursat       > > pseudo-elliptic.       > >       > > Does your favorite symbolic integrator then succeed in returning a (i)       > > real, (ii) elementary, and (iii) continuous result for:       > >       > > INT(SQRT(1 - x^4)/(1 + x^4), x)       > >       > > ? Yes, a simple such a solution does exist.       > >       > > Martin.       > >       > > PS: Derive 6.10 just returns the integral with its denominator factored       > > into real quadratics.       > >       >       > Maple 2016.2       > ===============       >       > int( sqrt(1 - x^4)/(1 + x^4), x);       >       > -(1/4)*arctan((-x^4+1)^(1/2)/x+1)+(1/4)*       > arctan(-(-x^4+1)^(1/2)/x+1)-(1/8)*       > ln(((1/2)*(-x^4+1)/x^2-(-x^4+1)^(1/2)/x+1)/((1/2)*(-x^4+1)/x^2       (-x^4+1)^(1/2)/x+1))              This is (i) real and (ii) elementary, but (iii) discontinuous at x=0.       There are also imaginary discontinuities at x^2 = 1 + SQRT(2).              >       > Mathematica 11.0.1       > ==================       > Integrate[ Sqrt[1 - x^4]/(1 + x^4), x]       > x AppellF1[1/4, -(1/2), 1, 5/4, x^4, -x^4]              This is not elementary. I am not surprised.              >       > Rubi 4.10       > ===========       > ShowSteps = False;       > Int[ Sqrt[1 - x^4]/(1 + x^4), x]       >       > -EllipticF[ArcSin[x], -1] + EllipticPi[-I, ArcSin[x], -1] +       > EllipticPi[I, ArcSin[x], -1]              This is neither real nor elementary. I am surprised.              >       > FriCAS 1.3       > ==========       > integrate( sqrt(1 - x^4)/(1 + x^4), x)       >       > (log((2*x*((-1)*x^4+1)^(1/2)+((-1)*x^4+2*x^2+1))/(x^4+1))+((-2)*atan((((-1)*       > x^4+1)^(1/2)+x)/x)+2*atan(x^2)))/4              This is (i) real and (ii) elementary, but (iii) discontinuous at x=0.       There are no other discontinuities.              -              Only Maple 2016.2 and FriCAS 1.3 return elementary answers, but both are       discontinuous at x=0.              More work is therefore needed to make the output of these integrators       Davenport compatible.              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