The Debye function

Mirrors the C# Debye.Function: a piecewise approximation of D(x) = (n/x^n) * integral_0^x t^n / (e^t - 1) dt, vectorized over x.

Usage

debye(x)

Arguments

x

a non-negative numeric vector.

Value

a numeric vector the same length as x.

Examples

debye(0.5)
#> [1] 0.8249631
debye(c(0.1, 1, 10))
#> [1] 0.96299994 0.67441569 0.01929577