[xlev,vlev,nlev] = feval(optim_func, "ff", "dff", xinit) ;

Package: optim