fp=fopen('prtcl.txt ','w');
XX=load('dan_vuz.txt');
size(XX);
X=XX(:,3:13);
R=corr(X);
[vect,lambda]=eig(X'*X);
Sobst=diag(lambda);
fprintf(fp,'Eigenvalues:\n %f \n',Sobst);
fprintf(fp,'\n');
SobMax=Sobst(end);
GlComp=vect(:,end);
Delt=100*SobMax/sum(Sobst);
fprintf(fp,'Delta= %d \n ',round(Delt));
Res=X*GlComp;
fprintf(fp,' Results \n ');
fprintf(fp,'%d  %f \n ',[XX(:,1),Res] ');
MinRes= min(Res);
MaxREs= max(Res);
MeanRes= mean(Res);
stdRes= std(Res);
fprintf(fp,'Min=%f\n', MinRes);
fprintf(fp,'Max=%f\n', MaxRes);
fprintf(fp,'Mean=%f\n', MeanRes);
fprintf(fp,'Std=%f\n', stdRes);
save res.mat Res -mat;
hist(Res,20);
xlabel('Results ');
ylabel('Number of Unis ');
saveas(gcf, 'Hist.jpg ', 'jpg ');
CorFin=corr(Res,XX(:,2));
fprintf(fp,'Correlation of Results and Money = %f \n',CorFin);
fclose(fp);