close all
set(0,'DefaultAxesFontSize', 24)
load ../Data/mit_cmu_faces.mat

% run PCA
[PC, score, eigenvalues, tsquared, explained] = pca(DATA);