       Re: a suprising result from Integrate (Null appeared in the

• To: mathgroup at smc.vnet.net
• Subject: [mg74526] Re: a suprising result from Integrate (Null appeared in the
• From: Peter Pein <petsie at dordos.net>
• Date: Sat, 24 Mar 2007 05:19:01 -0500 (EST)
• References: <eu1r2o\$jd6\$1@smc.vnet.net>

```Bhuvanesh schrieb:
> Yes, this was reported and fixed quite a while back, and should be fine in the next version. Another example of the bad behavior was Integrate[Sqrt[Sin[x] + Cos[x]], x]. For your indefinite integrals I currently get, in the development build:
>
until then, a simple linear shift of variables helps:

f[x_] = (1 - Sin[x])^(1/4);
approx = NIntegrate[f[x], {x, 0, Pi/2, 2*Pi}, WorkingPrecision -> 32]
inexact = Integrate[f[x], {x, 0, Pi/2, 2*Pi}]

5.699347567674386465367
0

exact =
FullSimplify[Integrate[FullSimplify[f[x + Pi/2]], {x, -Pi/2, 0, 3*(Pi/2)}]]

8*2^(1/4)*EllipticE[Pi/4, 2]

Chop[N[exact] - approx]
0

as an alternative to FullSimplify, use TrigToExp:

F[z_] = Simplify[Integrate[TrigToExp[f[x + Pi/2]], {x, -Pi/2, z - Pi/2}]]

z*((1/z)*(4*I - 4*(1 + I)^(3/2)*Hypergeometric2F1[-(1/4), 1/2, 3/4, -I]) +
(2*2^(3/4)*((I*(-I + E^(I*z))^2)/E^(I*z))^(1/4)*(-1 - I*E^(I*z) +
2*Sqrt[1 + I*E^(I*z)]*Hypergeometric2F1[-(1/4), 1/2, 3/4, (-I)*E^(I*z)]))/
((-I + E^(I*z))*z))

Plot[Chop[F[z]], {z, 0, 2*Pi}]
[omitted]

alternateexact = FullSimplify[
Subtract @@ (Limit[F[z], z -> Pi/2, Direction -> #1] & ) /@ {1, -1}]

(8*2^(1/4)*Sqrt[Pi]*Gamma[3/4])/Gamma[1/4]

Chop[N[alternateexact] - approx]
0

Peter

```

• Prev by Date: Re: Applying a list of functions to a single number
• Next by Date: Re: Definite Integration in Mathematica
• Previous by thread: Re: a suprising result from Integrate (Null appeared in the
• Next by thread: Re: a suprising result from Integrate (Null appeared in the