function [A]=correl(C)
%--------------------------------------------------------------------------
%Esta función calcula la matriz de correlaciones de los vectores
%pertenecientes a la matriz C. Si esta es de NxM, entonces A es una matriz
%de NxN con unos en la diagonal principal
%--------------------------------------------------------------------------
G=size(C);
n=G(1,1);
t=G(1,2);
M=mean(C,1);
S=std(C,1);
for i=1:t
    for j=1:t
        vecaux1=C(:,i)-M(i);
        vecaux2=C(:,j)-M(j);
        vec1=(vecaux1.'*vecaux2)/n;
        A(i,j)=vec1/(S(i)*S(j));
    end
end