ss="comment">%for i=1:25</span>
<span class="comment">%  subplot(5,5,i)</span>
<span class="comment">%  f = reshape(X(i,:), [h w]);</span>
<span class="comment">%  imagesc(f);</span>
<span class="comment">%  axis off</span>
<span class="comment">%end</span>
<span class="comment">%colormap gray</span>

Selected_Face = 125;<span class="comment">%312;%255;</span>

[N,D]=size(X);
mean_face = mean(X);
<span class="comment">%X = X - repmat(mean_face,N,1);</span>
fprintf(<span class="string">'Performing PCA.... st