Пошаговый вывод метода Гаусса.

[Leshii]

Собственно нужна помощь, в пошаговом выводе метода Гаусса.
хотя бы до размера 5х5. Без корней, просто по порядку. Я сейчас не очень хорошо соображаю, прошу помощи. Заранее благодарен.

uses crt;
const n=5
var a:array[1..n,1..n] of real;
    b:array[1..n] of real;
    n,i,j,k:integer;z,r,g:real;
begin
clrscr;
writeln('Введите коэффициенты системы и свободные члены');
for i:=1 to n do
for j:=1 to n do
read(A[i,j]);
readln;
clrscr;
writeln('Матрица коэффициентов и свободных членов');
for i:=1 to n do begin
 for j:=1 to n do
    write(A[i,j]:6:2);
    writeln;
end;
readln;
for i:=1 to n do
b[i]:=a[i,n];
for k:=1 to n do begin
       for j:=k+1 to n do
          begin
            r:=a[j,k]/a[k,k];
            for i:=k to n do
               begin
                 a[j,i]:=a[j,i]-r*a[k,i];
               end;
            b[j]:=b[j]-r*b[k];
          end;
     end;
for i:=1 to n do begin
 for j:=1 to n do
     write(A[i,j]:6:2);
     writeln;
end;
readln;
end.