numerical_deriv function

Compute numerical derivatives