function [c vox_comuni] = costo_stat(cams,cam1,cam2)
%
% INPUT  : le due telecamere che confronto cam1 e cam2
%          cam1 e cam2 rappresentano le etichette delle telecamerea da
%          confrontare
%
% OUTPUT : matrice costo c tra le due telecamere
%          sarebbe utile creare tabella che mostri tutte le c


%   Detailed explanation goes here
vox1=cams(cam1).voxels;
vox2=cams(cam2).voxels;
fine1=size(vox1,2);
fine2=size(vox2,2);
c=0;
fine=0;
i=1;
j=1;
vox_comuni=[];
while fine==0
    a=vox1(i);
    b=vox2(j);
    if a==b
        c=c+1;
        vox_comuni=[vox_comuni a];
        i=i+1;
        j=j+1;
    end
    if a>b
        vox_comuni=[vox_comuni b];
        j=j+1;
    end
    if a<b
        vox_comuni=[vox_comuni a];
        i=i+1;
    end
    if i>fine1 || j>fine2
        fine=1;
    end
end
cams(cam1).versore;
cams(cam2).versore;
phi=acos(cams(cam1).versore'*cams(cam2).versore)/(pi/2);
c=c*abs(sin(phi));

end
