qtrap.txt

插值 1.拉格朗日插值 2.有理函数插值 3.三次样条插值 4.有序表的检索法 5.插值多项式 6.二元拉格朗日插值 7.双三次样条插值

procedure QTRAP(A, B:real; var S:real);
const
    EPS = 0.000001;    JMAX = 20;
var
    OLDS:real;  J:integer;
begin
    OLDS:= -1E+30;
    For J:= 1 To JMAX do
    begin
        TRAPZD(A, B, S, J);
        If (Abs(S - OLDS) < (EPS * Abs(OLDS))) Then Exit;
        OLDS:= S;
    end;
    ShowMessage('Too many steps.');
end;