chder.txt

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

procedure CHDER(A,B:real;C:array of real;var CDER:array of real;N:integer);
var
    J:integer;    CON:real;
begin
    CDER[N]:=0; 
    CDER[N - 1]:=2 * (N - 1) * C[N];
    If N >= 3 Then
    begin
        For J:=N - 2 DownTo 1 do
            CDER[J]:=CDER[J + 2] + 2 * J * C[J + 1];
    end;
    CON:=2/(B - A);
    For J:=1 To N do
        CDER[J]:=CDER[J] * CON;
end;