unit2.pas

来自「1、说明： 本书中所有的常用数值算法子过程按书中的章数分别放在以C开头」· PAS 代码 · 共 40 行

unit Unit2;

interface
uses
  unit1;
procedure EULSUM(var SUM:real; TERM:real; JTERM:integer);

implementation
procedure EULSUM(var SUM:real;TERM:real;JTERM:integer);
var
    J:Integer;  TMP,DUM:real;
    WKSP:array[0..41] of real;
begin
    If JTERM = 1 Then
    begin
        NTERM^:=1;
        WKSP[1]:=TERM;
        SUM:=0.5 * TERM;
    end
    Else
    begin
        TMP:=WKSP[1];
        WKSP[1]:=TERM;
        For J:=1 To NTERM^ do
        begin
            DUM:=WKSP[J + 1];
            WKSP[J + 1]:=0.5 * (WKSP[J] + TMP);
            TMP:=DUM;
        end;
        If Abs(WKSP[NTERM^ + 1]) <= Abs(WKSP[NTERM^]) Then
        begin
            SUM:=SUM + 0.5 * WKSP[NTERM^ + 1];
            NTERM^:=NTERM^ + 1;
        end
        Else
            SUM:=SUM + WKSP[NTERM^ + 1];
    end;
end;
end.