function [indices, q] = findClusters( newdataset , k )
%
% University of Texas at El Paso
% PCAT
%
% do k-means clustering

indices = kmeans(newdataset, k, 'display', 'off', 'replicates', 10);

% evaluate the clustering quality

silh = silhouette(newdataset, indices);
q = mean(silh);