Skalarfeld

h:(x1^2)*(x2^3)+1;

Gradient

jacobian([h],[x1,x2]);

Hessematrix

hessian(h,[x1,x2]);