function simmethodb

  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
    L = N/c;
    momsR = (1/2) .* ones(1,8);
    
    X = randn(N,L);
    eigenvals = eig((1/L) .* X * X' * R);
    momsapprox = momentsfromeigenvals(eigenvals,8);
    
    momsexact = convolvemp(momsR,c);
    mmse8 = [mmse8 (momsexact-momsapprox) * transpose(momsexact-momsapprox)];
    
    newmoms = convolvemp(momsR(1:4),c);
    mmse4 = [mmse4 (momsexact(1:4)-momsapprox(1:4)) * transpose(momsexact(1:4)-momsapprox(1:4))];
  end
    
  plot(Nvalues,mmse8,'k+');
  xlabel('N')
  ylabel('MMSE')
  axis([8 512 0 50]);
  print -deps ../images/simmethodb8.eps;
  
  plot(Nvalues,mmse4,'k+');
  xlabel('N')
  ylabel('MMSE')
  axis([8 512 0 50]);
  print -deps ../images/simmethodb4.eps;