curmin = 10
global m n
m=3
save tmpdat A curmin

global HH PP Ainf Binf delta Lmf Lpd
global myopt
myopt = optimset('Diagnostics','off');
myopt = optimset('Display','off');
warning off

rho = .5001

global phi01 domask
domask = 0;
phi01 = 0;

for n=3
  [Ainf,Binf] = lipab(m,n);
  Lmf = getLmf(m,n);
  curmin = 10
  for t=1:50
    %[A,v]=minmaxvar(m,n,.1)
    %[A,v] = minmaxvar_hcon(m,n,n-.1);
    [A,v] = minmaxvar_hcon(m,n,rho);
    %[A,v]=minmaxvar_incr(m,n,.1,CTRX{end})
    
    if v<curmin
      curmin = v;
      Amin = A;
      fprintf('n=%d iter = %d  curmin = %1.9f \n',n,t,curmin);
      save tmpdat Amin curmin CTRX
    end
    if curmin < 1e-7
      break
    end
  end
  CTRX{n}=Amin;
end