MathGroup Archive 2007

[Date Index] [Thread Index] [Author Index]

Search the Archive

Re: Definite Integration in Mathematica

  • To: mathgroup at smc.vnet.net
  • Subject: [mg74537] Re: Definite Integration in Mathematica
  • From: Peter Pein <petsie at dordos.net>
  • Date: Sat, 24 Mar 2007 05:25:05 -0500 (EST)
  • References: <etqo3f$10i$1@smc.vnet.net> <ett6i6$gd2$1@smc.vnet.net> <eu1qr3$ih3$1@smc.vnet.net>

dh schrieb:
> Hello Dimitris,
> 
> Yes, you are right, Mathematica chooses an antiderivative function that 
> 
> has branch cuts at approx. +/-2. However, this is not the only possible 
> 
> antiderivative, there are others with different branch cuts. But it is 
> 
> not too easy to find them.
> 
> Daniel
> 

But it is quite easy to do this in Mathematica:

In[1]:= << "DiscreteMath`"
In[2]:= f[x_] = (4 + 2*x + x^2) / (17 + 2*x - 7*x^2 + x^4);
In[3]:= F[x_] = Simplify[Assuming[k >= 0 && k \[Element] Integers,
  PowerSum[Simplify[ToRadicals[(x/(k + 1))*SeriesTerm[f[x], {x, 0, k}]]],
   {x, k}]]]
Out[3]=
(-(1/2))*I*(Log[34 + (1 + 4*I)*x - Sqrt[257 - 60*I]*x] +
   Log[34 + ((1 + 4*I) + Sqrt[257 - 60*I])*x] -
   Log[34 - ((-1 + 4*I) + Sqrt[257 + 60*I])*x] -
   Log[34 + ((1 - 4*I) + Sqrt[257 + 60*I])*x])

In[4]:= FullSimplify[F[4] - F[0]]
Out[4]= Pi - (1/2)*ArcTan[2752/825]
In[5]:= N[%]
Out[5]= 2.5018228707631676
In[6]:= Plot[F[x], {x, -4, 4}]
...

 This will propably work for "simple enough" (whatever that means) rational
functions.


Peter


  • Prev by Date: Re: a suprising result from Integrate (Null appeared in the
  • Next by Date: Re: Superimposing 3D plots with different colours
  • Previous by thread: Re: Definite Integration in Mathematica
  • Next by thread: Re: Definite Integration in Mathematica