Maxima Function
hessian (f,vars)
Return the hessian matrix of f with respect to the variables in the list vars. The i,j entry of the hessian matrix is diff(f vars[i],1,vars[j],1).