СЛАУ методом Гаусса

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.

uses  crt;
const n=3;
var a:array[1..n,1..n] of real;
    b,x:array[1..n] of real;
    i,j,k:integer;
    z,r,g:real;
begin
clrscr;
a[1,1]:=0.62; a[1,2]:=0.56; a[1,3]:=-0.43;
a[2,1]:=1.32; a[2,2]:=-0.88; a[2,3]:=1.76;
a[3,1]:=0.73; a[3,2]:=1.43; a[3,3]:=-0.34;
b[1]:=1.16; b[2]:=2.07; b[3]:=2.18;
writeln('Расширенная матрица');
for i := 1 to n do
 begin
  for j := 1 to n do
  write(a[i,j]:7:2);
  writeln(b[i]:9:2);
 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
    a[j,i]:=a[j,i]-r*a[k,i];
    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
writeln('x[', i, ']=', x[i]:0:2);
readln
end.