function suphisketch(m0,n0) global m n randomize m = m0; n = n0; global pstr pstr = mfilename global lstr rstr global Ainf Binf mylinopt mylinopt = optimset('Diagnostics','off'); mylinopt = optimset('Display','off'); [Ainf,Binf] = lipab(m,n); global mind mind = 1e6; while 1 x = randn(m^n,1); y = fmincon(@myobj,x,[],[],[],[],0*x-10,0*x+10); fprintf('%s: PARAMS: n=%d m=%d; \n',pstr,n,m); fprintf('DIFF _________________________ mind = %1.9f \n',mind); end return function y = myobj(K) global m n global Ainf Binf mylinopt global pstr [x,fval] = linprog(-K,Ainf,Binf,[],[],0*K+1,0*K+n,0,mylinopt); lhs = -fval; [x,fval] = linprog(-K,Ainf,Binf,[],[],0*K,0*K+n,0,mylinopt); rhs = -fval + sum(K); y = rhs - lhs; global mind if (y < mind) global mindK mindK = K; mind = y; fprintf('%s: PARAMS: n=%d m=%d; \n',pstr,n,m); fprintf('DIFF _________________________ mind = %1.9f \n',mind); if y < -.01 save BADESUPH m n K y end end return