Simpson求积公式.docVIP

实验求积公式 【】1．掌握Simpson数值积分的基本原理； 2．编写Matlab对数值积分的实现程序【实验内容】 用复合Simpson公式求 ，误差不超过的近似值，已知定积分的准确值为-12.070346316。 【】【实验方法或步骤】 个节点，若被积函数依次在两个小区间上逐次利用Simpson法，则 若最后剩下三个小区间，就利用Simpson法，即： 则。 建立一个以l.m命名的m文件，程序如下： function I = Simps_v(f,h) n=length(f)-1; if n==1, fprintf(Data has only one interval),return; end if n==2, I = h/3*(f(1) + 4*f(2) + f(3)); return;end if n==3, I = 3/8*h*(f(1) + 3*f(2) + 3*f(3) + f(4)); return;end I=0; if 2*floor(n/2)~=n I = 3/8*h*(f(n -2) + 3*f(n -1) + 3*f(n) + f(n+1)); m=n -3 else m=n end I = I+ (h/3)*( f(1)+ 4*sum(f(2:2:m)) + f(m+1)); if m2, I = I