 For an equidistant grid one gets a Toeplitz matrix.The 2D case shows all the characteristics of the more general nD case.Each second partial derivative needs to be approximated similar to the 1D case 
which is usually given by the following stencil

Consistency of the above-mentioned approximation can be shown for highly regular functions, such as formula_62.