p = 3;
k = 3;

%X = toydata(p, 200, 1);
X = toydata(p,10,1);

figure;
hold all

y = kmeans(X, k, 1);
for i = 1:k
    plot(X(y==i,1), X(y==i,2), '+');
end
