code: def dg_series(z, n):
      """Symbolic expansion of digamma(z) in z=0 to order n.
  
      See https://dlmf.nist.gov/5.7.E4 and with https://dlmf.nist.gov/5.5.E2
      """
      k = symbols("k")
      return (
          -1 / z
          - EulerGamma
          + sympy.summation(
              (-1) ** k * zeta(k) * z ** (k - 1), (k, 2, n + 1)
          )
      )
  
file_path: catboost/contrib/python/scipy/py3/scipy/special/_precompute/wright_bessel.py
imports: from argparse import ArgumentParser, RawTextHelpFormatter
  import numpy
  from scipy.integrate import quad
  from scipy.optimize import minimize_scalar, curve_fit
  from time import time
label_desc: Use image node_id 2 for calling a global function with example usage: dg_series(z, n) and returns: unknown