function [] = vecgaus(rho)
n=1000;
u1=rand(1,n);
u2=rand(1,n);
g1=sqrt(-2*log(u1)).*cos(2*%pi*u2);
g2=sqrt(-2*log(u1)).*sin(2*%pi*u2);
A=[1,0;rho,sqrt(1-rho^2)];
z=A*[g1;g2];
x=z(1,:);
y=z(2,:);
plot2d(x,y,-1)
endfunction


