function Lmf = getLmf(m,n)
mn = m^n;
Lmf = zeros(mn,m*n);
MN = repmat([m],[1 n]);
rind = 0;
for xind=1:mn
  rind = rind+1;
  x = ind2coord(xind,MN);
  for t=1:n
    k = x(t);
    cind = coord2ind([k t],[m n]);
    Lmf(rind,cind) = 1;
  end
end
return