debye

mpsci.fun.debye(x, *, n, method='quad')

Compute the Debye function D_n(x).

\(D_n(x) = \frac{n}{x^n} \int_{0}^{x}\frac{t^n}{e^t - 1} dt\)

If method == 'quad', mp.quad is used to evaluate the integral form of the function, so it might be necessary to set mp.dps to a value much larger than the desired output precision. For example, mp.dps = 150 is required to get an accurate double precision result from mp_debye1(1e100, 1).

If method == 'nsum', mp.nsum is used to evalute the series form of the function.