function o = kmeans_demo(sigma)
k = 3;
ps = [];
for i = 1:k
    p = rand(2,1);
    ps = [ps, p*ones(1,20) + sigma*randn(2,20)];
end
kmeans(ps,k);