hookjeeves.m

来自「various matlab programs to slove various」· M 代码 · 共 12 行

function [x,minf,nit]=hookejeeves(n,h,f,x0,toll)
x = x0; minf = eval(f); nit = 0;
while h > toll
  [y] = explore(h,n,f,x);
  if y == x,  h = h/2; else
    x = 2*y-x;     [z] = explore(h,n,f,x);
    if z == x,   x = y;  else,   x = z;  end
  end
  nit = nit +1;
end
minf = eval(f);
return