First time here? Check out the FAQ!
answered 10 years ago
I guess
H = f.hessian().function(x,y) H(z[0,0],z[1,0])
would also work