% nacteni obrazku
for i = 1 : 13
   a = imread (sprintf ('subject%02d.tif', i));
   m (i,:) = double (a (:)')/255;
end;

% zobrazeni nactenych obrazku
figure;
for i = 1 : 13
   subplot (4, 4, i);
   imshow (reshape (m (i, :), 220, 170));
end;
   




% serazeni vlastnich vektoru podle vlastnich cisel

[v d] = eig (c);              % vypocet vl. vektoru a vl. cisel pro matici c
d = diag (d);                 % d byla diagonalni matice (vl. cisla na diagonale)
                              % ted je d jednorozmerny vektor
[d i] = sortrows (d);         % vytvoreni indexu serazenych vl. cisel, vzestupne
i = i (end:-1:1);             % otoceni indexu - nejvetsi napred
d = d (end:-1:1);             % otoceni i vektoru vlastnich cisel 
v = v (:, i);                 % zprehazeni vl. vektoru podle poradi vl. cisel
