code: def _harmonic_number(b):
      b = np.asarray(b)
      return sp_special.psi(b + 1) - sp_special.psi(1)
  
imports: import numpy
  from scipy import special
  from scipy import stats
  import tensorflow.compat.v1
  import tensorflow.compat.v2
  from tensorflow_probability.python.distributions import kumaraswamy
  from tensorflow_probability.python.internal import test_util
label_desc: Use image node_id 3 for calling a global function with example usage: _harmonic_number(b) and returns: unknown