eigsrt.txt

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

Procedure EIGSRT(var D:array of real; var V:matrx2; N:integer);
var
    I,K,J:integer;  P:real;
begin
    For I:=1 To N - 1 do
    begin
        K:=I;
        P:=D[I];
        For J:=I + 1 To N do
        begin
            If D[J] >= P Then
            begin
                K:=J;
                P:=D[J];
            end ;
        end; 
        If K <> I Then
        begin
            D[K]:=D[I];
            D[I]:=P;
            For J:=1 To N do
            begin
                P:=V[J, I];
                V[J, I]:=V[J, K];
                V[J, K]:=P;
            end; 
        end;
    end; 
end;