如何在ODE中插入列表：

Question

我正在尝试以数值方式解决以下ODE

有边界条件

y(x=0)=0;y(x=20)=1

其中f（x）是我已经数值求解的另一个ODE的解。 我正在尝试使用https://docs.scipy.org/doc/scipy-0.18.1/reference/generation/scipy.integrate.solve_bvp.html解决我的ODE，但是我遇到的问题是f（x）被写为清单

 coordinates x = np.arange(0,20.01,0.05)

而且我不知道如何在scipy.integrate.solve_bvp中插入列表以获得ODE解决方案。

提前致谢

Answer 1

最好的方法是求解耦合系统，即将DE与f集成在一起。 如果这是您先前在math.SE上遇到的问题的变体，那么这种耦合方法是不可避免的。

至于文字解，请使用插值法

def f(x): return numpy.interp(x, x_samples, f_samples)

然后在ODE函数中调用此函数

def y_ODE(x,y): return [y[1], -0.5*f(x)*y[1]]