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 9,253 of 10,432   
   clicliclic@freenet.de to oldk1331@gmail.com   
   Re: is your integrator output Davenport    
   31 Jan 17 17:39:51   
   
   oldk1331@gmail.com schrieb:   
   >   
   > > >   
   > > > 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.   
   >   
   > Hi Martin, FriCAS does an extra simplification, causes the result to be   
   > discontinuous at x=0.  This is caused by the same bug you mentioned   
   > in the previous thread.  I think there is a good enough fix, and the   
   > result after patching this bug is continuous:   
   >   
   > (log((2*x*((-1)*x^4+1)^(1/2)+((-1)*x^4+2*x^2+1))/(x^4+1))+2*atan((x^2*((-1)*   
   >   x^4+1)^(1/2)+(x^3+x))/(((-1)*x^4+1)^(1/2)+((-1)*x^3+x))))/4   
   >   
      
   This can still be simplified considerably:   
      
     1/2*ATANH(x*(1 - x^2)/SQRT(1 - x^4))   
     + 1/2*ATAN(x*(1 + x^2)/SQRT(1 - x^4))   
      
   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