function simgirko

  cd Z:/Education/cvs/uio/postdoc/luminy/matlab
  set(gca,'FontSize',24);
  
  c =0.5;
  
  mmse4 = [];
  mmse8 = [];
  Nvalues = 8:4:512;
  for N=Nvalues
    sprintf('N=%i\n',N)
    
    R = zeros(N,N);
    for (k=1:N/2) 
      R(k,k) = 1;
    end
    momsR = (1/2) .* ones(1,8);
    L = N/c;
    
    X = randn(N,L);
    eigenvals = eig((1/L) .* X * X' * R);
    moms = momentsfromeigenvals(eigenvals,8);
    
    newmoms = deconvolvemp(moms,c);
    mmse8 = [mmse8 (newmoms-momsR) * transpose(newmoms-momsR)];
    
    
    newmoms = deconvolvemp(moms(1:4),c);
    mmse4 = [mmse4 (newmoms-momsR(1:4)) * transpose(newmoms-momsR(1:4))];
  end
    
  plot(Nvalues,mmse8,'k+');
  xlabel('N')
  ylabel('MMSE')
  ax = axis;
  bx = ax(4);
  axis([8 512 0 1]);
  print -deps ../images/simgirko8.eps;
  
  plot(Nvalues,mmse4,'k+');
  xlabel('N')
  ylabel('MMSE')
  axis([8 512 0 1]);
  print -deps ../images/simgirko4.eps;