hessian function

Calculate hessian for a fitted model