如何求解具有边界条件的微分方程，其中之一是使用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

这个动作应该是什么取决于问题。希望这会有所帮助。