code: def _compute_c(x, z, eps=0.00001):
      """
      this function returns the mle of the shape parameter if a 1D gamma
      density
      """
      eps = 1.0e-7
      y = np.dot(z, np.log(x)) / np.sum(z) - np.log(
          np.dot(z, x) / np.sum(z)
      )
      if y > -eps:
          c = 10
      else:
          c = _psi_solve(y, eps=0.00001)
      return c
  
imports: import numpy
  import scipy.special
  import scipy.stats
label_desc: Use image node_id 3 for calling a global function with example usage: _compute_c(x, z, eps) and returns: c