function [U,D]=eigenface_space_extractor(sample_vectors,space_dimension)

covariance_matrix = sample_vectors' * sample_vectors;
[V,D] = eig(covariance_matrix);

complete_space_dimension=size(covariance_matrix,1);

for i=complete_space_dimension-space_dimension+1:complete_space_dimension
    U(:,i-(complete_space_dimension-space_dimension)) = sample_vectors*V(:,i);
end