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,583 of 10,432    |
|    clicliclic@freenet.de to Albert Rich    |
|    Re: Rubi 4.5 released    |
|    22 Jun 14 14:15:34    |
      Albert Rich schrieb:       >       > On Friday, June 20, 2014 4:56:06 PM UTC-10, Nasser M. Abbasi wrote:       >       > >       > > Int[ArcSin[Sqrt[x + 1] - Sqrt[x]], x]       > >       > > Out[68]= (-x)*ArcSin[Sqrt[x] - Sqrt[1 + x]] +       > > Subst[Int[Sqrt[1 - x^2 + x*Sqrt[-1 + x^2]], x], x, Sqrt[1 + x]]/Sqrt[2]       >       > As the above output shows, using integration by parts Rubi is able to       > find one term of the antiderivative (ie (-x)*ArcSin[Sqrt[x]). However       > it is unable to integrate the resulting algebraic integrand       >       > Sqrt[-x + Sqrt[x]*Sqrt[1 + x]]/Sqrt[1 + x]       >       > to get the antiderivative       >       > 1/2*(Sqrt[x] + 3*Sqrt[1 + x])*Sqrt[-x + Sqrt[x]*Sqrt[1 + x]] -       > 3*ArcSin[Sqrt[x] - Sqrt[1 + x]]/(2*Sqrt[2])       >       > Mathematica is able to integrate it but gets an antiderivative 3 times       > the above size. I do not know how to integrate this algebraic       > integrand. Tell me the integration steps you would use to integrate       > it, and I will be happy to teach Rubi how to do it.       >              Perhaps the best way is to manipulate the complex roots, as we did when       Rubi 2 was taught to handle Bondarenko's 'An exact 1-D integration       challenge - 58 - (sqrt)' in 2010. For arbitrary complex x we have:              SQRT(x + 1) = SQRT(1 - #i*SQRT(x))*SQRT(1 + #i*SQRT(x))              SQRT(SQRT(x)*(SQRT(x + 1) - SQRT(x))) =        SQRT(SQRT(x))*SQRT(SQRT(x + 1) - SQRT(x))              SQRT(2)*SQRT(SQRT(x + 1) - SQRT(x)) =        SQRT(#i)*SQRT(1 + #i*SQRT(x)) + SQRT(-#i)*SQRT(1 - #i*SQRT(x))              The first identity is well known, the second holds since       |PHASE(SQRT(x))| <= pi/2 and |PHASE(SQRT(x+1) - SQRT(x))| < pi/2, and       the third was derived in 2010. So we can write:              INT(SQRT(SQRT(x)*(SQRT(x + 1) - SQRT(x)))/SQRT(x + 1), x)        = INT(SQRT(SQRT(x))*SQRT(SQRT(x + 1) - SQRT(x))/SQRT(x + 1), x)        = INT(SQRT(SQRT(x))*(SQRT(#i)*SQRT(1 + #i*SQRT(x))        + SQRT(-#i)*SQRT(1 - #i*SQRT(x)))        / (SQRT(2)*SQRT(1 - #i*SQRT(x))*SQRT(1 + #i*SQRT(x))), x)        = INT(SQRT(#i)*x^(1/4)/(SQRT(2)*SQRT(1 - #i*SQRT(x))), x)        + INT(SQRT(-#i)*x^(1/4)/(SQRT(2)*SQRT(1 + #i*SQRT(x))), x)              where the final integrals should be no problem for Rubi. The alternative       would be to use variable substitution; one way is described by       Charlwood, another is to apply the substitution t = SQRT(x + 1) -       SQRT(x), SQRT(x) = (1 - t^2)/(2*t), dx = (t^4 - 1)/ (2*t^3)*dt to the       integral:               INT(SQRT(SQRT(x))*SQRT(SQRT(x + 1) - SQRT(x))/SQRT(x + 1), x)              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