answered 11 years ago
Does it answers your question ?
sage: x, y = var('x y') sage: f = log(x+y) sage: f.hessian() [-1/(x + y)^2 -1/(x + y)^2] [-1/(x + y)^2 -1/(x + y)^2]