Вычислить значение функции

const
  n=3;
type
  d = array[1..n,1..n] of real;
const
  a: d = ((  1.2, -1.3,  2.1),
          (  0.2, -0.4,  0.9),
          (  2.1,  0.5,  4.2));
  b: d = ((  0.3, -2.1,  4.2),
          ( -1.4,  0.5,  3.4),
          (  0.7, -1.5, -3.1));
 
  
procedure output(a: d);
var i,j: integer;
begin
  writeln('Полученная матрица:');
  for i:=1 to n do
  begin
    for j:=1 to n do
      write(a[i,j]:7:1);
    writeln;
  end;
  writeln;
end;
 
procedure multiply(a,b: d; var c: d);
var i,j,p: integer;
begin
  for i:=1 to n do
  begin
    for j:=1 to n do
    begin
      c[i,j]:=0;
      for p:=1 to n do
        c[i,j]:=c[i,j]+a[i,p]*b[p,j];
    end;
  end;
end;
 
function minim(c: d; var im,jm: integer): boolean;
var i,j: integer;
    f: boolean;
begin
  f:=true;
  for i:=1 to n do
    for j:=1 to n do
      if abs(c[i,j])<2 then
      begin
        f:=false;
        im:=i;
        jm:=j;
        break;
      end;
 
  for i:=1 to n do
    for j:=1 to n do
      if (abs(c[i,j])<2) and (c[i,j]<c[im,jm]) then
      begin
        im:=i;
        jm:=j;
      end;
  minim:=f;
end;
 
function avarage(c: d): real;
var i,j: integer;
    s,r: real;
begin
  s:=1;
  for i:=1 to n do
    for j:=1 to n do
      s:=s*c[i,j];
  r:=exp(1/n/n*ln(abs(s)));
  if (odd(n*n)) and (s<0) then r:=-r;
  avarage:=r;
end;
 
var
  c: d;
  k,im,jm: integer;
  y,p,r: real;
begin
  multiply(a,b,c);
  output(c);
  if minim(c,im,jm) then writeln('Нет элементов, удовлетворяющих условию') else
  begin
    p:=c[im,jm];
    r:=avarage(c);
    writeln('p = ',p:0:5,'; r = ',r:0:5);
    y:=1/5*ln(1+p*p+r*r)/ln(5);
    writeln('y = ',y:0:5);
  end;
  readln;
end.