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

uses
  crt;
 
const
  e = 0.000001;
 
var
  a,a1: array[1..10, 1..10] of real;
  b,b1: array[1..10] of real;
  x: array[1..10] of real; 
  n, i, j, k: integer;z, r, g: real;
 
begin
  clrscr;
  write('Порядок: ');
  readln(n);
  writeln('Введите коэффициенты системы и свободные члены');
  for i := 1 to n do
  begin
    for j := 1 to n do
    begin
      write('a[', i, ',', j, ']= ');
      readln(a[i, j]);
      a1[i,j]:=a[i,j];
    end;
    Writeln('b[', i, ']= ');
    readln(b[i]);
  end;
  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 k := n downto 1 do { обратный ход Гаусса }
  begin
    r := 0;
    for j := k + 1 to n do
    begin
      g := a[k, j] * x[j];
      r := r + g;
    end;
    x[k] := (b[k] - r) / a[k, k];
  end;
  writeln('Корни системы:');
  for i := 1 to n do
    write('x[', i, ']=', x[i]:0:2, '   ');
  writeln;
  writeln('Проверка в соответствии с матрицей B: ');
 
  for i:=1 to n do
   begin 
    for j:=1 to n do  
      b1[i]:=b1[i]+a1[i,j]*x[j];     
     writeln(b1[i]:3:2,' ');
    end; 
  readln;
end.

uses
  crt;
 
const
  e = 0.000001; { точность вычислений }
var
 
  a,a1: array[1..10, 1..10] of real;
  b,b1: array[1..10] of real;
  x: array[1..10] of real;
  n, i, j, k: integer;z, r, g: real;
 
begin
  clrscr;
  n:=3;
 
  A[1, 1] := 1.54; A[1, 2] := 1.70; A[1, 3] := 1.62;
  A[2, 1] := 3.69; A[2, 2] := 3.73; A[2, 3] := 3.59;
  A[3, 1] := 2.45; A[3, 2] := 2.43; A[3, 3] := 2.25;
  b[1] := -1.97; b[2] := -3.74; b[3] := -2.26;
 
  writeln('Матрица A: ');
  for i := 1 to n do
  begin
    for j := 1 to n do
    begin
      write(a[i,j]:5:2);
      a1[i,j]:=a[i,j];
    end;
     writeln;
     end;
   writeln;
    Writeln('Матрица B: ');
   for i:=1 to n do
    writeln(b[i]);
  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 k := n downto 1 do { обратный ход Гаусса }
  begin
    r := 0;
    for j := k + 1 to n do
    begin
      g := a[k, j] * x[j];
      r := r + g;
    end;
    x[k] := (b[k] - r) / a[k, k];
  end;
  writeln('Корни системы:');
  for i := 1 to n do
    write('x[', i, ']=', x[i]:0:2, '   ');
end.