function [mu,phi1,K1,maxs] = boundsumnl(K,m,n)


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

K1 = getK1(K,m,n);

LB = 1+zeros(m^(n-1),1);
UB = n+zeros(m^(n-1),1);
[x,fval] = linprog(-K1,[],[],[],[],LB,UB);

phi1 = x'-min(x);
 

K0 = getK0(K,m,n);
LB =   zeros(m,1);
UB = 1+zeros(m,1);
[x,fval] = linprog(-K0,[],[],[],[],LB,UB);


mu = x'-min(x);

maxs = phi1*K1 + mu*K0;

return