code: def test_stieltjes_evalf():
      assert abs(stieltjes(0).evalf() - 0.577215664) < 1e-9
      assert abs(stieltjes(0, 0.5).evalf() - 1.963510026) < 1e-9
      assert abs(stieltjes(1, 2).evalf() + 0.072815845) < 1e-9
  
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 14 for calling a global function with example usage: test_stieltjes_evalf() without return types