python - 动态积分方程系统

我有一个用于odeint Python模块的ODE系统：

def equations(p, t, lmbds):
    return np.array([-p[0] * lmbds[0] + p[1] * lmbds[1],
                     -p[1] * lmbds[1] - p[1] * lmbds[2] + p[0] * lmbds[0]])

该系统对应于具有两个状态，三个强度的马尔可夫链。

方程组可以包含不同数量的方程。能否根据用户输入的方程式数量随时生成它们？

最佳答案

请参阅以下创建n个方程的方法：

def get_equation(p, t, lmbds, n):
    """ Returns equation for provided n. """
    # Create and return equation

def equations(p, t, lmbds, eq_num):
    equations = []
    for n in range(eq_num):
        equation = get_equation(p, t, lmbds, n)
        equations.append(equation)
    return np.array(*equations)

关于python - 动态积分方程系统，我们在Stack Overflow上找到一个类似的问题：https://stackoverflow.com/questions/51896674/