Помогите написать программу для Вычисления обратной матрицы методом Гаусса.
Вот что смогла, обратный ход не умею
Program obr_matr;
const n=4;
var A,B: array [1..n,1..n] of real;
i,j,k,m: integer;
c: Real;
function Max(k: integer):integer;
var m:integer;
i: integer;
begin
m:=k;
for i:=k+1 to n do begin
if abs(A[i,k])>abs(A[m,k]) then m:=i;
end;
result:=m;
end;
procedure Perestanovka(x,y: integer);
var l:real;
j:integer;
begin
for j:= 1 to n do begin
l:=A[x,j];
A[x,j]:=A[y,j];
A[y,j]:=l;
end;
end;
begin
randomize;
for i:=1 to n do begin
for j:=1 to n do begin
A[i,j]:=random(-5,5);
end;
end;
writeln('Матрица A:');
for i:=1 to n do begin
for j:=1 to n do begin
write(A[i,j]:6:3 );
end;
writeln;
end;
writeln;
for k:=1 to n-1 do begin
m:=max(k);
if a[m,k]=0 then continue;
if m<>k then perestanovka (j,k);
for i:=k+1 to n do begin
c:=a[i,k]/a[k,k];
for j:=k to n do
a[i,j]:=a[i,j]-a[k,j]*c;
end;
end;
writeln('после преобразования к треугольному виду');
for i:=1 to n do begin
for j:=1 to n do begin
write(a[i,j]:6:3,' ');
end;
writeln;
end;