function phi1 = maxlip(K,phi,m,n) dims = repmat([m],[1 n]); pind = find(K>0); posx = ind2coord(pind,dims); K1 = getK1(K,m,n); phi1 = zeros(m^(n-1),1); for x1nd = 1:m^(n-1) x1 = ind2coord(x1nd,dims(2:end)); ff = zeros(m,1); for k=1:m x = [k x1]; xind = coord2ind(x,dims); if K(xind)>0 ff(k) = max(0,phi(xind)-1); else ff(k) = phi(xind); end end if K1(x1nd) > 0 phi1(x1nd) = max(ff); else phi1(x1nd) = min(ff); end end return