相关代码 |
|
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;