我刚刚开始学习计算机科学以将其应用于物理学，我正在尝试编写一个代码来解决谐波振荡器的薛定谔方程（设置$V=\frac{x^2}{2}$) 一维。

我知道如何做到这一点的基础知识，但我想知道是否有一些资源可用，或者是否有人可以帮助我如何编写我的程序（我对 Python 更有信心，但 C、C++ 或任何其他语言也可以正常工作）。

我知道我应该：

1. 选择最大角动量的值$\ell$和量子数$n$.
2. 寻求能量区间$E_i, E_{i-1}$我可以在哪里更准确地寻找解决方案。
3. 通过算法找到能量的最佳近似值（例如：割线公式）。
4. 以一定的精度比较能量间隔$\epsilon$成立之初。

预先感谢您的帮助！