Introduction - If you have any usage issues, please Google them yourself
Scalar function gradient procedures using Richardson extrapolation method, in the function richarson (fun, x0), h is the initial step, j is the number of iterations, precision is 2* j+2 order.