def METHOD_NAME(self):
    Rs = 0.2
    alpha_Rs = 0.1
    r_trunc = 1000000000000 * Rs
    x = np.linspace(0.1 * Rs, 5 * Rs, 100)
    y = np.linspace(0.2, 1, 100)
    xxt, xyt, yxt, yyt = self.tnfw.hessian(x, y, Rs, alpha_Rs, r_trunc)
    xx, xy, yx, yy = self.nfw.hessian(x, y, Rs, alpha_Rs)
    np.testing.assert_almost_equal(xy, xyt, 4)
    np.testing.assert_almost_equal(yy, yyt, 4)
    np.testing.assert_almost_equal(xy, xyt, 4)
    np.testing.assert_almost_equal(yxt, xyt, 8)
    Rs = 0.2
    r_trunc = 5
    xxt, xyt, yxt, yyt = self.tnfw.hessian(Rs, 0, Rs, alpha_Rs, r_trunc)
    xxt_delta, xyt_delta, yxt_delta, yyt_delta = self.tnfw.hessian(
        Rs + 0.000001, 0, Rs, alpha_Rs, r_trunc
    )
    npt.assert_almost_equal(xxt, xxt_delta, decimal=6)