Python odeint用法

Question

谁能告诉我如何将odeint函数的参数放入python中，以获得积分。

我试图得到e^x从0到2的积分，但我不确定如何输入参数。文档对我来说不是很清楚。

谢谢

Answer 1

如果你只想计算积分，而不是解微分方程，你也可以使用

from numpy import exp
from scipy.integrate import quad

def f(x):
    return exp(x)

result, error = quad(f, 0, 2)

Answer 2

下面是一个假设y(0) = 1的示例。

import scipy
import scipy.integrate


def integrateExp(y0, a, b):
  limits = [a, b]
  integral = scipy.integrate.odeint(lambda y, t : scipy.exp(t), y0, limits)
  return integral[1]

print integrateExp(1, 0, 2)

第一个参数应为(y，t)并返回相应的导数。我在这里使用了lambda，因为d/dt (T)是微不足道的。