function y = myF2(x) global Ainl Binl Ain1 Bin1 GG iif m n mylinopt Aing Bing K = x(:); %K' %n1 = n-1; %K1 = getK1(K,m,n); [x,fval] = linprog(-K,Ainl,Binl,[],[],0*K,0*K+n,0,mylinopt); phi = x'; %[x,fval] = linprog(-K1,Ain1,Bin1,[],[],0*K1,0*K1+n1,0,mylinopt); %phi1 = x'; y = phi*K - sum(pl(K)); return