function [x]=explore(h,n,f,x0)x = x0; f0 = eval(f);for i=1:n  x(i) = x(i) + h(i);  ff = eval(f);  if ff < f0,    f0 = ff;   else    x(i) = x0(i) - h(i);   ff = eval(f);    if ff < f0,   f0 = ff;   else,   x(i) = x0 (i);   end  endendreturn