PyTorch: Higher order derivates of bivariate function

Question

Suppose I have a bivariate function f(x,y), and I have two sets of inputs X = [x1,x2,...,xm] and Y = [y1,y2,...,yn] and I would like to c