Accelerating the pace of engineering and science

# Documentation Center

• Trials

### Contents

Numerically evaluate triple integral

## Description

q = triplequad(fun,xmin,xmax,ymin,ymax,zmin,zmax) evaluates the triple integral fun(x,y,z) over the three dimensional rectangular region xmin <= x <= xmax, ymin <= y <= ymax, zmin <= z <= zmax. The first input, fun, is a function handle. fun(x,y,z) must accept a vector x and scalars y and z, and return a vector of values of the integrand.

Parameterizing Functions explains how to provide additional parameters to the function fun, if necessary.

q = triplequad(fun,xmin,xmax,ymin,ymax,zmin,zmax,tol) uses a tolerance tol instead of the default, which is 1.0e-6.

## Examples

Pass function handle @integrnd to triplequad:P

`Q = triplequad(@integrnd,0,pi,0,1,-1,1);`

where the file integrnd.m is

```function f = integrnd(x,y,z)
f = y*sin(x)+z*cos(x);```

Pass anonymous function handle F to triplequad:

```F = @(x,y,z)y*sin(x)+z*cos(x);

This example integrates y*sin(x)+z*cos(x) over the region 0 <= x <= pi, 0 <= y <= 1, -1 <= z <= 1. Note that the integrand can be evaluated with a vector x and scalars y and z.