
k=50;
x=grand(1,k,'nor',0,1);
y=gsort(x,'c','i');
z=[];
s=0;
for i=1:k
    z=[z,[y(i),y(i);s,s+1/k]];
    s=s+1/k;
end
clf();
plot(z(1,:),z(2,:),'bl');
