Services & Resources / Wolfram Forums
-----
 /
MathGroup Archive
2004
*January
*February
*March
*April
*May
*June
*July
*August
*September
*October
*November
*December
*Archive Index
*Ask about this page
*Print this page
*Give us feedback
*Sign up for the Wolfram Insider

MathGroup Archive 2004

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

Search the Archive

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


  • Prev by Date: Using NDSolve in NonlinearFit
  • Next by Date: Re: Integrating UnitSteps
  • Previous by thread: (N)Integrate help - convergence problem
  • Next by thread: 3D plots from lists