Произведение двух матриц с комплексными коэфициентами представленные в алгебраической форме записи

const max=10;
type complex=record//алгебраическая форма комплексного числа
             re,im:real;
             end;
     mtr=array[1..max,1..max] of complex;
function mult(a,b:complex):complex;//умножение чисел в алгебраической форме
begin
result.re:=a.re*b.re-a.im*b.im;
result.im:=a.re*b.im+b.re*a.im;
end;
function add(a,b:complex):complex;//сложение чисел в алгебраической форме
begin
result.re:=a.re+b.re;
result.im:=a.im+b.im;
end;
procedure mult_mtr(n:integer;a,b:mtr;var c:mtr);
var i,j,k:byte;
begin
for i:=1 to n do
for j:=1 to n do
 begin
  c[i,j].re:=0;
  c[i,j].im:=0;
  for k:=1 to n do
  c[i,j]:=add(c[i,j],mult(a[i,k],b[k,j]));
 end;
end;
procedure vvod(n:integer;var a:mtr);
var i,j:byte;
begin
for i:=1 to n do
for j:=1 to n do
 begin
  a[i,j].re:=-2+4*random;
  a[i,j].im:=-2+4*random;
 end;
end;
procedure vyvod(n,k:integer;a:mtr;s:string);
var i,j:byte;
begin
writeln(s);
for i:=1 to n do
 begin
  for j:=1 to n do
   begin
    write(a[i,j].re:k:1);
    if a[i,j].im<0 then write(a[i,j].im:k:1,'i ')
    else write('+',a[i,j].im:k-1:1,'i ');
   end;
  writeln;
 end;
end;
var n:integer;
    a,b,c:mtr;
begin
repeat
write('Введите размер матриц от 2 до ',max,' n=');
readln(n);
until n in [2..max];
vvod(n,a);
vyvod(n,4,a,'Исходная матрица А');
vvod(n,b);
vyvod(n,4,b,'Исходная матрица B');
mult_mtr(n,a,b,c);
vyvod(n,5,c,'Матрица С=А*B');
end.

procedure vvod(n:integer;var a:mtr);
var i,j:byte;
begin
writeln('Введите значения действительной и мнимой частей числа');
for i:=1 to n do
for j:=1 to n do
 begin
  readln(a[i,j].re);
  readln(a[i,j].im);
 end;
end;