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,259 of 10,432   
   Albert Rich to clicl...@freenet.de   
   Re: is your integrator output Davenport    
   01 Feb 17 21:52:30   
   
   From: Albert_Rich@msn.com   
      
   On Tuesday, January 31, 2017 at 6:40:08 AM UTC-10, clicl...@freenet.de wrote:   
   > 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*   
   tan((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.   
      
   Nice!  Is it possible to Riobooize the arctan term so the result is continuous   
   at x^2=1?   
      
   Albert   
      
   --- 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