Умножение двух матриц представленных в показательной форме записи

uses crt;
type pok=record
         z,fi:real;
         end;
     alg=record
         re,im:real;
         end;
     mtr=array[1..10,1..10] of pok;
var n:byte;
 
function poktoalg(a:pok):alg;//перевод показательной в алгебраическую
var b:alg;
begin
b.re:=a.z*cos(a.fi);
b.im:=a.z*sin(a.fi);
poktoalg:=b;
end;
 
function algtopok(a:alg):pok;//перевод алгебраической а показательную
var b:pok;
begin
b.z:=sqrt(sqr(a.re)+sqr(a.im));
b.fi:=arctan(a.im/a.re);
algtopok:=b;
end;
 
function mult(a,b:alg):alg;
var c:alg;
begin
c.re:=a.re*b.re-a.im*b.im;
c.im:=a.re*b.im+b.re*a.im;
mult:=c;
end;
 
function add(a,b:alg):alg; 
var c:alg;
begin
c.re:=a.re+b.re;
c.im:=a.im+b.im;
add:=c;
end;
 
procedure mult_mtr(a,b:mtr;var c:mtr);
var a1,b1,c1:array[1..10,1..10] of alg;
    i,j,k:byte;
begin
for i:=1 to n do
for j:=1 to n do
 begin
  a1[i,j]:=poktoalg(a[i,j]);
  b1[i,j]:=poktoalg(b[i,j]);
 end;
for i:=1 to n do
for j:=1 to n do
 begin
  c1[i,j].re:=0;
  c1[i,j].im:=0;
  for k:=1 to n do
  c1[i,j]:=add(c1[i,j],mult(a1[i,k],b1[k,j]));
 end;
for i:=1 to n do
for j:=1 to n do
c[i,j]:=algtopok(c1[i,j]);
end;
 
procedure vvod(var a:mtr; c:char);
var i,j:byte;
begin
writeln('Введите элементы матрицы ',c);
for i:=1 to n do
for j:=1 to n do
 begin
  write('Введите модуль числа ',c,'[',i,',',j,'].z=');
  readln(a[i,j].z);
  write('Введите аргумент числа ',c,'[',i,',',j,'].fi=');
  readln(a[i,j].fi);
 end;
end;
 
procedure vyvod(a:mtr;s:string);
var i,j:byte;
begin
writeln(s);
for i:=1 to n do
 begin
  for j:=1 to n do
  write(a[i,j].z:4:1,'*e^i*(',a[i,j].fi:4:1,') ');
  writeln;
 end;
end;
 
var a,b,c:mtr;
begin
clrscr;
write('Введите размер матриц n=');
readln(n);
vvod(a,'A');
vvod(b,'B');
clrscr;
vyvod(a,'Исходная матрица А');
vyvod(b,'Исходная матрица B');
mult_mtr(a,b,c);
vyvod(c,'Матрица С=А*B');
readln
end.

  for i:=0 to m-1 do for j:=0 to n-1 do {with a[i,j] do}
    Read(a[i,j].re,a[i,j].im);