a = randint(3,3,[1 9])
[n,m] = size(a);
for i = 1:m
ai = a(:,i);
mi = mean(ai);
for j = 1:m
aj = a(:,j);
mj = mean(aj);
r(i,j) = sum((ai-mi).*(aj-mj))/(n-1);
end
end
r
a =
5 9 4
6 7 9
8 2 9
r =
2.3333 -5.5000 3.3333
-5.5000 13.0000 -7.5000
3.3333 -7.5000 8.3333