Hessian matrix is an N x N matrix containing second partial derivatives of a function with respect to its input variables.
Hessian matrices are used extensively in second-order optimization methods such as Newton's method. They give information about the curvature around a point x0 on a function surface.