Решение системы линейных уравнений методом Гаусса

var
  Form1: TForm1;
  N,M,i,j,k :integer;
 
  a,b,x,y :array of array of double;
 
  mn:double;
 
 
 
implementation
 
{$R *.dfm}
 
procedure TForm1.Button1Click(Sender: TObject);
 
begin
 if not tryStrToInt(Edit2.Text, M) or not tryStrToInt(Edit1.Text, N)
 then
 begin
 ShowMessage('Введите корректные данные.');
 exit
 end;
 
N:=StrToInt(Edit1.Text);
M:=StrToInt(Edit2.Text);
if (N<0) or (M<0) then
 begin
 ShowMessage('Количество строк или столбцов матриц не может быть равно отрицательному числу.');
 exit
 end;
 if (N = 0) or (N = 0) then
  begin
 ShowMessage('Количество строк или столбцов матриц не может быть равно 0.');
 exit
 end;
 
SetLength(A,N+1,M+1);
 
   with StringGrid1 do begin
   RowCount:=N+1; ColCount:=M+1;
   Cells[0,0]:=' N/N';
   for j:=1 to M do Cells[j,0]:=IntToStr(j);
    for i:=1 to N do Cells[0,i]:=IntToStr(i);
 
    for i:= 1 to N do
    for j:= 1 to M do begin
    Cells[j,i]:=FloatToStr(Random(100));
    A[i,j]:=StrToFloat(Cells[j,i]);
    end;
   end;
 
end;
procedure Perest;
var
Prom: double;
begin
 for i:= k+1 to N do
 begin
 if A[i,k] <> 0 then
   begin
   for j:= 1 to N do
   begin
    Prom:=A[k,j];
    A[k,j]:= A[i,j];
    A[i,j]:=Prom;
   end;
 
  end;
  end;
end;
procedure TForm1.Button2Click(Sender: TObject);
begin
 
  for k:= 1 to N-1 do begin
    for i:= k + 1 to N do begin
        if A[k,k] = 0 then begin
                           Perest;
                           end;
        Mn := A[i,k]/A[k,k];
        for j:= k to M do begin
 
           A[i,j] := A[i,j] - Mn * A[k,j];
        end;
    end;
  end;
 
  with StringGrid2 do begin
    RowCount:=N+1; ColCount:=M+1;
    for j:=1 to M do Cells[j,0]:=IntToStr(j);
    for i:=1 to N do Cells[0,i]:=IntToStr(i);
    for i:= 1 to N do
    for j:= 1 to M do Cells[j,i]:=FloatToStrF(A[i,j],ffFixed,5,1);
  end;
end;
 
 
procedure TForm1.StringGrid1KeyPress(Sender: TObject; var Key: Char);
begin
if not (key in ['0'..'9','.']) then
key:=#0;
end;
 
end.