function babysketch(a0,n0)
global a b
a = a0;
n = n0;
global m n

m=2
%n=2

global saf raf
saf = 1;
raf = n;
%raf = 10;

global P U
C = ones(2,1)/2;

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

%a = .2;
rr = [];
global maxr
  %F = randlip;
  %F = rand(1,m^n);
  F = [2 1 1 1]

for b=0:.01:1
  A(1,1) = a; A(2,1) = 1-a;
  A(1,2) = b; A(2,2) = 1-b;
  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