function cp    = mm2vox(mm,mat);
cp = inv(mat)*[mm; ones(1,size(mm,2))];
cp = cp(1:3,:);