function y = fnorm(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 = phi*K;
%y = fval;
return