function babyplot
global a b

m=2
n=2

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

x = rand(4,1);

global P U
C = ones(2,1)/2;
plot3(0,0,0);
for a=0:.01:1
  for b=0:.01:a
        
    A(1,1) = a; A(2,1) = 1-a;
    A(1,2) = b; A(2,2) = 1-b;
    P = markfillprob(A,C,2);
    U = prodmeas(P,2,2);
    
    lhs = P'*x;
    rhs = U'*x;
    rhs = rhs*(1+abs(a-b));
    
    F = rhs - lhs;
    hold on
    plot3(a,b,F,'.b');
    drawnow
  end
end

return