splin2.txt

来自「《Delphi常用数值算法集》的配书源码」· 文本 代码 · 共 18 行

procedure SPLIN2(X1A, X2A:array of real; YA, Y2A:matrx2;
                              M, N:integer;var X1, X2, Y:real);
var
    YTMP,Y2TMP,YYTMP:array[0..100] of real;
    J,K:integer;
begin
    For J:= 1 To M do
    begin
        For K:= 1 To N do
        begin
            YTMP[K]:= YA[J, K];
            Y2TMP[K]:= Y2A[J, K];
        end;
        SPLINT(X2A, YTMP, Y2TMP, N, X2, YYTMP[J]);
    end;
    SPLINE(X1A, YYTMP, M, 1E+30, 1E+30, Y2TMP);
    SPLINT(X1A, YYTMP, Y2TMP, M, X1, Y);
end;