使用triplequad函数可以对三重积分进行数值计算。该函数的语法格式如下。
q=triplequad(fun,xmin,xmax,ymin,ymax,zmin,zmax)：调用quad函数对fun(,x,y)函数进行三重积分运算，积分区间为xmin≤x≤xmax，ymin≤y≤ymax，zmin≤z≤zmax。fun参数是一个M文件函数或匿名函数的句柄。fun(x,y,z)函数必须接受向量x和标量y和z，并返回积分值组成的向量。
q=triplequad(fun,xmin,xmax,ymin,ymax,zmin,zmax,tol)：用误差容限tol代替默认的误差容限1.0e-6。
q=triplequad(fun,xmin,xmax,ymin,ymax,zmin,zmax,tol,method)：用method参数指定积分方法，代替默认的quad函数。method还可以被指定为@quadl或用户指定的函数句柄。该函数必须与quad和quadl具有相同的语法格式。
【例27】用triplequad函数求下面三重积分的数值近似。
\[\int_{-1}^{1} {\int_{0}^{1} {\int_{0}^{\pi} {(y\operatorname{sin}{(}x)+z\operatorname{cos}{(}x))} d xdydz}}\]
首先，建立函数的M文件integrnd.m，代码为
code.matlab
function z = integrnd(x, y,z)
z = y*sin(x)+z*cos(x);然后，在命令窗口中输入
code.matlab
>> Q = triplequad(@integrnd,0,pi,0,1,-1,1);
Q=
2.0000