《数值分析复化梯形公式,复化Simpson公式MATLAB程序》由会员分享,可在线阅读,更多相关《数值分析复化梯形公式,复化Simpson公式MATLAB程序(1页珍藏版)》请在人人文库网上搜索。
1、分别用复化梯形公式、复化Simpson公式计算定积分,取n=2,4,8,16分别验证结果(精确值I=4.006994)。复化梯形公式求定积分:function I=tquad(x,y)%复化梯形求积公式,其中,%x为向量,被积函数自变量的等距结点;%y为向量,被积函数在结点处的函数值;n=length(x);m=length(y);%积分自变量的结点数应与它的函数值的个数相同h=(x(n)-x(1)/(n-1);a=1 2*ones(1,n-2) 1;I=h/2*sum(a.*y);复化Simpson公式求定积分:function I=squad(x,y)%复化Simpson求积公式,其中,%x为向量,被积函数自变量的等距结点;%y为向量,被积函数在结点处的函数值;n=length(x);m=length(y);%积分自变量的结点数应与它的函数值的个数相同if rem(n-1,2)=0I=tquad(x,y);return;endN=(n-1)/2;h=(x(n)-x(1)/N;a=zeros(1,n);for k=1:Na(2*k-1)=a(2*k-1)+1;a(2*k)=a(2*k)+4;a(2*k+1)=a(2*k+1)+1;endI=h/6*sum(a.*y。