MathGroup Archive 2005

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

Search the Archive

Re: About FindRoot?


stevenfield schrieb:
> Why not work of following notes? Thanks in advance!
> --------------------------
> hbar=1;m=1;V[x]=(1/2)x^2;stepEnergy[n_] = -5(1 - n/10);
> Sol[Energy_] := NDSolve[{-hbar^2/(2m)D[Ï?'[x], x] + V[x]Ï?[x]== Energy Ï?[x],Ï?[-15] == 0,Ï?'[-15] == 1},Ï?, {x, -15, 15}][[1]];
> Ï?end[Energy_] := Ï?end[Energy] = (Ï?[15] /. Sol[Energy]);
> Energyshooting[n_] := Energy /. FindRoot[Ï?end[Energy]==0, {Energy, 
>     1.1 stepEnergy[n], stepEnergy[n]}][[1]];
> Energyshooting[5]
> ----------------------------------
> stevenfield
> 

Hi Steven,

I uploaded a notebook (~55kB) to 
http://people.freenet.de/Peter_Berlin/Mathe/stevenfield/Deq.nb .
It has almost no comments but I hope it will nevertheless be helpful.

Peter


  • Prev by Date: Re: Saving Points Selected from a ListPlot
  • Next by Date: Re: Re: How to View Mathematica and Hardcopy Books
  • Previous by thread: About FindRoot?
  • Next by thread: Making a phase plot