function babysketch(m0,n0)
global m n
m=m0;n=n0;

C = ones(m,1)/m;

  global Ainf Binf mylinopt
  mylinopt = optimset('Diagnostics','off');
  mylinopt = optimset('Display','off');
  [Ainf,Binf] = lipab(m,n);

while 1
    
  P = markfillprob(A,C,n);
  U = prodmeas(P,m,n);
  
  x = rand(m^n,1);
  %F = randlip;
  
  lhs = F*P;
  %rhs = (1+abs(a-b))*F*U;
  tha = abs(a-b);
  Hn = sum(tha.^[0:n-1]);
  rhs = Hn*F*U;
  
  y = rhs - lhs;
  
  %maxr = 0;
  %x = fmincon(@radrat,x,[],[],[],[],0*x+saf,0*x+raf);
  %x = babyaux(x);
  %if maxr > 1.01
  %  keyboard
  %end
  %rr(end+1) = y;
  hold on
  %plot(b,maxr,'.b');
  plot(b,y,'.b');
  drawnow
end

return



function phi = randlip
global Ainf Binf m n mylinopt
global saf raf
K = genK(m,n);
[x,fval] = linprog(-K,Ainf,Binf,[],[],0*K+saf,0*K+raf,0,mylinopt);
phi = x';
return