MathGroup Archive 2008

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

Search the Archive

Re: Incorrect integral in Mathematica

  • To: mathgroup at smc.vnet.net
  • Subject: [mg91436] Re: Incorrect integral in Mathematica
  • From: Bhuvanesh <BhuvaneshBhatt at gmail.com>
  • Date: Thu, 21 Aug 2008 05:56:44 -0400 (EDT)

This was a bug. Essentially, the following should have Sqrt[Pi]/Sqrt[a], but it was incorrectly giving zero in version 5.0.x:

Integrate[Exp[-a*(y - y0)^2], {y, -Infinity, Infinity}, GenerateConditions->False]

Here are a couple of workarounds:

$Version      = 5.0.1 for Linux (November 18, 2003)

In[1]:= Integrate[Exp[-x^2/(2w^2)]Exp[-(y - y0)^2/(2z^2)],
{x, -Infinity, Infinity}, {y, -Infinity, Infinity}, GenerateConditions->True] /. If[_,res_,_]:>res        //InputForm

Integrate::gener: Unable to check convergence.

Out[1]//InputForm=
Sqrt[2*Pi]*Sqrt[z^2]*If[Re[w^2] > 0, Sqrt[2*Pi]*Sqrt[w^2],
  Integrate[E^(-x^2/(2*w^2)), {x, -Infinity, Infinity},
   Assumptions -> Re[w^2] <= 0]]

In[2]:= Integrate[Exp[-x^2/(2w^2)]Exp[-(y - y0)^2/(2z^2)], {x, -Infinity, Infinity}, {y, -Infinity, Infinity}, Assumptions->z>0] //InputForm

Out[2]//InputForm=
Sqrt[2*Pi]*z*If[Re[w^2] > 0, Sqrt[2*Pi]*Sqrt[w^2],
  Integrate[E^(-x^2/(2*w^2)), {x, -Infinity, Infinity},
   Assumptions -> z > 0 && Re[w^2] <= 0]]

The integral works correctly in version 5.1 and newer versions. Sorry for the inconvenience.

Bhuvanesh,
Wolfram Research


  • Prev by Date: Re: Multifile search
  • Next by Date: Re: Installing package documentation
  • Previous by thread: Re: Incorrect integral in Mathematica
  • Next by thread: Fastest way to check inequality