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
