function plotvandermonde2(N)
  
  L=N;
  
  angl = pi/4;
  d = 1;
  lambda = 10;
  
  rhos = [];
  capacitiesV = [];
  
  
  for (rho=0.1:0.1:10) 
    rho
    for (numsamples=1:1:10)
      V = zeros(N,L);
      for (k=1:L) 
        theta = 2*angl*(rand(1)-0.5);
        V(:,k) = ((exp(2*pi*j*sin(theta)*d/lambda)).^(0:(N-1)))';
      end
      
      V = (1/sqrt(N))*V;
      capacitiesV = [ capacitiesV (1/N)*log2(det( eye(N,N) + rho * V * V' )) ];
      
      rhos = [rhos rho];
    end
  end
  
  plot(rhos,capacitiesV,'k.');
  xlabel('\rho');
  ylabel('Capacity');
  axis([0 10 0 3]);
  print -deps ../images/vandermondesamples2.eps
  
  
  
  
 