Plain run, just to make sure ######################################
 Minimum wrt 'x' is y0
 [xlev,vlev,nlev] = feval (optim_func, "ff", "dff", {x0,y0,1});
 ctl.df = "dff";

Package: optim