How to vectorize this five-point difference code to fast calculate derivatives of a matrix?

Question

I am trying to calculate the partial first derivatives with respect to each of the two dimensions in a 2D matrix, i.e dF/dx and dF/dy, using the five point method. I have ma