Python: определенный интеграл

Задача №2: Расчитать значение определенного интегралла известной функции на диапазоне [a,b]

Суть задачи сводится к нахождению площади фигуры ограниченной подынтегральной функцией и осью абцисс. Проиллюстрирован метод гвадратурных формул Гаусса и Симпсона.

import scipy.integrate as spint
def f(x):
    return x*x
a = 0
b = 1
print spint.quad(f,a,b)

Если функция задана таблично, то нужно использовать другой метод, например Симпсона:

import scipy.integrate as spint
import numpy as np
x =  np.arange(0,1,0.1)
y = x**2 #массив значений функции в точках, для примера он посчитан, но можно и задать вручную
print spint.simps(f,a,b)

Добавить комментарий

Ваш e-mail не будет опубликован. Обязательные поля помечены *