Использование Python odeint

Question

Может кто-нибудь сказать мне, как вставить параметры для функции odeint для питона, чтобы получить интеграл.

Я пытаюсь получить интеграл от e ^ x от 0 до 2, но я не уверен, как ввести параметры. Документация мне не совсем понятна.

Спасибо

Answer 1

Вот пример, в котором предполагается, что 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) и возвращать соответствующую производную. Я использовал лямбду здесь, поскольку d / dt exp (t) тривиально.

Answer 2

Если вы хотите просто вычислить интегралы, а не решать дифференциальные уравнения, вы также можете использовать

from numpy import exp
from scipy.integrate import quad

def f(x):
    return exp(x)

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