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