To compute the gradient of a function, it can be made with matematica using Grad(f(x,y),{x,y})

Suprisingly, it’s not possible to compute the Hessian matrix (something as HessianMatrix(f(x,y),{x,y}), doesn’t exist (I check in Wolfram website here, but nothing relevant). Is it possible that mathematica can’t compute the Hessian of a function ? And if no (what I guess), how can I do ?