%sigma=var(cluster_center');
% sigma=zeros(10,20);
% for i=1:10
%     ind=find(cluster_idx==i);
%     V=[];
%     for j=1:size(ind)
%         V=[V, pca(:,j)];
%     end
%     sigma(i,:)=var(V');
% end

for t=4:-1:1
disp('toc');
end