function [Ainf,Binf] = lipabd(m,n,dd)
mn = m^n;
L = mn;
Ainf = zeros(L,mn);
Binf = zeros(L,1);
MN = repmat([m],[1 n]);
rind = 0;
for ind1 = 1:mn
  for ind2 = ind1+1:mn
    x = ind2coord(ind1,MN);
    y = ind2coord(ind2,MN);
    d = ddist(dd,x,y);
    for s = [-1,1]
      rind = rind+1;
      Ainf(rind,ind1) =  s;
      Ainf(rind,ind2) = -s;
      Binf(rind,1) = d;
    end
  end
end
return