dttixing.m

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

%递推法运用梯形法求得函数得积分,此函数求得所要满足条件的值可能所需的步骤为多.在需要高精度是可能较为慢,一般求
%低精度的积分,f1是所需进行积分的函数,a,b是积分区间,a为上限,b为下限,e是精度值,是求得值与精确值的差小于e,这只是
%充分条件,而不是必要条件
function dttixing(f1,a,b,e)
n=1;
f=inline(f1);
h=(b-a);
T=h/2*(f(a)+f(b));
S=1/2*(T+h/2*f((a+b)/2));
k=1;
fprintf('%d%13.8f\n',k,S);
while abs(T-S)>e
    
    T=S;
    n=2*n;
    h=h/2;
    s1=0;
    for i=0:n-1
        s1=s1+f(a+0.5*h+i*h);
    end
    S=1/2*(T+h*s1);
    k=k+1;
    fprintf('%d%13.8f\n',k,S);
end