function [dx,nev] = bs_gradient (f, args, narg) error("bs_gradient: bs_gradient has been replaced by numgradient"); endfunction