hessian function

Compute Hessian matrix at the MLE