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

uses crt;
const n=3;
      a:array[1..n,1..n] of real=((4,0.24,-0.08),
                                  (0.09,3,-0.15),
                                  (0.04,-0.08,4));
      b:array[1..n] of real=(8,9,20);
var x:array[1..3] of real;
    i,j,k:integer;z,r,g:real;
begin
clrscr;
{приведение матрицы к треугольному виду}
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
write('x[',i,']=',x[i]:0:2,'   ');
readln;
end.

uses crt;
const n=3;
var a1,b1,c1,d1,
    a2,b2,c2,d2,
    a3,b3,c3,d3,
    x1,x2,x3:real;
begin
clrscr;
{приведение матрицы к треугольному виду - прямой ход Гаусса}
a1:=4;b1:=0.24;c1:=-0.08;d1:=8;
a2:=0.09;b2:=3;c2:=-0.15;d2:=9;
a3:=0.04;b3:=-0.08;c3:=4;d3:=20;
{обнуляем коэффициенты при х1 во 2 и 3 уравнениях}
b2:=b1+b2*(-a1/a2); c2:=c1+c2*(-a1/a2); d2:=d1+d2*(-a1/a2);
b3:=b1+b3*(-a1/a3); c3:=c1+c3*(-a1/a3); d3:=d1+d3*(-a1/a3);
{обнуляем коэффциет при х2 в третьем уравнении}
c3:=c3+c2*(-b3/b2); d3:=d3+d2*(-b3/b2);
{ищем корни, начиная с последнего - обратный ход Гаусса}
x3:=d3/c3;
x2:=(d2-c2*x3)/b2;
x1:=(d1-b1*x2-c1*x3)/a1;
writeln('Корни системы:');
write('x1=',x1:0:2,'  x2=',x2:0:2,'  x3=',x3:0:2);
readln;
end.

program Gauss_Classic;
uses crt;
 
type vec=array[1..10] of real;
mas=array[1..10] of vec;
 
var a:mas;b,x:vec;
h:real;
i,j,k,n:integer;
 
Begin
Clrscr;
write('Vvedite kol-vo uravneniy: ');read(n);
for i:=1 to n do begin
writeln('Uravnenie # ',i);
writeln('Vvedi koefficienti');
for j:=1 to n do begin
write('a[',i,',',j,']=');read(a[i,j]);end;
writeln('Vvedi pravu chast uravneni');
write('b[',i,']=');read(b[i]); end;
Clrscr;
writeln('Ishodnaya matrica');
for i:=1 to n do begin
write(i);
for j:=1 to n do
write(a[i,j]:10:5);
writeln(b[i]:10:5) end;
writeln; writeln; writeln;
{Pryamoi hod - iskl peremenih}
for i:=1 to n-1 do
for j:=i+1 to n do begin
a[j,i]:=-a[j,i]/a[i,i];
for k:=i+1 to n do
a[j,k]:=a[j,k]+a[j,i]*a[i,k];
b[j]:=b[j]+a[j,i]*b[i] end;
x[n]:=b[n]/a[n,n];
{Obr hod  - nahozhdenie kornei}
for i:=n-1 downto 1 do begin
h:=b[i];
for j:=i+1 to n do h:=h-x[j]*a[i,j];
x[i]:=h/a[i,i]; end;
writeln('Korni uravneniya');
for i:=1 to n do writeln('x(',i,')=',x[i]:10:5);
readkey
end.