code: def tange_debyetemp(v, v0, gamma0, a, b, theta0):
      
      x = v / v0
      gamma = tange_grun(v, v0, gamma0, a, b)
      if isuncertainties([v, v0, gamma0, a, b, theta0]):
          theta = theta0 * np.power(x, (-1. * (1. - a) * gamma0)) * \
              unp.exp((gamma0 - gamma) / b)
      else:
          theta = theta0 * np.power(x, (-1. * (1. - a) * gamma0)) * \
              np.exp((gamma0 - gamma) / b)
      return theta
docs: calculate Debye temperature for the Tange equation
  
      :param v: unit-cell volume in A^3
      :param v0: unit-cell volume in A^3 at 1 bar
      :param gamma0: Gruneisen parameter at 1 bar
      :param a: volume-independent adjustable parameters
      :param b: volume-independent adjustable parameters
      :param theta0: Debye temperature at 1 bar in K
      :return: Debye temperature in K