function [y] = kmeans(X, k)
[n,d] = size(X);
% Pick k points at random as initial centers. Do this by randomly
% permuting the points and then taking the first k in the permutation
centers = X(randperm(n),:);
centers = centers(1:k,:);
for i = 1:100
% Fill in the k-means updates here
end