compute_grad function

Compute the Numerical Gradient of the Log-Likelihood