function plotvandermonde3(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)
      V1 = zeros(N,L);
      V2 = zeros(N,L);
      for (k=1:L) 
        theta = 2*angl*(rand(1)-0.5);
        V1(:,k) = ((exp(2*pi*j*sin(theta)*d/lambda)).^(0:(N-1)))';
        theta = 2*angl*(rand(1)-0.5);
        V2(:,k) = ((exp(2*pi*j*sin(theta)*d/lambda)).^(0:(N-1)))';
      end
      V1 = (1/sqrt(N))*V1;
      V2 = (1/sqrt(N))*V2;
      capacitiesV = [ capacitiesV (1/N)*log2(det( eye(L,L) + rho * V1' * V2 * V2' * V1 )) ];
      
      rhos = [rhos rho];
    end
  end
  
  plot(rhos,capacitiesV,'k.');
  xlabel('\rho');
  ylabel('Capacity');
  axis([0 10 0 3]);
  print -deps ../images/vandermondesamples3.eps
  
  
  
  
 