Re: Bug 3! Gaussian integration in 5.0 is broken!
- To: mathgroup at smc.vnet.net
- Subject: [mg43777] Re: Bug 3! Gaussian integration in 5.0 is broken!
- From: Paul Abbott <paul at physics.uwa.edu.au>
- Date: Sat, 4 Oct 2003 02:04:55 -0400 (EDT)
- Organization: The University of Western Australia
- References: <blj66o$1i8$1@smc.vnet.net>
- Sender: owner-wri-mathgroup at wolfram.com
In article <blj66o$1i8$1 at smc.vnet.net>,
Konstantin L Kouptsov <kouptsov at wsu.edu> wrote:
> Check this:
>
> Integrate[Exp[(I/h)*(x*Sqrt[B/2 + D] + x0)^2], {x, -Infinity, Infinity}]
>
> and observe the Erf[] functions that should not be there. Neither should x0.
However, note the assumptions that Integrate makes. On the other hand,
Assuming[x0 > 0 && B > 0 && D > 0 && h > 0,
Integrate[E^((I*(Sqrt[B/2 + D]*x + x0)^2)/h),
{x, -Infinity, Infinity}]]
yields a reasonable result.
Cheers,
Paul
--
Paul Abbott Phone: +61 8 9380 2734
School of Physics, M013 Fax: +61 8 9380 1014
The University of Western Australia (CRICOS Provider No 00126G)
35 Stirling Highway
Crawley WA 6009 mailto:paul at physics.uwa.edu.au
AUSTRALIA http://physics.uwa.edu.au/~paul