load irises1.dat

features=irises1;
%features=d;

[m, s, d] = normalize(features, [], []);
[datas, proj] = pca(d, [], 3);


plot3(datas(1:50,1),datas(1:50,2),datas(1:50,3),'*'); hold
plot3(datas(51:100,1),datas(51:100,2),datas(51:100,3),'+r')
plot3(datas(101:150,1),datas(101:150,2),datas(101:150,3),'.g')
rotate3d