MathGroup Archive 2012

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

Search the Archive

Finding Converging Function

I have found a solution to a function (e.g. the radial function of hydrogen) using NDSolve and would now like to vary a parameter such that the function converges to zero at infinity.
I am thinking of writing an iterative process that imposes convergence conditions.
Can Mathematica do this for me in an easier way? Or, is it maybe even possible to impose convergence to my function in NDSolve?

Here's an example:
The function is bound when a=-1/18

eps = $MachineEpsilon;
a = -0.055;
func = r D[r R[r], {r, 2}] + (2 (r + a r^2) - 2) R[r] == 0;
nrad = NDSolve[{func, R[eps] == 0, R'[eps] == 1}, 
  R[r], {r, eps, 100}]
Plot[R[r] /. nrad, {r, 0, 100}, PlotRange -> {{0, 50}, {-0.15, 0.3}}]


  • Prev by Date: Re: new functional operator
  • Next by Date: Re: new functional operator
  • Previous by thread: Re: Using Throw/Catch in MeshFunctions doesn't work? (8.0.4)
  • Next by thread: notebook interface problem -- erratic line wrapping