Make run again by using a FD approximation to compute the gradient and hessian

I didn't test the changes because I don't have a proper parameter set. Please run and test before merging.