function y = myF(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 = sum(pl(K)) + phi1*K1 - phi*K; return