在python中求解微分方程组

问题描述

我正在尝试用 python 求解一个微分方程组。 我有一个由两个方程组成的系统，其中有两个变量 A 和 B。 初始条件是A0=1e17 和B0=0，它们同时变化。 我使用 ODEINT 编写了以下代码：

import numpy as np
from scipy.integrate import odeint

def dmdt(m,t):
    A,B = m

    dAdt = A-B
    dBdt = (A-B)*A

    return [dAdt,dBdt]

# Create time domain
t = np.linspace(0,100,1)

# Initial condition
A0=1e17
B0=0

m0=[A0,B0]

solution = odeint(dmdt,m0,t)

显然我得到的输出与预期的不同，但我不明白错误。 有人能帮我吗？ 谢谢

解决方法

暂无找到可以解决该程序问题的有效方法，小编努力寻找整理中！

如果你已经找到好的解决方法，欢迎将解决方案带上本链接一起发送给小编。

小编邮箱:dio#foxmail.com (将#修改为@）