[x,v,niter] = feval (optim_func, "testfunc","dtestf", xinit);

Package: optim