function y = phinorm(K)
global Ainf Binf m n mylinopt

[x,fvalp] = linprog( K,Ainf,Binf,[],[],0*K,0*K+n,0,mylinopt);
[x,fvaln] = linprog(-K,Ainf,Binf,[],[],0*K,0*K+n,0,mylinopt);
y = max(abs(fvalp),abs(fvaln));
return



function y = phinorm0(K)
global Ainl Binl Ain1 Bin1 GG iif m n mylinopt Aing Bing

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

% should really consider max over +/- K

y = abs(phi*K);
%y = fval;
return