code: def test_issue_8404():
      i = Symbol("i", integer=True)
      assert (
          Abs(
              Sum(1 / (3 * i + 1) ** 2, (i, 0, S.Infinity)).doit().n(4)
              - 1.122
          )
          < 0.001
      )
  
imports: from sympy.concrete.summations import Sum
  from sympy.core.function import expand_func
  from sympy.core.numbers import Float, I, Rational, nan, oo, pi, zoo
  from sympy.core.singleton import S
  from sympy.core.symbol import Symbol
  from sympy.functions.elementary.complexes import Abs, polar_lift
  from sympy.functions.elementary.exponential import exp, exp_polar, log
  from sympy.functions.elementary.miscellaneous import sqrt
  from sympy.functions.special.zeta_functions import dirichlet_eta, lerchphi, polylog, riemann_xi, stieltjes, zeta
  from sympy.series.order import O
  from sympy.core.function import ArgumentIndexError
  from sympy.functions.combinatorial.numbers import bernoulli, factorial, genocchi, harmonic
  from sympy.testing.pytest import raises
  from sympy.core.random import test_derivative_numerically, random_complex_number, verify_numerically
label_desc: Use image node_id 10 for calling a global function with example usage: test_issue_8404() without return types