手続き hessian(expr, vars) は expr の vars に関するヘッシアン行列を計算します。

計算される行列は、n x n です。ただし、n は vars の長さです。計算される行列の (i,j) 成分は diff(expr, vars[i], vars[j]) です。

コマンド with(linalg,hessian) により、このコマンドの省略形を使うことができます。

with(linalg):

Warning, the protected names norm and trace have been
redefined and unprotected

hessian(x*y*z, [x,y,z]);

hessian(x^2*y + 3*x*y^2, [x,y]);