Решить систему алгебраических уравнений методом Гаусса

uses crt;
const
m=3;
n=4;
var d,b,db,c,dbc:array[1..n,1..n]of integer;
    i,j,p,x:integer;
begin
clrscr;
randomize;
writeln('Matrica D:');
for i:=1 to m do
   begin
     for j:=1 to n do
       begin
         d[i,j]:=random(10)-4;
         write(d[i,j]:4);
       end;
     writeln;
   end;
writeln('Matrica B:');
for i:=1 to m do
   begin
     for j:=1 to n do
       begin
         b[i,j]:=random(10)-4;
         write(b[i,j]:4);
       end;
     writeln;
   end;
{сложение}
writeln('Матрица DB=D+B:');
for i:=1 to m do
   begin
     for j:=1 to n do
       begin
         db[i,j]:=d[i,j]+b[i,j];
         write(db[i,j]:4);
       end;
     writeln;
   end;
readln;
writeln('Matrica C:');
for i:=1 to n do
   begin
     for j:=1 to m do
       begin
         c[i,j]:=random(10)-4;
         write(c[i,j]:4);
       end;
     writeln;
   end;
{умножение}
for p:=1 to m do
for j:=1 to m do
   begin
     dbc[p,j]:=0;
     for i:=1 to n do
       begin
         dbc[p,j]:=dbc[p,j]+db[p,i]*c[i,j];
       end;
   end;
writeln('Matrica DBC=(D+B)*C:');
for i:=1 to m do
    begin
      for j:=1 to m do
      write(dbc[i,j]:5);
      writeln;
    end;
readln;
{транспонирование}
for i:=1 to m do
for j:=1 to m do
if j<=i then
 begin
  x:=dbc[i,j];
  dbc[i,j]:=dbc[j,i];
  dbc[j,i]:=x;
 end;
writeln('Matrica DBC_T:');
for i:=1 to m do
    begin
      for j:=1 to m do
      write(dbc[i,j]:5);
      writeln;
    end;
readln;
{добавление единичной матрицы}
writeln('Rezultat:');
for i:=1 to m do
    begin
      for j:=1 to m do
       begin
        if i=j then dbc[i,j]:=dbc[i,j]+1;
        write(dbc[i,j]:5);
       end;
      writeln;
    end;
readln
end.