xy_new=fminsearch (@(v) f([x' y'],d,v,k),[x(k) y(k)]); x(k)=xy_new(1); y(k)=xy_new(2);