I was wondering if there was any way to have mathematica find the solution to Schrodinger's equation [(−h^2/2m)(d^2ψ/dx^2) + kx^2ψ=Eψ] for a particle centered at the origin. I keep getting a Tag Plus Protected error when I try to enter and evaluate the equation.
If you type
eqn = (-h^2/2 m) D[\[Psi][x], {x, 2}] + k x^2 \[Psi][x] == e \[Psi][x]
DSolve[eqn, \[Psi][x], x]
Mathematica will return
\[Psi](x)->Subscript[c, 1] Subscript[D, (Sqrt[2] e-h Sqrt[k] Sqrt[m])/(2 h Sqrt[k] Sqrt[m])]
((2^(3/4) Power[k, (4)^-1] x)/(Sqrt[h] Power[m, (4)^-1]))+Subscript[c, 2]
Subscript[D, (-Sqrt[2] e-h Sqrt[k] Sqrt[m])/(2 h Sqrt[k] Sqrt[m])]((I 2^(3/4)
Power[k, (4)^-1] x)/(Sqrt[h] Power[m, (4)^-1]))
which is the solution keeping in mind that D stands for ParabolicCylinderD and Subscript[c, 1] and Subscript[c, 2] are constant of integration.
If you love us? You can donate to us via Paypal or buy me a coffee so we can maintain and grow! Thank you!
Donate Us With