code: def speziale_debyetemp(v, v0, gamma0, q0, q1, theta0):
      
      if isuncertainties([v, v0, gamma0, q0, q1, theta0]):
          f_vu = np.vectorize(uct.wrap(integrate_gamma),
                              excluded=[1, 2, 3, 4, 5, 6])
          integ = f_vu(v, v0, gamma0, q0, q1, theta0)
          theta = unp.exp(unp.log(theta0) - integ)
      else:
          f_v = np.vectorize(integrate_gamma, excluded=[1, 2, 3, 4, 5, 6])
          integ = f_v(v, v0, gamma0, q0, q1, theta0)
          theta = np.exp(np.log(theta0) - integ)
      return theta
docs: calculate Debye temperature for the Speziale 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 q0: logarithmic derivative of Gruneisen parameter
      :param q1: logarithmic derivative of Gruneisen parameter
      :param theta0: Debye temperature at 1 bar in K
      :return: Debye temperature in K