доделать прогу

program pr53;
const nm=20;
var i,j,k,l,n:integer;
e,a,ais,ap:array [1..nm,1..nm] of real;
begin
writeln ('Input n');
readln (n);
writeln;
writeln('Inputed matrix');
randomize;
    for i:=1 to n do
    for j:=1 to n do
     ais[i,j]:=random(9)+1;
     a:=ais;
      for i:=1 to n do
 begin
 for j:=1 to n do
 write (ais[i,j]:2:0);
 writeln;
 end;
  readln;
 for k:=1 to n do
 begin
 for j:=1  to n do
 if (j<>k) then
 ap[k,j]:=-A[k,j]/a[k,k];
 for i:=1 to n do
 if (i<>k) then
 ap[i,k]:=a[i,k]/a[k,k];
 for i:=1 to n do
             for j:=1 to n do
            if (i<>k) and (i<>k) then
            ap[i,j]:=a[i,j]-a[i,k]*a[k,j]/a[k,k];
            ap[k,k]:=1/a[k,k];
            a:=ap;
            end;
 writeln ('obr matrix');
 for i:=1 to n do
 begin
 for j:=1 to n do
 write (ap[i,j]:2:0);
 writeln;
 end;
 readln;
 
  for i:=1 to n do
 begin
 for j:=1 to n do
 begin
 E[i,j]:=0;
 for l:=1 to n do
 E[i,j]:=E[i,j]+ap[i,l]*ais[l,j]
 end;
 writeln ('Delaem proverku');
 writeln;
 Writeln ('If A*A(obr)=E then wright');
 writeln;
 for i:=1 to n do
 begin
 for j:=1 to n do
 write (e[i,j]:2:0);
 writeln;
 end;
 readln;
 end;
 end.