n=10000
    
r=sqrt(rand(1,n));
theta=2*pi*rand(1,n);

x=r.*sin(theta);
y=r.*cos(theta);

scatter(x,y)