Re: About FindRoot?
- To: mathgroup at smc.vnet.net
- Subject: [mg62548] Re: About FindRoot?
- From: Peter Pein <petsie at dordos.net>
- Date: Sun, 27 Nov 2005 02:41:42 -0500 (EST)
- References: <dm9549$6nl$1@smc.vnet.net>
- Sender: owner-wri-mathgroup at wolfram.com
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