maxL = 0;
maxA = zeros(m,m,n-1);
maxF = zeros(m^n);
for t=1:10000
  m=4
  n=3
  h=.1;
  [A,phi,fval] = gamlip(m,n,h);
  if fval > maxL
    maxL = fval;
    maxA = A;
    maxF = phi;
  end
  maxA
  maxL
end