Метод Гаусса, вывести все этапы преобразования!

[PuHrBuH]

Помогите вывести все этапы преобразования матрицы, в поле StringGrid. Так же на форме должна быть кнопка, при нажатии на которую, отображается следующий этап преобразования.

unit gauss;

interface

uses
  Windows, Messages, SysUtils, Variants, Classes, Graphics, Controls, Forms,
  Dialogs, StdCtrls, Grids, ExtCtrls;

  const Max = 10;

type
    Vector = array[1..Max] of Double;
    Matrix = array[1..Max] of Vector;
    TForm1 = class(TForm)
    StringGrid1: TStringGrid;
    StringGrid2: TStringGrid;
    LabeledEdit1: TLabeledEdit;
    Button1: TButton;
    Button2: TButton;
    Memo1: TMemo;
    Button3: TButton;
    Panel1: TPanel;
    procedure LabeledEdit1Change(Sender: TObject);
    procedure Button3Click(Sender: TObject);
    procedure Button1Click(Sender: TObject);
  private
    { Private declarations }
  public
    { Public declarations }
  end;

var
  Form1: TForm1;

implementation

{$R *.dfm}


procedure TForm1.LabeledEdit1Change(Sender: TObject);
begin
  with StringGrid1, labelededit1 do
  begin
    ColCount := StrToIntDef(Text, 3);
    RowCount := StrToIntDef(Text, 3);
  end;
  with StringGrid2, labelededit1 do
    RowCount := StrToIntDef(Text, 3)
end;

procedure TForm1.Button3Click(Sender: TObject);
begin
 memo1.Clear;
end;

procedure TForm1.Button1Click(Sender: TObject);
  var a: Matrix;
      b,x: Vector;
      h: Double;
      i,j,k,n:integer;
begin
  n := StrToIntDef(Text, StringGrid1.ColCount);
  for j := 0 to n - 1 do
    for i := 0 to n - 1 do
      a[i + 1, j + 1] := StrToFloatDef(StringGrid1.Cells[j, i], 0);
  for I := 0 to n - 1 do b[i + 1] := StrToFloatDef(StringGrid2.Cells[0, i], 0);
  //Прямой ход
  for i:=1 to n-1 do
    for j:=i+1 to n do
    begin
      a[j,i]:=-a[j,i]/a[i,i];
      for k:=i+1 to n do
        a[j,k]:=a[j,k]+a[j,i]*a[i,k];
        b[j]:=b[j]+a[j,i]*b[i]
    end;
    x[n]:=b[n]/a[n,n];
        //Обратный ход
    for i:=n-1 downto 1 do
    begin
      h:=b[i];
      for j:=i+1 to n do h:=h-x[j]*a[i,j];
      x[i]:=h/a[i,i]
    end;
Memo1.Lines.Add('Решение:');
Memo1.Lines.Add('--------------------------------------------------------');
    for i:=1 to n do
    memo1.lines.add('x(' + IntToStr(i) + ')=' + FloatToStr(x[i]));
Memo1.Lines.Add('--------------------------------------------------------');
end;

end.

[PuHrBuH]

HELP me please...