trapz (MATLAB Functions)

Trapezoidal numerical integration

Syntax

Z = trapz(Y)
Z = trapz(X,Y)
Z = trapz(...,dim)

Description

Z = trapz(Y) computes an approximation of the integral of Y via the trapezoidal method (with unit spacing). To compute the integral for spacing other than one, multiply Z by the spacing increment. Input Y can be complex.

If Y is a vector, trapz(Y) is the integral of Y.

If Y is a matrix,trapz(Y) is a row vector with the integral over each column.

If Y is a multidimensional array, trapz(Y) works across the first nonsingleton dimension.

Z = trapz(X,Y) computes the integral of Y with respect to X using trapezoidal integration. Inputs X and Y can be complex.

If X is a column vector and Y an array whose first nonsingleton dimension is length(X), trapz(X,Y) operates across this dimension.

Z = trapz(...,dim) integrates across the dimension of Y specified by scalar dim. The length of X, if given, must be the same as size(Y,dim).

Examples

Example 1

The exact value of

is 2.

To approximate this numerically on a uniformly spaced grid, use

```
X = 0:pi/100:pi;
Y = sin(x);
```

Then both

```
Z = trapz(X,Y)
```

and

```
Z = pi/100*trapz(Y)
```

produce

```
Z =
    1.9998
```

Example 2

A nonuniformly spaced example is generated by

X = sort(rand(1,101)*pi);
Y = sin(X);
Z = trapz(X,Y);

The result is not as accurate as the uniformly spaced grid. One random sample produced

```
Z =
    1.9984
```

Example 3

This example uses two complex inputs:

z = exp(1i*pi*(0:100)/100);

trapz(z, 1./z)
ans =
   0.0000 + 3.1411i

See Also

cumsum, cumtrapz

trace treelayout