Now I found it but I have to zoom it. Actually I want to do a more complicated job. In that case I probably don't know what the scale Phi would be(I have to normalise Phi later). Do you have some idea to avoid problem like it?
Can't get Quantum Harmonic Oscillator Wave Function by ode45
1 view (last 30 days)
Show older comments
I have just tested if I can get correct Quantum Harmonic Oscillator Wave Function by ode45 But the result doesn't looks like what it should be. What's problem? (I haven't used ode45 for a while, maybe some basic silly mistakes)
function Wavefcn
[x,Phi] = ode45(@WaveFunction,[0 10],[1 0]);
plot(x,Phi(:,1));
function dPhi = WaveFunction(x,phi)
n = 0; %grand state
eps = 2*n+1;
dPhi = zeros(2,1);
dPhi(1) = phi(2);
dPhi(2) = (x^2-eps)*phi(1);
end
end
Answers (0)
See Also
Categories
Find more on Quantum Mechanics in Help Center and File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!