n1 = n-1;
phi = zeros(m^n,1);
phi1 = zeros(m^n1,1);

LBK = -ones(size(phi));
UBK =  ones(size(phi));

LBg =   zeros(size(phi));
%UBg = n+zeros(size(phi));
%LBg = -n*ones(size(phi));
UBg =  n*ones(size(phi));

LBg1 =    zeros(size(phi1));
UBg1 = n1+zeros(size(phi1));

[A,B] = lipab(m,n);
[A1,B1] = lipab(m,n1);