function [] = gaussian_vector(rho)
  if abs(rho) >1
	disp('la corrélation doit être comprise entre -1 et 1!')
	disp('aborting...')
	return
  end
  xbasc();
  n=1000;
  u1=rand(1,n);
  u2=rand(1,n);
  g1= ZZ // à compléter 
  g2= ZZ // à compléter
  A=[1,0;rho,sqrt(1-rho^2)];
  z=A*[g1;g2];
  x=z(1,:);
  y=z(2,:);
  plot2d(x,y,-1)
endfunction