function r=rho(x,y)
r = sum(x~=y);
return