Анализ гауссовской квадратурной функции

Question

Попытка заставить Python решить уравнение f (x) = ((2pi) ^ (- 1/2)) * e ^ ((- x ^ (2)) / 2), используя квадратуру Гаусса из x = -С 4 по 4.

Попытка выяснить, как сказать Python интегрироваться относительно линии, пересекающей график при x = -4 и 4.

Не знаю, как вообще начать решатьэта проблема, кодов ошибок пока нет.

Answer 1

Вы можете использовать библиотеку Scipy. субпакет scipy.integrate включает гауссову квадратуру , которую вы можете использовать.

from scipy import integrate
from numpy import exp,pi

f = lambda x: ((2*pi)**(-1/2))*exp(-0.5*x**2)
print(integrate.quadrature(f, -4.0, 4.0))