6.1 数值积分基本原理求解定积分数值方法各种多样，如简单梯形法、辛普生(Simpson)法、牛顿－柯特斯(Newton-Cotes)法等都是经常采取方法。它们基本思想都是将整个积分区间[a, b]分成n个子区间[xi, xi+1]，i=1, 2, …, n，其中x1 = a，xn+1 = b。这么求定积分问题就分解为求和问题。
6.2 数值积分实现方法6.2.1 变步长辛普生法基于变步长辛普生法，MATLAB给出了quad函数来求定积分。该函数调用格式为：  [I, n] = quad('fname', a, b, tol, trace)其中fname是被积函数名。a和b分别是定积分下限和上限。tol用来控制积分精度，缺省时取tol=0.001。trace控制是否展现积分过程，若取非0则展现积分过程，取0则不展现，缺省时取trace=0。返回参数I即定积分值，n为被积函数调用次数。
例1 求定积分：(1) 建立被积函数文件fesin.m。function f=fesin(x)f=exp(-0.5*x).*sin(x+pi/6);(2) 调用数值积分函数quad求定积分。[S,n]= ...