tixing.m

来自「基于Matlab的程序」· M 代码 · 共 20 行

function [I,n]=Tixing(a,b,errors)
h=b-a;
T=h*(f(a)+f(b))/2;
n=0;
while 1
    H=0;
    x=a+h/2;
    while x<b
        H=H+f(x);
        x=x+h;
    end
    T1=(T+h*H)/2;
    if abs(T1-T)<errors
        I=T1;
        break;
    end
    h=h/2;
    n=n+1;
    T=T1;
end