function [x1,x2,xind1,xind2] = randrho1(m,n)
% rho(x1,x2) = 1
dims = repmat([m],[1 n]);
xind1 = ceil(randinseg(0,m^n));
x1 = ind2coord(xind1,dims);
t = ceil(randinseg(0,n));
mm = setdiff(1:m,x1(t));
rp = randperm(m-1);
k = mm(rp(1));
x2 = x1;
x2(t) = k;
xind2 = coord2ind(x2,dims);
return