code: def _debye(alpha):
      # EPSILON = np.finfo(np.float32).eps
      EPSILON = np.finfo(np.float64).eps * 100
  
      def integrand(t):
          return np.squeeze(t / (np.exp(t) - 1))
  
      _alpha = np.squeeze(alpha)
      debye_value = (
          integrate.quad(integrand, EPSILON, _alpha)[0] / _alpha
      )
      return debye_value
  
imports: import sys
  import numpy
  from scipy import stats, integrate, optimize
  from .None import transforms
  from .copulas import Copula
  from statsmodels.tools.rng_qrng import check_random_state
label_desc: Use image node_id 1 for calling a global function with example usage: _debye(alpha) and returns: debye_value, np