2016春计算方法实验指导-作业

I=trapz(x,y)

x和y是同维向量或矩阵，返回积分结果． 3. 复化Simpson公式数值积分quad

quad(fun,a,b,tol)

quad采用自适应步长的Simpson求积法，是低阶法数值积分的函数．它自动变换、选择步长，以满足精度要求，实现变步长复合抛物形积分计算．其中，fun为积分函数，可用字符表达式，内联函数或M文件的函数形式标识，[a,b]为积分区间，tol 为积分绝对误差，默认为1E-6． 4. quadl函数

quadl(fun,a,b,tol)

采用自适应递推步长复合Lobatto 数值积分法计算积分． 5. dblquad函数

I=dblquad(fun,a,b,c,d,tol)

在矩形区域上求二重积分，fun为二元函数，可用inline定义或写成M文件函数的形式，[a,b]为变量x的上下限，[c,d]是变量y的上下限，tol 为精度要求，缺省值为1E-6． 6. triplequad函数

I=triplequad(fun,a,b,c,d,e,f,tol)

用于三重积分在立体区域上计算三重积分结果，fun为三元函数，可用内联函数定义或写成M文件函数形式，[a,b]是变量x的积分上下限，[c,d]是变量y的积分上下限，[e,f]为变量z的积分上下限，tol为积分精度，缺省值为1E-6． 7. 复化的8阶Newton-Cotes公式quad8

I=quad8(fun,a,b,tol) 参数意义同上．

6.3 实验作业

实验名称：数值积分．

实验目的：学会利用Matlab实现数值积分. 实验内容：用不同数值方法计算积分

10?4xlnxdx??.

9实验要求：取不同的步长h，分别用复合梯形公式及复合辛普森公式计算积分，给出误差中关于h的函数，并与积分精确值比较两个公式的精度.

49