def METHOD_NAME(self):
    gamma = 2.5
    rc = 0.01
    rho0 = 1.0
    R = 5.0
    args = (rho0, rc, gamma)
    mass_numerical = quad(self._mass_integrand3d, 0, R, args=args)[0]
    mass_analytic = self.profile.mass_3d(R, rho0, rc, gamma)
    npt.assert_almost_equal(mass_analytic, mass_numerical)
    gamma = 2.0
    args = (rho0, rc, gamma)
    mass_numerical = quad(self._mass_integrand3d, 0, R, args=args)[0]
    mass_analytic = self.profile.mass_3d(R, rho0, rc, gamma)
    npt.assert_almost_equal(mass_analytic, mass_numerical)
    gamma = 3.0
    args = (rho0, rc, gamma)
    mass_numerical = quad(self._mass_integrand3d, 0, R, args=args)[0]
    mass_analytic = self.profile.mass_3d(R, rho0, rc, gamma)
    npt.assert_almost_equal(mass_analytic, mass_numerical)
    gamma = 1.4
    args = (rho0, rc, gamma)
    mass_numerical = quad(self._mass_integrand3d, 0, R, args=args)[0]
    mass_analytic = self.profile.mass_3d(R, rho0, rc, gamma)
    npt.assert_almost_equal(mass_analytic, mass_numerical)
    sigma0 = rho0 * rc
    mass_analytic_from_sigm0 = self.profile.mass_3d_lens(R, sigma0, rc, gamma)
    npt.assert_almost_equal(mass_analytic_from_sigm0, mass_numerical)