chebpc.txt

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

procedure CHEBPC(C:array of real;var D:array of real; N:integer);
var
    DD:array[0..50] of real;  J,K:integer;
    SV:real;
begin
    For J:=1 To N do
    begin
        D[J]:=0 ;
        DD[J]:=0; 
    end;
    D[1]:=C[N];
    For J:=N - 1 DownTo 2 do
    begin
        For K:=N - J + 1 DownTo 2 do
        begin
            SV:=D[K];
            D[K]:=2* D[K - 1] - DD[K];
            DD[K]:=SV;
        end;
        SV:=D[1];
        D[1]:=-DD[1] + C[J];
        DD[1]:=SV;
    end;
    For J:=N DownTo 2 do
        D[J]:=D[J - 1] - DD[J];
    D[1]:=-DD[1] + 0.5 * C[1];
end;