我一直在研究量子波动方程的径向求解器（薛定谔和狄拉克）。在这两种情况下，建议似乎都是采用“射击方法”，采用各种精度的集成方案，例如 Numerov 之一。

现在，对于一种拍摄方法，需要从边界条件强加起始值。许多文本提到了渐近行为，例如附近的贝塞尔函数$r=0$等等，因为我的兴趣是解决基本库仑势的扰动版本，并且这个问题的解决方案是已知的，但是，我不能从已知的类氢解决方案开始并从那里开始工作吗？或者有什么原因没有说明这一点？