MathGroup Archive 1997

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

Search the Archive

Re: Another Bug in Mathematica 3.0.0 definite integration

  • To: mathgroup at smc.vnet.net
  • Subject: [mg9215] Re: Another Bug in Mathematica 3.0.0 definite integration
  • From: weber at math.uni-bonn.de (Matthias Weber)
  • Date: Fri, 24 Oct 1997 01:00:30 -0400
  • Organization: RHRZ - University of Bonn (Germany)
  • Sender: owner-wri-mathgroup at wolfram.com

In article <62hfha$m06$4 at dragonfly.wolfram.com>, "Gregor Overney"
<overney at worldnet.att.net> wrote, commenting on :

> luca ciotti wrote in message <624fv1$les at smc.vnet.net>...
> >Dear Users,
> >
> >unfortunately I found another erroneous result in  a definite integral
> >in Mathematica 3.0.0
> >
> >Let
> >
> >        a=Integrate[1/Sqrt[Sin[x]+Cos[x]], {x,0,Pi/2}]
> >


> Mathematica 3.0.1.1x would give you at least a warning, suggesting to
> carefully check the convergence.
> 
> your input produces:
> 
> Integrate::gener: Unable to check convergence
> 
> and N[a] gives the obviously wrong value of -3.0123622967174799.
> 
> GTO
> 

Mathematica 3.0 gives the same messages already. This were also no
improvement, because for this integral there is no convergence to
check. And, in a previous post, somebody had an obviously finite
integral which Mathemetica declined to evaluate because of that
´convergence test´. A response by Wolfram recommended to turn this test
off (in general? in that case?). At least here, it is totally nonsense.


Concerning the actual problem: The primitive of the integrand
  1/Sqrt[Sin[x]+Cos[x]
involves hypergeometric functions (after Simplify):

(Cot[(Pi - 4*x)/4]*Hypergeometric2F1[1/4, 1/2, 5/4, 
     (1 + Sin[2*x])/2]*Sqrt[2 - 2*Sin[2*x]])/
  Sqrt[Cos[x] + Sin[x]]


and these being multivalued in general, it is conceivable that
Mathematica is ignoring the branching behaviour while computing the
definite integral. This is supported by the Plot of the above
expressiion which clearly has a jump at x=Pi/4. If you take this into
account by hand (no comments), you get the explicit value

(2*2^(3/4)*Sqrt[Pi]*Gamma[5/4])/Gamma[3/4] - 
  2*2^(3/4)*HypergeometricPFQ[{1/4, 3/4}, {5/4}, -1]

which has the correct numerical value.

Is there anybody who has received information on this bug by Wolfram, 
the impact of it on other integrals, possible workarounds and whether
they are willing to provide a fix?

Matthias Weber


  • Prev by Date: Re: Help - Mathematica bombs out on loading <<Calculus`LaplaceTransform`
  • Next by Date: Re: Another Bug in Mathematica 3.0.0 definite integration
  • Previous by thread: Re: Another Bug in Mathematica 3.0.0 definite integration
  • Next by thread: Re: Another Bug in Mathematica 3.0.0 definite integration