function xyz_new = testt(A, xyz)

% function that outputs the predicted xyz values in response to a
% bunch of input RGB test values. These RGB values are gamma
% corrected first, and then multiplied by the transformation
% matrix A.

rgb = [];
xyz = xyz(:,13:24);
for i = 4:6
  for j = 1:4
    [im,m] = imread(['scanner_data/' num2str(i) num2str(j) '.bmp']);
    [r,g,b] = getRGB(im);
    rgb = [rgb [r;g;b]];
  end
end
rgb = gammaCorrect(rgb);
if (size(A)==[3 6])
 rgb =[rgb(1,:);rgb(2,:);rgb(3,:);rgb(1,:).*rgb(2,:);rgb(2,:).*rgb(3,:);rgb(1,:).*rgb(3,:)];
end
rgb
xyz_new = A*rgb;



