Решение уравнения

unit Unit1;
 
interface
 
uses
  Windows, Messages, SysUtils, Variants, Classes, Graphics, Controls, Forms,
  Dialogs, StdCtrls, Grids;
 
type
  TForm1 = class(TForm)
    Button1: TButton;
    StringGrid1: TStringGrid;
    Label1: TLabel;
    Label2: TLabel;
    StringGrid2: TStringGrid;
    Label3: TLabel;
    StringGrid3: TStringGrid;
    procedure Button1Click(Sender: TObject);
  private
    { Private declarations }
  public
    { Public declarations }
  end;
const n=3;
var
  Form1: TForm1;
 
implementation
 
{$R *.dfm}
 
procedure TForm1.Button1Click(Sender: TObject);
var i,j,k:byte;
    a,a1:array[1..n,1..n] of real;
    b,b1,x,y,m:array[1..n] of real;
    h,g:real;
begin
//исходная матрица коэффициентов
 
for i:=1 to n do
for j:=1 to n do
a[i,j]:=Random(100);
a1:=a; //запомним ее для проверки, первуая изменится
//то же свободные члены
 
for j:=1 to n do
b[j]:=Random(100);
b1:=b;
for i:=1 to n do
for j:=1 to n do
StringGrid1.Cells[j-1,i-1]:=FloattostrF(a[i,j],ffFixed,0,2);
for i:=0 to n-1 do
 begin
  StringGrid2.Cells[0,i]:='X'+inttostr(i+1)+'=';
  StringGrid1.Cells[n,i]:=FloattostrF(b[i+1],ffFixed,0,2);
 end;
{приведение матрицы к треугольному виду}
for k:=1 to n do
 begin
  for j:=k+1 to n do
   begin
    h:=a[j,k]/a[k,k];
    for i:=k to n do
    a[j,i]:=a[j,i]-h*a[k,i];
    b[j]:=b[j]-h*b[k];
   end;
 end;
{вычисление корней}
for k:=n downto 1 do
 begin
  h:=0;
  for j:=k+1 to n do
   begin
    g:=a[k,j]*x[j];
    h:=h+g;
   end;
  x[k]:=(b[k]-h)/a[k,k];
 end;
 
for i:=1 to n do
StringGrid2.Cells[1,i-1]:=FloattostrF(x[i],ffFixed,0,2);
//проверка
for i:=1 to n do
 begin
  y[i]:=0;
  for j:=1 to n do
  y[i]:=y[i]+a1[i,j]* x[j];
 end;
for i:=1 to n do
m[i]:=abs(b1[i]-y[i]); //невязка, у меня все по нолям
for i:=1 to n do
StringGrid3.Cells[0,i-1]:=FloattostrF(m[i],ffFixed,0,6);
end;
 
end.

  a[i,j]:=Random(100);
...
  b[j]:=Random(100);

unit Unit1;
 
interface
 
uses
  Windows, Messages, SysUtils, Variants, Classes, Graphics, Controls, Forms,
  Dialogs, StdCtrls, Grids;
 
type
  TForm1 = class(TForm)
    StringGrid1: TStringGrid;
    StringGrid2: TStringGrid;
    StringGrid3: TStringGrid;
    Label1: TLabel;
    Label2: TLabel;
    Label3: TLabel;
    Button1: TButton;
    procedure Button1Click(Sender: TObject);
  private
    { Private declarations }
  public
    { Public declarations }
  end;
 
var
  Form1: TForm1;
 
const
  n = 3;
 
implementation
 
{$R *.dfm}
 
procedure TForm1.Button1Click(Sender: TObject);
var
  I, J, K : Integer;
  a, a1    : array[1..n, 1..n] of Single;
  b, b1,
  x, y, m  : array[1..n] of Single;
  h, g     : Single;
begin
  //исходная матрица коэффициентов
  for I := 1 to n do
    for J := 1 to n do
      a[I, J] := Random(100) + 1;
 
  a1 := a;                              // запомним ее для проверки, первуая изменится
                                        // то же свободные члены
 
  for J := 1 to n do
    b[J] := Random(100) + 1;
 
  b1 := b;
 
  for I := 1 to n do
    for J := 1 to n do
      StringGrid1.Cells[J - 1, I - 1] := Format('%.2f', [a[I, J]]);
 
  for I := 0 to n - 1 do
  begin
    StringGrid2.Cells[0, I] := 'X' + IntToStr(I + 1) + '=';
    StringGrid2.Cells[1, I] := Format('%.2f', [b[I + 1]]);
  end;
 
  { приведение матрицы к треугольному виду }
  for K := 1 to n do
  begin
    for J := k + 1 to n do
    begin
      h := a[J, K] / a[K, K];
      for I := K to n do
        a[J, I] := a[J, I] - H * a[K, I];
      b[J] := b[J] - H * b[K];
    end;
  end;
 
  { вычисление корней }
  for K := n downto 1 do
  begin
    H := 0;
    for J := k + 1 to n do
    begin
      g := a[K, J] * x[J];
      h := h + g;
    end;
    x[K] := (b[K] - H) / a[K, K];
  end;
 
  for I := 1 to n do
    StringGrid2.Cells[1, I - 1] := Format('%.2f', [x[I]]);
 
  // проверка
  for I := 1 to n do
  begin
    y[I] := 0;
    for J := 1 to n do
      y[I] := y[I] + a1[I, J] * x[J];
  end;
 
  for I := 1 to n do
    m[I] := abs(b1[I] - y[I]);                  // невязка, у меня все по нолям
 
  for I := 1 to n do
    StringGrid3.Cells[0, I - 1] := Format('%.6f', [m[I]]);
end;
 
end.

StringGrid1.RowCount:=n;
StringGrid1.ColCount:=n;