Проверка для метода Гаусса

//Метод Гаусса
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,  m1, k1: integer;
  z, r, g, m, s: real;
 
begin
  clrscr;
  n := 4;
  
  A[1, 1] := -5;A[1, 2] := 2;A[1, 3] := 7;A[1, 4] := -5;
  A[2, 1] := 1;A[2, 2] := 1;A[2, 3] := -1;A[2, 4] := 9;
  A[3, 1] := 3;A[3, 2] := 5;A[3, 3] := 4;A[3, 4] := 6;
  A[4, 1] := -3;A[4, 2] := 2;A[4, 3] := 4;A[4, 4] := 7;
  b[1] := -4;b[2] := 5;b[3] := 4;b[4] := 6;
  
  writeln('Матрица A: ');
  for i := 1 to n do
  begin
    for j := 1 to n do
    begin
      write(a[i, j]:5:0);
      a1[i, j] := a[i, j];
    end;
    writeln;
  end;
  writeln;
  Writeln('Матрица B: ');
  for i := 1 to n do
    writeln(b[i]);
  writeln;
  {Выбор главного элемента по диагонали}
  
  for k := 1 to n do
  begin
    k1 := k + 1;
    m := (A[k, k]);
    m1 := k;
    for i := k1 to n do 
    begin
      r := a[i, k];
      if abs(r) > abs(m) then begin
        m := r;
        m1 := i;
      end;
    end;
 
  
    for j := 1 to n do 
    begin
      z := A[k, j];
      A[k, j] := A[m1, j];
      A[m1, j] := z;
    end;
    
  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:5, '   ');
  
end.