pcshft.txt

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

procedure PCSHFT(A, B:real;var D:array of real; N:integer);
var
    FAC,CONST1:real;  J,K:integer;
begin
    CONST1:=2/ (B - A);
    FAC:=CONST1;
    For J:=2 To N do
    begin
        D[J]:=D[J] * FAC;
        FAC:=FAC * CONST1;
    end;
    CONST1:=0.5 * (A + B);
    For J:=1 To N - 1 do
        For K:=N - 1 DownTo J do
            D[K]:=D[K] - CONST1 * D[K + 1];
end;