% function [rgb3] = Model3Rgb(rgb);
% Returns matrix with r,g,b, and rg,gb,br,rr,gg,bb rows
% input: rgb, 3xN matrix
% output: rgb3, 9xN matrix
function [rgb3] = Model3Rgb(rgb);

rgb3 = zeros(9,size(rgb,2));
rgb3(1:3,:)=rgb(:,:);
rgb3(4,:)=rgb(1,:).*rgb(2,:);
rgb3(5,:)=rgb(2,:).*rgb(3,:);
rgb3(6,:)=rgb(3,:).*rgb(1,:);
rgb3(7,:)=rgb(1,:).*rgb(1,:);
rgb3(8,:)=rgb(2,:).*rgb(2,:);
rgb3(9,:)=rgb(3,:).*rgb(3,:);