Schuster_hessian function

Computes the hessian matrix of second-order partial derivatives of log(L).