% Scatter plot

n = 1000;
x = pi * randn(n,1);
y = pi * randn(n,1);

c = cos(x) .* cos(y);

scatter(x,y,4,c);
colorbar
