%遞推法運用梯形法求得函數得積分,此函數求得所要滿足條件的值可能所需的步驟為多.在需要高精度是可能較為慢,一般求
%低精度的積分,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