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

Question

可以设置像 psi(z>L)< eps 这样的边界条件吗？

我有一个系统:

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

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