Computation of Hessian matrices