N = 1000;
a = -2; b = +2;

v_unif = a + (b-a) * rand(N,2);
% N bidimensional random variables, with unif. distribution in [-1 +1]

figure;
plot(v_unif(:,1), v_unif(:,2),...
    'd','MarkerEdgeColor','b',...
        'MarkerFaceColor','b',...
        'MarkerSize',10);
grid on
hold on;

v_gauss = randn(N,2);
% gaussian random varaible, with expected value 0 and variance 1
plot(v_gauss(:,1), v_gauss(:,2),...
        'o','MarkerEdgeColor','r',...
            'MarkerFaceColor','r',...
            'MarkerSize',10);
axis square