Table of Contents
Numerical integration is a fundamental technique in engineering for approximating the value of integrals that cannot be solved analytically. SciPy, a Python library, provides powerful tools to perform these calculations efficiently. This article presents step-by-step examples to help engineers understand how to use SciPy for numerical integration.
Basic Numerical Integration with SciPy
The most common function for numerical integration in SciPy is quad. It computes the definite integral of a function over a specified interval.
First, import the necessary module and define the function to integrate.
Example:
import scipy.integrate as integrate
def func(x):
return x**2
Then, perform the integration over the interval [0, 1].
result, error = integrate.quad(func, 0, 1)
The variable result contains the approximate value of the integral, and error estimates the error.
Integrating Multiple Functions
SciPy can handle more complex functions and multiple intervals. For example, integrating a sine function over [0, π].
Define the function:
import numpy as np
def sine_func(x):
return np.sin(x)
Calculate the integral:
area, error = integrate.quad(sine_func, 0, np.pi)
The result approximates the area under the sine curve between 0 and π.
Using Simpson’s Rule with SciPy
SciPy also provides simpson for Simpson’s rule, which is useful for equally spaced data points.
Example with data points:
x = np.linspace(0, 1, 5)
y = x**2
Calculate the integral:
area = integrate.simpson(y, x)
This method is suitable when data points are sampled at regular intervals.