dtsimpson.asv

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

%此函数是功能是求数的积分,此程序的设计思想是运用辛甫生公式而编成的,它能较为精确的求的函数在区间的积分
%在此程序中好处为运行步骤较为少就能得到要成立的精度,此精度的值是所求得的值与精确值的充分
function dtSimpson(f1,a,b,e)
f=inline(f1);
n=1;
h=(b-a)/2;
s2=0;
s1=f(a+h);
T=h/3*(f(a)+f(b)+4*s1);
s2=s2+s1;
s1=f(a+h/2)+f(a+3*h/2);
R=h/6*(f(a)+f(b)+2*s2+4*s1);
k=1;
fprintf('%d%13.8f\n',k,R);
while abs(R-T)>e
    n=2*n;
    h=h/2;
    s2=s2+s1;
    T=R;
    s3=0;
    for i=0:2*n-1
        s3=s3+f(a+h/2+h*i);
    end
    s1=s3;
    R=h/6*(f(a)+f(b)+4*s1+2*s2);
    k=k+1;
    fprintf('%d%13.8f\n',k,R);
end