Написать программу решения по методу Гаусса системы линейных уравнений

program slau;
 
{$APPTYPE CONSOLE}
 
uses SysUtils, Windows;
 
var i, j, k: integer;
    t: real;
    s: string;
    x: array [1..3] of real;
    a: array [1..3, 1..4] of real = ((4.0, 0.24, 0.08, 8.0),
                                     (0.09, 3.0, 0.15, 9.0),
                                     (0.04, 0.08, 4.0, 20.0));
 
procedure print(m: string);
begin
  writeln(m);
  for i := 1 to 3 do
    begin
      for j := 1 to 4 do write(a[i, j]:10:4);
      writeln
    end
end;
 
begin
  SetConsoleCp(1251);
  SetConsoleOutputCp(1251);
  print('Исходная матрица системы:');
  for k := 1 to 3 do
    begin
      for i := k + 1 to 3 do
        begin
          t := a[i, k] / a[k, k];
          for j := k to 4 do a[i, j] := a[i, j] - a[k, j] * t
        end
    end;
  print('Треугольная матрица системы:');
  for i := 1 to 3 do x[i] := 0;
  for i := 3 downto 1 do
  begin
    t := 0;
    for j := i to 3 do t := t + a[i, j] * x[j];
    x[i] := (a[i, 4] - t) / a[i, i];
  end;
  writeln('Корни системы:');
  for i := 1 to 3 do writeln('x[', i, '] = ', x[i]:10:4);
  readln
end.

unit Unit1;
 
interface
 
uses
Windows, Messages, SysUtils, Variants, Classes, Graphics, Controls, Forms,
Dialogs, StdCtrls, Grids;
 
type
TForm1 = class(TForm)
label1: TLabel;
label2: TLabel;
edit1: TEdit;
edit2: TEdit;
label3: TLabel;
button1: TButton;
button2: TButton;
stringrid1: TStringGrid;
procedure button2Click(Sender: TObject);
procedure button1Click(Sender: TObject);
private
{ Private declarations }
public
{ Public declarations }
end;
 
var
Form1: TForm1;
 
implementation
 
{$R *.dfm}
 
procedure TForm1.button2Click(Sender: TObject);
begin
Close;
end;
 
procedure TForm1.button1Click(Sender: TObject);
 
var a:array[1..100,1..100] of integer;
n,m,i,j,s:integer;
sr:real;
 
begin
n:=StrToInt(Edit1.Text);
m:=StrToInt(Edit2.Text);
stringrid1.RowCount:=n;
stringrid1.ColCount:=m;
s:=0;
randomize;
for i:=1 to n do
for j:= 1 to m do
begin
a[i,j]:=random(10)-random(10);
stringrid1.Cells[j-1,i-1]:=IntToStr(a[i,j]);
s:=s+a[i,j];
end;
sr:=s/(m*n);
label3.Caption:='sr='+FloatToStrF(sr,fffixed,4,2);
end;
end.