Заменить прямой метод Гаусса

const n=2;
type
Matr = array[1..n,1..n] of real;
Vect = array[1..n] of real;
var a:Matr;
    b,x,d,b2:Vect;
    a2:Matr;
    i,j,m,k,z,p:integer;
    max,r,g:real;
begin
     randomize;
     for i:=1 to n do
         begin
              for j:=1 to n do
              a[i,j]:=random;
              b[i]:=random;
         end;
         a2:=a;
         b2:=b;
         for k:=1 to n do
             begin {Прямой ход}
                  b[k]:=b[k]/a[k,k];
                  for j:=n downto k do a[k,j]:=a[k,j]/a[k,k];
                  for j:=k+1 to n do
                      begin
                           r:=a[j,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 i:=1 to n do
         begin
         write('b=',b[i] );
         writeln;
         for j:=1 to n do begin
         write(a[i,j]:2:1,'  ');
 
         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('Korni:');
         for i:=1 to n do
             write('x[',i,']=',x[i]:3:2,' ');
         writeln;
         {погрешность}
         writeln;
         for i:=1 to n do
             begin
                  d[i]:=0;
             end;
         for i:=1 to n do
             begin
                  for j:=1 to n do
                      begin
                           d[i]:=d[i]+a2[i,j]*x[j];
                      end;
             end;
         max:=abs(d[1]-b2[1]);
         m:=1;
         for m:=2 to n do
             for m:=2 to n do
                 begin
                      if abs(d[i]-b2[i])>max then
                         begin
                              max:=abs(d[i]-b2[i]);
                              m:=i;
                         end;
                 end;
         writeln(max,' ',m);
end.

for k:=1 to n do
 begin {Прямой ход}
   b[k]:=b[k]/a[k,k];
   for j:=n downto k do a[k,j]:=a[k,j]/a[k,k];
   for j:=k+1 to n do
    begin
     for i:=k to n do a[j,i]:=a[j,i]-a[j,k]*a[k,i];
     b[j]:=b[j]-r*b[k];
    end;
 end;