X=randn(100,5); Xn=X-mean(X); isequal(((Xn'*Xn)./100)./cov(X,1),ones(5))