如何用边界条件求解微分方程，其中之一是使用python的不等式

Question

可以设置psi（z> L）我有一个系统：

dD/dz = ....  
dPsi/dz = ....   

边界条件：

D(0) = 1,6*Pi  
Psi(z>L) < eps

我了解了scipy.integrate.solve_bvp和odeint的知识，但是所有示例都具有如下边界条件：y(0) = y(1) = 0

有人知道如何设置这种条件吗？

Answer 1

如果您提供了有关该问题的更多信息，则回答起来会更容易。 如果我理解正确，则L是积分的极限，因此边界条件实际上是

Psi(L)=eps

如果积分超过L，则可能需要有一个条件才能检查每次迭代（新点等）。大致情况：

if z>L and Psi(z)>=eps:
    do the action

此操作应采取何种措施取决于问题。希望这会有所帮助。