Ny[0, 0] = -(z[1, 0] - z[0, 0]) / csz