function x=grm(a,b,n)
  r=(sqrt(5)-1)/2;
  c=a+(1-r)*(b-a);
  d=a+r*(b-a);
  for i=1:n 
    if fv(c)>fv(d);
      a=c;
      b=b;
      c=d;
      d=a+r*(b-a);
    else 
      a=a;
      b=d;
      d=c;
      c=a+(1-r)*(b-a);
    endif
  endfor
  x=(a+b)/2
endfunction