在3D和2D库仑势的情况下，我需要数值计算径向薛定谔方程的特征值和本征函数。微分方程包含R（r）的一阶和二阶导数。
问题是scipy.integrate.ode或scipy.integrate.odeint需要＆＃34;初始值＆＃34;对于函数及其一阶导数 - 我在网上看到的每个例子都使用了t中的微分方程，因此指定系统的初始条件是微不足道的，是任意的。
由于我在考虑太空中的ODE而不是时间，我应该注意到＆＃34;初始条件＆＃34;应该更准确地称为边界条件＆＃34;。然而，对于我的情况，函数R（r）必须仅在原点处是有限的，但是在r = 0时它没有特定的值。此外，径向波函数的导数R＆＃39;在物理上是无意义的，因此约束（或特别指定）其在原点的值是无意义的。系统中唯一明确的边界条件是对于非常大的r，函数必须指数衰减到零。
在这种情况下，我正在考虑，不是让linspace从零计数到一些大数，而是让它从一个大值倒数到零，为此我可以设置我的＆＃34;初始条件为r - ＆GT;无穷＆＃34;这样R和R＆＃39;是零。
其他人有这个问题，你找到了哪些解决方法？
谢谢！