Uses CRT;
Type
 Mas1=array[1..9,1..5] of real;
 Mas2=array[1..5] of real;
 Mas3=array[1..5,1..5] of real;

var
 a:Mas1;
 Mu:Mas2;
 Dz:Mas3;
 p,q:Integer;

procedure fill_massive(var FM:Mas1);
var i,j:Integer;
begin
 for i:=1 to 9 do
  for j:=1 to 5 do
   begin
    write('a[',i,',',j,']= ');
    readln(FM[i,j]);
   end;
end;

procedure average_meaning(AM1:Mas1; var AM2:Mas2);
var i,j:Integer;
begin
 for j:=1 to 5 do AM2[j]:=0;
 for j:=1 to 5 do
  for i:=1 to 9 do
   AM2[j]:=AM2[j]+AM1[i,j]/9;
end;

procedure Correlation(C1:Mas1; C2:Mas2; var C3:Mas3);
var i,j,k:Integer;
begin
 for i:=1 to 5 do
  for j:=1 to 5 do
   C3[i,j]:=0;
 for j:=1 to 5 do
  for k:=1 to 5 do
   for i:=1 to 9 do
    C3[j,k]:=C3[j,k]+(C1[i,j]-C2[j])*(C1[i,k]-C2[k])/9;
end;

Begin
 ClrScr;
 fill_massive(a);
 average_meaning(a,Mu);
 Correlation(a,Mu,Dz);
 ClrScr;
 Writeln('Initial matrix:');
 For p:=1 to 9 do
  begin
   For q:=1 to 5 do write(a[p,q]:8:3);
   Writeln;
  end;
 Writeln;
 Writeln('Correlation matrix:');
 For p:=1 to 5 do
  begin
   For q:=1 to 5 do write(Dz[p,q]:10:3);
   Writeln;
  end;
 ReadKey
End.