Nx[1:n-1, m-1] = (2*(z[1:n-1, m-2] - z[1:n-1, m-1]) + z[0:n-2, m-2]
                          - z[0:n-2, m-1] + z[2:n, m-2] - z[2:n, m-1]) / csz