第一题.矩阵法-梯形法积分(共6页).doc

梯形法数值积分A算法说明： 梯形法数值积分采用的梯形公式是最简单的数值积分公式，函数在区间a,b上计算梯形法数值积分表达式为： 由于用梯形公式来求积分十分粗糙，误差也比较大，后来改进后提出了复合梯形公式：,其中，n为积分区间划分的个数；h为积分步长。在MATLAB中编程实现的复合梯形公式的函数为：Combine Traprl.功能：复合梯形公式求函数的数值积分。调用格式：I,step=CombineTraprl(f,a,b,eps).其中，f为函数名； a为积分下限； b为积分上限； eps为积分精度； I为积分值； Step为积分划分的区间个数B流程图开始abs ( I2 I1 ) eps n=n+1; h=(b-a) / n; I1=I2; I2=0; for i=0:n-1 x=a+h*i; x1=x+h;I2=I2+( h/2 ) * ( subs (