home bbs files messages ]

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,950 of 10,432   
   Nasser M. Abbasi to oldk1331@gmail.com   
   Re: if result of dsolve contains integra   
   08 Jan 16 13:50:55   
   
   From: nma@12000.org   
      
   On 1/6/2016 3:54 AM, oldk1331@gmail.com wrote:   
   >   
   > In my opinion, many integration results given by Mathematica are   
   > consisted of complicated non-elementary functions or general   
   > solutions with large expression of symbolic radicals, when   
   > there exist simpler elementary results. I don't think Mathematica   
   > solves those integrals.   
   >   
      
   Yes, I noticed this. Also Mathematica sometimes gives an   
   inverseFunction as result when there seems to be simpler   
   result, here is an example:   
      
   DSolve[-Sqrt[Abs[y[x]]] + y'[x]==0,y[x],x]   
      
   {{y[x] -> InverseFunction[   
               -((2*2^(3/4)*   
                Hypergeometric2F1[1/4, 3/4,   
                 7/4, (2*Abs[Im[#1]] + I*(-2 + Conjugate[#1] +   
                   #1))/(4*Abs[Im[#1]])]*(Abs[Im[#1]] +   
                   I*(1 - Re[#1]))^(1/4)*(1 + I*Abs[Im[#1]] - Re[#1])*(1 -   
   #1))/(3*   
                Abs[Im[#1]]^(1/4)*   
                   (Im[#1]^2 + (1 - Re[#1])^2))) - (2*2^(3/4)*   
              Hypergeometric2F1[1/4, 3/4, 7/4,   
                   (2*Abs[Im[#1]] +   
                  I*(2*Abs[#1]^2 - Conjugate[#1] - #1))/(4*Abs[Im[#1]])]*   
                   ((-I)*Abs[Im[#1]] + Im[#1]^2 - Re[#1] + Re[#1]^2   
   *((Abs[Im[#1]] +   
                   I*(-Im[#1]^2 + (1 - Re[#1])*Re[#1]))/(Im[#1]^2 + Re[#1]^2))^(1/   
                 4)*(1 - #1))/(3*   
              Abs[Im[#1]]^(1/4)*(Im[#1]^2 + (1 - Re[#1])^2)) & ][x + C[1]]}}   
      
   Maple:   
      
   dsolve(diff(y(x),x)-abs(y(x))^(1/2) = 0,y(x))   
     x-piecewise(y(x) <= 0, -2*sqrt(-y(x)), 0 < y(x), 2*sqrt(y(x)))+_C1 = 0   
      
   I think Mathematica anaswer here is a bit too complicated. May be   
   due to the how they implement Risch algorithm? I do not know.   
      
   > On the other hand, when FriCAS returns an integral unevaluated,   
   > it means that this integral has been proven to have no elementary   
   > result by Risch algorithm.  While when other systems return an   
   > integral unevaluated, it is possible that elementary answer   
   > exist but the system can't find it.   
   > Eg. Integrate[(5*x^9-5*x^6+9*x^5-5*x^4-x^2+1)/(Sqrt[x^5+x-1]*(   
   ^5-x^2+x-1)),x]   
   >   
      
   I wonder what does FriCAS give for the above ode?   
      
   --Nasser   
      
   --- 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