Re: (N)Integrate help - convergence problem
- To: mathgroup at smc.vnet.net
- Subject: [mg48795] Re: (N)Integrate help - convergence problem
- From: Paul Abbott <paul at physics.uwa.edu.au>
- Date: Wed, 16 Jun 2004 07:48:56 -0400 (EDT)
- Organization: The University of Western Australia
- References: <ca928h$cc$1@smc.vnet.net>
- Sender: owner-wri-mathgroup at wolfram.com
In article <ca928h$cc$1 at smc.vnet.net>, Patrick Jemmer <patrick.jemmer at unn.ac.uk> wrote: > (* WAVEFUNCTION *) > > Fit5g = > > Sqrt[ > 0.0711280497855832 / E^(16*r^2) + > 0.0950242005793304 / E^(4*r^2) + > 0.09001525518121072 / E^r^2 + > 0.011332742959622262 / E^(r^2/4) - > 0.0003397835348217242 / E^(r^2/16) > ]; If you try to plot this wavefunction over, say, 0 < r < 5, I think you'll see what the problem is. The radial wavefunction is not always real. This is because the last term in the square root, -0.0003397835348217242 / E^(r^2/16) dominates for sufficiently large r, leading to complex values. 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