ay tuned\n'</span>);
K = 100;
[B, mu, Xp, Xrecon] = pcaPmtk(X, K);
<span class="comment">%tic; [B] = pcaPmtk(X,K); toc</span>
Xproj = X*B;

<span class="comment">% visualize basis functions (eigenfaces)</span>
figure;
subplot(2,2,1)
f = reshape(mean_face, [h w]);
imagesc(f);
axis <span class="string">off</span>
<span class="keyword">for</span> i=1:3
  subplot(2,2,i+1)
  f = reshape(B(:,i), [h w]);
  imagesc(f);
  axis <span class="string">off</span>
<span class="keyword">end</span>
colormap <span class="str