数值积分matlab程序.docVIP

第二章 数值积分2.1. 复化 Simpson 公式功能：利用复化 Simpson 公式计算被积函数 f(x)在给定区间上的积分值-----------------------------------------function S=FSimpson(f,a,b,n)% f:被积函数句柄% a,b:积分区间的两个端点% n:子区间个数% S:用复化 Simpson 法求得的积分值h=(b-a)/n;fa=feval(f,a);fb=feval(f,b);S=fa+fb;x=a;for i=1:Nx=x+h/2;fx=feval(f,x);S=S+4*fx;x=x+h/2;fx=feval(f,x);S=S+2*fx;end S=h*S/6;----------------------------------------------------------------------------------------------------------------------附：函数值为向量形式的 simpson 求积法function I=simpson_h(f,h)%调用格式 I＝simpson(f,h)%f 为一向量，指定已知节点处的函数值%h 为步长n=length(f)-1;if n==1fprintf(Data h