Исправить код. Тема: умножение матриц методом Штрассена

Const Max=500;    //верхний предел для размерности матриц
type ms{,aa,bb,cc:}= array [1..Max,1..Max] of Longint;
Var  
  aa,bb,cc: ms;
  ii,jj,i,j: Integer;   //нумерация по матрицам
  n: String;   // размерность матриц
  Num: Longint; Code,nn,nnn: Integer;    //для проверки введённог значения размерности матриц
  //ms: array [1..Max,1..Max] of Shortint;
procedure rand(var aa:ms; nn:integer);
var
  i,j: Integer;
begin
  for i:=1 to nn do
  for j:=1 to nn do begin
    aa[i,j]:= 50-random(100);
    end;
end;
 
procedure form(aa,bb:ms; nnn:integer; var cc:ms; nn:integer);
var
 p1,p2,p3,p4,p5,p6,p7,p8,p9,p10,p11,p12,p13,p14: ms;
 i1,i2,j1,j2: Integer;
begin   //накосячил с див
 if (nn>1) then begin
   i1:=1; i2:=(nn div 2)+1;
   while (i1<=nn div 2) and (i2<=nn) do begin
     j1:=1; j2:=(nn div 2)+1;
     while (j1<=nn div 2) and (j2<=nn) do begin2
       p1[i1,j1]:=aa[i1,j1]+aa[i2,j2];
       p2[i1,j1]:=bb[i1,j1]+bb[i2,j2];
       p3[i1,j1]:=aa[i2,j1]+aa[i2,j2];
       p4[i1,j1]:=bb[i1,j1];
       p5[i1,j1]:=aa[i1,j1];
       p6[i1,j1]:=bb[i1,j2]-bb[i2,j2];
       p7[i1,j1]:=aa[i2,j2];
       p8[i1,j1]:=bb[i2,j1]-bb[i1,j1];
       p9[i1,j1]:=aa[i1,j1]+aa[i1,j2];
       p10[i1,j1]:=bb[i2,j2];
       p11[i1,j1]:=aa[i2,j1]-aa[i1,j1];
       p12[i1,j1]:=bb[i1,j1]+bb[i1,j2];
       p13[i1,j1]:=aa[i1,j2]-aa[i2,j2];
       p14[i1,j1]:=bb[i2,j1]+bb[i2,j2];
       j1:=j1+1; j2:=j2+1;
     end;
     i1:=i1+1; i2:=i2+1;
    end;
   form(p1,p2,nnn,cc,nn div 2);
   form(p3,p4,nnn,cc,nn div 2);
   form(p5,p6,nnn,cc,nn div 2);
   form(p7,p8,nnn,cc,nn div 2);
   form(p9,p10,nnn,cc,nn div 2);
   form(p11,p12,nnn,cc,nn div 2);
   form(p13,p14,nnn,cc,nn div 2);
  end;
 i:=1; j:=1;
 while i<=nnn do
   begin
     j:=1;
     while j<=nnn do
       begin
       cc[i,j]:=cc[i,j]+p1[i,j]*p2[i,j]+p7[i,j]*p8[i,j]-p9[i,j]*p10[i,j]+p13[i,j]*p14[i,j];
       cc[i,j+1]:=cc[i,j+1]+p5[i,j]*p6[i,j]+p9[i,j]*p10[i,j];
       cc[i+1,j]:=cc[i+1,j]+p3[i,j]*p4[i,j]+p7[i,j]*p8[i,j];
       cc[i+1,j+1]:=cc[i+1,j+1]+p1[i,j]*p2[i,j]+p5[i,j]*p6[i,j]-p3[i,j]*p4[i,j]+p11[i,j]*p12[i,j];
       j:=j+2;
       end;
     i:=i+2;
   end;
end;
Begin
Randomize;
  Writeln('Введите размерность для двух квадратных матриц: ');
  repeat
      ReadLn(n);
  Val(n, Num, Code);
    if Code <> 0 then //проверка на число
    Write('Введено не числовое значение. перепроверьте: ');
  until Code = 0;
  nn:= StrToInt(n);   //преобразование введённой размерности после проверки, в целочисленное значение
  nnn:=nn;
  if (nn mod 2 = 0) then begin
    {if (nn=68) then begin 
    rand(aa,nn); rand(bb,nn);
    form(aa,bb,nnn,cc,nn);
    for i:=1 to nn do begin
    for j:=1 to nn do write(cc[i,j]:5); writeln; end;
    end else begin}
    for ii:= 1 to nn do    //запрос у пользователя каждый элемент для матрицы aa*
    for jj:= 1 to nn do begin    //циклит до момента пока не буду введены все элементы от одного до N*
        write('Введите A[',ii,',',jj,']:  ');
        readln(aa[ii,jj]);
      end;
      
    writeln('Введите B матрицу :  ');
    for ii:= 1 to nn do    //запрос у пользователя каждый элемент для матрицы bb*
    for jj:= 1 to nn do begin    //циклит до момента пока не буду введены все элементы от одного до N*
        write('Введите B[',ii,',',jj,']:  '); 
          readln(bb[ii,jj]);
      end;
 
form(aa,bb,nnn,cc,nn);
writeln('Результат умнажения матриц:  ');
  for ii:=1 to nn do begin
       for jj:= 1 to nn do
          write(cc[ii,jj]:5);
          writeln;
      end;
    end
  else WriteLn('Введённое значение не является кратным размерности степени два.');
end.