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

可以设置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

有谁知道如何设置这种条件?

参考方案

如果您提供了有关该问题的更多信息,它将更容易回答。
如果我理解正确,L是积分的极限,所以边界条件实际上是

Psi(L)=eps

如果积分超出L,则可能需要有一个条件才能检查每次迭代(新点等)。大致情况:

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

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

Python sqlite3数据库已锁定 - python

我在Windows上使用Python 3和sqlite3。我正在开发一个使用数据库存储联系人的小型应用程序。我注意到,如果应用程序被强制关闭(通过错误或通过任务管理器结束),则会收到sqlite3错误(sqlite3.OperationalError:数据库已锁定)。我想这是因为在应用程序关闭之前,我没有正确关闭数据库连接。我已经试过了: connectio…

Python pytz时区函数返回的时区为9分钟 - python

由于某些原因,我无法从以下代码中找出原因:>>> from pytz import timezone >>> timezone('America/Chicago') 我得到:<DstTzInfo 'America/Chicago' LMT-1 day, 18:09:00 STD…

用大写字母拆分字符串,但忽略AAA Python Regex - python

我的正则表达式:vendor = "MyNameIsJoe. I'mWorkerInAAAinc." ven = re.split(r'(?<=[a-z])[A-Z]|[A-Z](?=[a-z])', vendor) 以大写字母分割字符串,例如:'我的名字是乔。 I'mWorkerInAAAinc”变成…

如何打印浮点数的全精度[Python] - python

我编写了以下函数,其中传递了x,y的值:def check(x, y): print(type(x)) print(type(y)) print(x) print(y) if x == y: print "Yes" 现在当我打电话check(1.00000000000000001, 1.0000000000000002)它正在打印:<…

在Mac上的终端中停止python - python

在Mac上的终端中使用python,键入ctrl-z 将停止python,但不退出它,给出如下输出:>>> [34]+ Stopped python 如您所见,我已经停止了34个python调用。虽然我可以用>>> exit() 退出python,问题是:是否有一个快捷键可以真正在终端中退出(而不只是停止)python?而…