کد برنامه گاوس جردن

خودم خیلی تو اینترنت دنبال کدش گشتم اما نبود آخرشم مجبور شدم خودم کدشو بنویسم
کد برنامه حل سه معادله ی سه مجهولی به روش گاوس جردن به زبان پاسکال.
می دونم این زبان اصلاً به درد نمی خوره اما خوب هنوز خیلی دانشجوها لازمه برنامه نویسی رو با این زبان پاس کنن.

Program GJ;

uses crt;

var mat : array [1..3 , 1..4] of real;

    mata : array [1..4] of real;

    i : integer;

    a : real:

 begin

  

 clrscr;

 for i := 1 to 4 do

  begin

  write('zarib 'i' moadeleye 1 ra vared konid: ');

  readln(mat[1,i]);

  end;

for i := 1 to 4 do

  begin

  write('zarib 'i' moadeleye 2 ra vared konid: ');

  readln(mat[2,i]);

  end;

for i := 1 to 4 do

  begin

  write('zarib 'i' moadeleye 3 ra vared konid: ');

  readln(mat[3,i]);

  end;

  a := mat[1,1];

for i := 1 to 4 do

  begin

  mat[1,i] := ( mat[1,i] / a );

  end;

  a := ( mat[2,1] * -1 );

for i := 1 to 4 do

  begin

  mata[i] := ( mat[1,i] * a );

  mat[2,i] := ( mat[2,i] + mata[i] );

  end;

 a := ( mat[3,1] * -1 );

for i := 1 to 4 do

  begin

  mata[i] := ( mat[1,i] * a );

  mat[3,i] := ( mat[3,i] + mata[i] );

  end;

 a := mat[2,2];

for i := 1 to 4 do

  begin

  mat[2,i] := ( mat[2,i] / a );

  end;

a := ( mat[1,2] * -1 );

for i := 1 to 4 do

  begin

  mata[i] := ( mat[2,i] * a );

  mat[1,i] := ( mat[1,i] + mata[i] );

  end;

a := ( mat[3,2] * -1 );

for i := 1 to 4 do

  begin

  mata[i] := ( mat[2,i] * a );

  mat[3,i] := ( mat[3,i] + mata[i] );

  end;

 a := mat[3,3];

for i := 1 to 4 do

  begin

  mat[3,i] := ( mat[3,i] / a );

  end;

a := ( mat[1,3] * -1 );

for i := 1 to 4 do

  begin

  mata[i] := ( mat[3,i] * a );

  mat[1,i] := ( mat[1,i] + mata[i] );

  end;

a := ( mat[2,3] * -1 );

for i := 1 to 4 do

  begin

  mata[i] := ( mat[3,i] * a );

  mat[2,i] := ( mat[2,i] + mata[i] );

  end;

writeln ('javab haye masale :'mat[1,4]','mat[2,4]','mat[3,4]' ast.');

readln;

end.