N=1000;
x=matrix(runif(N),,1);
y=matrix(runif(N),,1)*x;
plot(x,y,pch=18,col="blue")