simpson.m

为配套王能超版本数值分析简明教程编写,方便学生教师演示!

%此程序是运用辛甫生公式求的函数的在区间中的积分,它是预先要输入一个在区间中n等分,其实程序是把区间分成了2n等分,
%因为进行辛甫生法是只能在区间2n等分进行的,此算法能较为精确的求得函数得积分值,f1是所需进行积分的函数,a,b是积分
%区间,a为上限,b为下限,n是把区间n等分(实际是2n等分)
function s=Simpson(f1,a,b,n)
f=inline(f1);
h=(b-a)/(2*n);
s=h/3*(f(a)+f(b));
s1=0;
s2=0;
for i=1:n-1
    s1=s1+2*h/3*f(a+2*i*h);
    s2=s2+4*h/3*f(a+(2*i-1)*h);
end
s2=s2+4*h/3*f(a+(2*n-1)*h);
s=s+s1+s2;
fprintf('%13.8f',s);