About FindRoot?
- To: mathgroup at smc.vnet.net
- Subject: [mg62503] About FindRoot?
- From: stevenfield <boxiong at 126.com>
- Date: Sat, 26 Nov 2005 02:47:00 -0500 (EST)
- Sender: owner-wri-mathgroup at wolfram.com
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