       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:= << "DiscreteMath`"
In:= f[x_] = (4 + 2*x + x^2) / (17 + 2*x - 7*x^2 + x^4);
In:= F[x_] = Simplify[Assuming[k >= 0 && k \[Element] Integers,
PowerSum[Simplify[ToRadicals[(x/(k + 1))*SeriesTerm[f[x], {x, 0, k}]]],
{x, k}]]]
Out=
(-(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:= FullSimplify[F - F]
Out= Pi - (1/2)*ArcTan[2752/825]
In:= N[%]
Out= 2.5018228707631676
In:= 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