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: [mg9254] Re: Another Bug in Mathematica 3.0.0 definite integration
  • From: carlos at mars.Colorado.EDU (Carlos A. Felippa)
  • Date: Fri, 24 Oct 1997 01:01:07 -0400
  • Organization: University of Colorado, Boulder
  • Sender: owner-wri-mathgroup at wolfram.com

In article <62hfha$m06$4 at dragonfly.wolfram.com> "Gregor Overney" <overney at worldnet.att.net> writes:
>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
>
>
>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}]
>>
>>(Note that the integrand is definite positive in the integration range)
>>
>>Mathematica3.0.0 returns
>>
>>     a= -2 2^(3/4) HypergeometricPFQ[{1/4,3/4},{5/4},-1]
>>
>>and
>>     N[a]=-3.01236...
>>
>>With the standard change of variable t=Tan[x/2] the integral can be
>>easily evaluated symbolically,  and then the numeric evaluation returns
>>
>>              1.3974.....
>>
>>
>>in perfect agreement with the result obtained performing directly
>>NIntegrate on the original integrand.
>>

This is a pattern I have observed in many integrals involving
trigonometric functions that lead to hypergeometric and/or elliptic
functions.  The time comparison of 2.2 vs 3.0 is interesting.

Results from Mathematica 2.2.1 on above test function:

(Pi^(1/2)*Gamma[1/4])/(2*2^(1/4)*Gamma[3/4]) - 
   2*Hypergeometric2F1[1/2, 1, 5/4, -1] + Hypergeometric2F1[3/4, 1, 3/2,
-1]

   1.397395299268851

Time: 0.95 seconds on Mac 8500/120.

Results from Mathematica 3.0.1, same machine:

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

  -3.012362296717479

Time: 10.4 seconds, a factor of 10.  

Both 3.0 and 3.0.1 give the warning about convergence.


  • Prev by Date: Re: Another Bug in Mathematica 3.0.0 definite integration
  • Next by Date: Finding "start directory" with Mathematica 3.0 ?
  • Previous by thread: Re: Another Bug in Mathematica 3.0.0 definite integration
  • Next by thread: Re: Re: Another Bug in Mathematica 3.0.0 definite integration