Python 2nd order finite different approximation for multivariate functions

Question

How do you calculate delf/delx and delf/dely numerically using the 2nd order finite difference approximation for a multivariate function in python?