function [images,imgDiff,diffVect,diffRecon,diffReconVect] = modelFaces(data,U,xbar,sizeim)

weights = projectImage(data(:,1),U,xbar);
recon = reconstruct(weights,U,xbar,sizeim);

images = [];
imgDiff = [];
diffRecon = [];
for i=1:size(data,2)
weights = projectImage(data(:,i),U,xbar);
image = reconstruct(weights,U,xbar,sizeim);
size(image)
size(data(:,1))
images = [images image]; 
imgDiff = [imgDiff (data(:,i) - image)];
diffRecon = [diffRecon (recon - image)];

diff = (data(:,i) - image).^2;
diffR = (recon - image).^2;

diffVect(i) = sum(diff);
diffReconVect(i) = sum(diffR);
end