How can one implement the numerical derivatives currently implemented in
Stata's ml maximizer in a setting outside outside of ml? I would like to
calculate numerical derivatives at given parameter values. The dydx function
does not appear to do this.
