Calculating Integrals in Matlab

MATLAB provides powerful tools for calculating integrals, both numerically and symbolically. This guide explains how to perform integration in MATLAB, including basic methods, advanced techniques, and practical examples.

Introduction to Integrals in MATLAB

Integration is a fundamental operation in calculus that finds the area under a curve or the accumulation of quantities. MATLAB offers several functions to perform integration, including numerical methods like integral and symbolic computation with the Symbolic Math Toolbox.

The definite integral of a function f(x) from a to b is:

∫[a,b] f(x) dx = F(b) - F(a)

where F(x) is the antiderivative of f(x).

MATLAB's integral function performs numerical integration using adaptive quadrature, which is efficient and accurate for many functions. The Symbolic Math Toolbox provides symbolic integration capabilities for exact results.

Basic Integration Methods

Numerical Integration with integral

The integral function computes the definite integral of a function over a specified interval. Here's a basic example:

Syntax: integral(fun, a, b)

where fun is a function handle, and a and b are the integration limits.

Example: Calculate the integral of sin(x) from 0 to π.

∫[0,π] sin(x) dx = 2

Symbolic Integration with int

If you have the Symbolic Math Toolbox, you can use the int function for symbolic integration. This provides exact results when possible.

Syntax: int(fun, x)

for indefinite integration, or int(fun, a, b) for definite integration.

Example: Find the antiderivative of x².

∫ x² dx = (1/3)x³ + C

Advanced Integration Techniques

Multiple Integrals

MATLAB can handle multiple integrals using nested integral calls or the Symbolic Math Toolbox.

Example: Calculate the double integral of e^(-x²-y²) over the unit square.

∫∫[0,1] e^(-x²-y²) dx dy ≈ 0.7468

Improper Integrals

For improper integrals with infinite limits, MATLAB's integral function can handle these cases by specifying Inf or -Inf as limits.

Example: Calculate the integral of 1/(1+x²) from 0 to ∞.

∫[0,∞] 1/(1+x²) dx = π/2

Practical Examples

Example 1: Area Under a Curve

Calculate the area under the curve of f(x) = x² from 0 to 1.

∫[0,1] x² dx = 1/3

Example 2: Work Done by a Variable Force

Calculate the work done by a force F(x) = x³ from x=0 to x=2.

Work = ∫[0,2] x³ dx = 4

Example 3: Probability Density Function

Calculate the probability that a normally distributed random variable X with mean 0 and standard deviation 1 falls between -1 and 1.

P(-1 ≤ X ≤ 1) = ∫[-1,1] (1/√(2π))e^(-x²/2) dx ≈ 0.6827

Frequently Asked Questions

What is the difference between integral and int?: The integral function performs numerical integration, while int performs symbolic integration. Use integral for numerical results and int when you need exact symbolic results.
How do I handle integration limits that are not finite?: For improper integrals, use Inf or -Inf as the integration limits in the integral function. MATLAB will handle the infinite limits appropriately.
Can I integrate functions with parameters?: Yes, you can integrate functions with parameters by using anonymous functions or function handles. For example, integral(@(x) x.^2 + a, 0, 1) integrates x² + a from 0 to 1.
What if my function is not continuous?: If your function has discontinuities, the integral function may still work, but you should check the results carefully. For functions with singularities, consider using symbolic integration or adjusting the integration limits.
How accurate are the numerical integration results?: The integral function uses adaptive quadrature, which provides high accuracy for most well-behaved functions. The default absolute and relative error tolerances are typically sufficient, but you can adjust them if needed.