How can I solve a non homogeneous diff equation numerically?
显示 更早的评论
I'm trying to solve a quantum harmonic oscillator numerically and I followed Matlab steps to crack the equation, however, it doesn't work. Any other idea for another method or where the mistake could be? the equation is d2psi/dx2=(constant-x^2)*psi Thanks for advanced
syms y(x)
k=23; m=9.1093821545*1e-31; w=sqrt(k/m);hb=(6.6260695729*1e-34)/(2*pi);E0=hb*w*.5;
[V] = odeToVectorField(diff(y, 2) == (hb*w*.5-.5*k*x.^2).*y);
x1=-sqrt(hb*w/k);
x2=-x1;
M = matlabFunction(V,'vars',{x,'Y'});
[m,n] = ode45(M,[x1 x2],[0. 0.]);
采纳的回答
Torsten
2017-11-27
k=23;
m=9.1093821545*1e-31;
w=sqrt(k/m);
hb=(6.6260695729*1e-34)/(2*pi);
E0=hb*w*.5;
M=@(x,y)[y(2);(hb*w*.5-.5*k*x^2)*y(1)];
x1=-sqrt(hb*w/k);
x2=-x1;
[m,n] = ode45(M,[x1 x2],[0.; 0.]);
Result will be that y(1) and y(2) are identically zero.
Best wishes
Torsten.
2 个评论
david feldstein
2017-11-27
Torsten
2017-11-27
But this is the code to solve the equation
d2y/dx2=(constant-x.^2)*y
y(1) will be y, y(2) will be dy/dx.
Best wishes
Torsten.
更多回答（0 个）
类别
在 帮助中心 和 File Exchange 中查找有关 Quantum Mechanics 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!
Translated by 
