simp3.txt

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

Procedure SIMP3(var A:matrx2; MP, NP, I1, K1:integer;var IP, KP:integer);
var
    PIV:real;  II,KK:integer;
begin
    PIV:=1  / A[IP + 1, KP + 1];
    For II:=1 To I1 + 1 do
    begin
        If II - 1 <> IP Then
        begin
            A[II, KP + 1]:=A[II, KP + 1] * PIV;
            For KK:=1 To K1 + 1 do
            begin
                If KK - 1 <> KP Then
                    A[II, KK]:=A[II, KK] - A[IP + 1, KK] * A[II, KP + 1];
            end;
        end;
    end;
    For KK:=1 To K1 + 1 do
        If KK - 1 <> KP Then A[IP + 1, KK]:=-A[IP + 1, KK] * PIV;
    A[IP + 1, KP + 1]:=PIV;
end;