function [phi,fv] = maxphi(K)
global Ainf Binf n myopt
if isempty(Ainf)
  error('Ainf not initialized')
end
[xp,fvalp] = linprog( K,Ainf,Binf,[],[],0*K,0*K+n,0,myopt);
[xn,fvaln] = linprog(-K,Ainf,Binf,[],[],0*K,0*K+n,0,myopt);
if abs(fvalp)>abs(fvaln)
  fv = fvalp;
  x = xp;
else
  fv = fvaln;
  x = xn;
end
phi = x';
fv = abs(phi*K);
return



function phi = maxphi0(K);
global Ainl Binl Ain1 Bin1 GG iif m n mylinopt

[x,fval] = linprog(-K,Ainl,Binl,[],[],0*K,0*K+n,0,mylinopt);
phi = round(x');
return