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,220 of 10,432   
   clicliclic@freenet.de to Albert Rich   
   Re: integration of sec(t)   
   02 Dec 16 18:17:23   
   
   Albert Rich schrieb:   
   >   
   > On Monday, November 21, 2016 at 6:35:36 PM UTC-10, Nasser M. Abbasi wrote:   
   > > It is interesting how different CAS outputs differ sometimes.   
   > > I was trying to verify book result for int(sec(t),t) and this   
   > > is what some CAS systems gave   
   > >   
   > > Mathematica:     ln(cos(t/2)+sin(t/2))-ln(cos(t/2)-sin(t/2))   
   > > Rubi:            arctanh(sin(t))   
   > > Fricas:          1/2 ( ln(1+sin(t)) - ln(1-sin(t)) )   
   > > Maple & Maxima:  ln(sec(t)+tan(t))   
   > >   
   > > [...]   
   > >   
   > > Was Wondering if there is any preference, math-wise, to any one   
   > > of the above results. From the Mathematica result, it seems   
   > > to hint that it used the Tangent half-angle substitution,   
   > > also called Weierstrass substitution.   
   > >   
   > > [...]   
   >   
   > Note that in addition to being the most compact, Rubi's antiderivative   
   > is always real when t is real, unlike the ones returned by   
   > Mathematica, Maple and Maxima. Also the other antiderivatives involve   
   > multiple instances of the integration variable.  This can result in   
   > catastrophic cancellation when these antiderivatives are numerically   
   > evaluated.   
   >   
   > However, if you still prefer an alternative antiderivative for sec(t),   
   > note that Rubi is an open-source, modular system.  So, all you have to   
   > do is change just one of Rubi's 6800+ integration rules to get the   
   > desired effect.  Specifically the rule   
   >   
   >     int(sec(c+d x),x) --> arctanh(sin(c+d x))/d   
   >   
      
   A zero imaginary part for Rubi's antiderivative on the entire real line   
   comes at the price of having two branch cuts terminate at each of its   
   pole locations t = n*pi/2 (n = +-1, +-3, +-5, ...) instead of just one   
   branch cut. One may consider this a reason to change the current Rubi   
   rule.   
      
   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