chint.txt

《Delphi常用数值算法集》的配书源码

procedure CHINT(A, B:real;C:array of real;var CINT1:array of real;N:integer);
var
    SUM,FAC,CON:real;  J:integer;
begin
    CON:=0.25 * (B - A);
    Sum:=0; 
    FAC:=1; 
    For J:=2 To N - 1 do
    begin
        CINT1[J]:=CON * (C[J - 1] - C[J + 1]) / (J - 1);
        Sum:=Sum + FAC * CINT1[J];
        FAC:=-FAC;
    end;
    CINT1[N]:=CON * C[N - 1] / (N - 1);
    Sum:=Sum + FAC * CINT1[N];
    CINT1[1]:=2* Sum;
end;