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

Const nmax=10;
Type mas=array[1..nmax] of real;
     matr=array[1..nmax] of mas;
Var A:matr;
    b,x:mas;
    i,j,n:integer;
Procedure Gaus(A:matr;B:mas;Var X:mas);
Var i,j,k,imax:integer;
    temp,pr,sum:real;
Begin
    For k:=1 to n-1 do
    Begin
     imax:=k;
     For i:=k+1 to n do
      if abs(A[i][k])>abs(A[imax][k]) then imax:=i;
     if imax<>k then
     Begin
       For j:=1 to n do
       Begin
         temp:=A[k][j];
         A[k][j]:=A[imax][j];
         A[imax][j]:=temp;
       End;
       temp:=B[k];
       B[k]:=B[imax];
       B[imax]:=temp;
     End;
     For i:=k+1 to n do
     Begin
        pr:=a[i][k]/a[k][k];
        For j:=1 to n do
         a[i][j]:=a[i][j]-a[k][j]*pr;
        b[i]:=b[i]-b[k]*pr;
     End;
     End;
     x[n]:=b[n]/a[n][n];
     For k:=n-1 downto 1 do
     Begin
        sum:=0;
        For i:=k+1 to n do
         sum:=sum+x[i]*a[k][i];
        x[k]:=(b[k]-sum)/a[k][k];
     End;
End;
Begin
    Repeat
     write('Размерность матрицы 1<=N<=',nmax,': ');readln(n);
    Until n in [1..nmax];
    writeln('Матрица A');
    For i:=1 to n do
    Begin
      For j:=1 to n do
      Begin
        A[i][j]:=random(10)+1;
        write(A[i][j]:6:2);
        B[i]:=B[i]+A[i][j];
      End;
      writeln(' | ',B[i]:6:2);
    End;
    Gaus(A,B,X);
    write('X = (',X[1]:6:2);
    For i:=2 to n do
     write(';',X[i]:6:2);
    writeln('  )');
    readln;
End.

Const nmax=10;
Type mas=array[1..nmax] of real;
     matr=array[1..nmax] of mas;
Var A:matr;
    b,x:mas;
    i,j,n:integer;
    flag:boolean;
Procedure Gaus(A:matr;B:mas;Var X:mas;Var flag:boolean);
Var i,j,k,imax:integer;
    temp,pr,sum:real;
Begin
    flag:=true;
    For k:=1 to n-1 do
    Begin
     imax:=k;
     For i:=k+1 to n do
      if abs(A[i][k])>abs(A[imax][k]) then imax:=i;
     if imax<>k then
     Begin
       For j:=1 to n do
       Begin
         temp:=A[k][j];
         A[k][j]:=A[imax][j];
         A[imax][j]:=temp;
       End;
       temp:=B[k];
       B[k]:=B[imax];
       B[imax]:=temp;
     End;
     For i:=k+1 to n do
     Begin
        pr:=a[i][k]/a[k][k];
        For j:=1 to n do
         a[i][j]:=a[i][j]-a[k][j]*pr;
        b[i]:=b[i]-b[k]*pr;
        if a[i][i]=0 then
        Begin
          flag:=false;
          Exit;
        End;
     End;
     End;
     x[n]:=b[n]/a[n][n];
     For k:=n-1 downto 1 do
     Begin
        sum:=0;
        For i:=k+1 to n do
         sum:=sum+x[i]*a[k][i];
        x[k]:=(b[k]-sum)/a[k][k];
     End;
End;
Begin
    Repeat
     write('Введите размерность 1<=N<=',nmax,': ');readln(n);
    Until n in [1..nmax];
    For i:=1 to n do
     For j:=1 to n do
     Begin
       write('A[',i,',',j,'] = ');
       readln(A[i][j]);
     End;
    For i:=1 to n do
    Begin
       write('B[',i,'] = ');
       readln(B[i]);
    End;
    writeln('Ìàòðèöà A');
    For i:=1 to n do
    Begin
      For j:=1 to n do
       write(A[i][j]:6:2);
      writeln(' | ',B[i]:6:2);
    End;
    Gaus(A,B,X,flag);
    if flag then
    Begin
      write('X = (',X[1]:6:2);
      For i:=2 to n do
       write(';',X[i]:6:2);
      writeln('  )');
    End
    else writeln('Система вырождена');
    readln;
End.