delphi数值算法集光盘源代码,用于delphi7科学出版社出版.

				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;