f</span>
    title(<span class="string">'Reconstructed Image with 10 bases'</span>);drawnow;
  <span class="keyword">elseif</span> i==100
      subplot(2,2,3)
      imagesc(reshape(X_Reconst',h,w)); axis <span class="string">off</span>
      title(<span class="string">'Reconstructed Image with 100 bases'</span>);drawnow;
  <span class="keyword">end</span>
  recon_err = [recon_err;sqrt(mean((X_Reconst - (X(Selected_Face,:) + mean_face) ).^2,2))];
  colormap <span class="string">gray</span>
  subplot(2,2,4)
 