% PROVA ORDINI (DA METTERE DENTRO IL CODICE)
% W = Yc*Yc';
% wmin=1;
% eig(W)
% W = W + wmin*eye(size(W));
% keyboard
% cond(W)
% 
% W = real(chol(W));
% Wp = real(pinv(W'));
% clear Yc ;
% clear Zc ;


%fpe sottospazi     


IVC(n) = log(det(Lambda))+log(N)/N*((p+m)*n+p*n);

dS = diag(S); (valori singolari)


SVC(n)=dS(n+1)+log(N)/N*((p+m)*n+p*n);