def METHOD_NAME(self):
    R = 1
    Rs = 3
    alpha_Rs = 1
    m_3d_nfw = self.nfw.mass_3d_lens(R, Rs, alpha_Rs)
    m_3d_cse = self.nfw_cse.mass_3d_lens(R, Rs, alpha_Rs)
    npt.assert_almost_equal(m_3d_nfw, m_3d_cse, decimal=8)