function [A,gam,fval] = gamrmaltlp(m,n,h,d)
% exploring |Ag|
global nstate
nstate=m;

T = 1;

%A = pnormdim(ones(m)+rand*h,1);
A = pnormdim(rand(m),1);
gam = rand(m)/1000;
while 1
  [A,fval2] = gamrmlp2(gam,h);
  [gam,fval1] = gamrmlp(A,d,n);
  if abs(fval1-fval2)<1e-8
    break
  end
end

fval = fval1;

return