function MU = findmu(K,phi0,m,n)


maxs = 0;
dims = repmat([m],[1 n]);

X = list_all_ind(dims);

[Ainf,Binf] = lipab(m,n-1);

MU = {};
global MINDS BINDS
MINDS = [];
BINDS = [];

for mind = 0:2^m-1
  mset = dec2bin(mind,m);
  muf = find(mset=='1');
  
  muv = ismember(X(:,1),muf);
  muv = reshape(muv,size(phi0));
  
  phi = phi0 - muv;
  
  phi1 = getK1(K.*phi',m,n);
  
  %g = all(Ainf*(phi1)<=Binf);
  g = all(Ainf*(phi1)<=Binf+1e-8);
  %gg(end+1) = max(abs(Ainf*(phi1)));
  if g,
    %break
    MU{end+1} = muv;
    MINDS(end+1) = mind;
else
    BINDS(mind+1) = max(abs(Ainf*(phi1)));
  end
  
  
end

return