matlab 梯形逐次分半法的代码怎么写呀，可以在我的代码的基础上改一下

%vari_Tn
format long
a=input('请输入a：');
b=input('请输入b：');
err=input('请输入err：');
n=1;
h=(b-a)/2;
T0=h*(fun1(a)+fun1(b));
F=0;
for i=1:n
    F=F+fun1(a+(2*i-1)*h);
    T=T0/2+h*F;
    if(T-T0<3*err) 
    break;
    else
    T0=T;
    n=2*n;
    h=h/2;
    end;  
end;