*Debye functions.@ 2017-03-24 9:20 Ed Smith-Rowland2017-03-25 21:24 ` maxgacode 0 siblings, 1 reply; 3+ messages in thread From: Ed Smith-Rowland @ 2017-03-24 9:20 UTC (permalink / raw) To: gsl-discuss Greetings, I've been looking at the Debye integrals D_n(x) = \frac{n}{x^n}\int_{0}^{x} \frac{t^n}{e^t - 1}dt The integrand is everywhere positive. The definite integral must be zero for x=0. The values returned by gsl debye functions start at one for x=0 and monotonically decrease. The definite integral of a positive functions must start at zero and monotonically increase. Is it possible that we have a complementary Debye integral? Perhaps scaled? In any case, the functions can't match the formulas in the manual. Thank you. Ed Smith-Rowland ^ permalink raw reply [flat|nested] 3+ messages in thread

*Re: Debye functions.2017-03-24 9:20 Debye functions Ed Smith-Rowland@ 2017-03-25 21:24 ` maxgacode2017-03-27 19:04 ` Ed Smith-Rowland 0 siblings, 1 reply; 3+ messages in thread From: maxgacode @ 2017-03-25 21:24 UTC (permalink / raw) To: gsl-discuss Il 24/03/2017 10:19, Ed Smith-Rowland ha scritto: > Greetings, > > I've been looking at the Debye integrals > > D_n(x) = \frac{n}{x^n}\int_{0}^{x} \frac{t^n}{e^t - 1}dt > > The integrand is everywhere positive. > > The definite integral must be zero for x=0. But the 1/x factor goes to zero and so you get a 0/0 indeterminate ratio. Computing the limit to zero returns 1.0! > > The values returned by gsl debye functions start at one for x=0 and > monotonically decrease. Please note the factor \frac{n}{x^n} That factor is the responsible of the observed behavior. > > The definite integral of a positive functions must start at zero and > monotonically increase. > > Is it possible that we have a complementary Debye integral? Perhaps scaled? > > In any case, the functions can't match the formulas in the manual. > I don't think so. Please try to multiply the result of gsl_sf_debye_n(x) by n/x^n and see. Moreover the Chapter 27 of Abramowitz and Stegun (page 998 of my ninth edition) is listing the values of the Debye functions, you can easily verify that GSL implementation is correct. Hope this helps Max ^ permalink raw reply [flat|nested] 3+ messages in thread

*2017-03-25 21:24 ` maxgacodeRe: Debye functions.@ 2017-03-27 19:04 ` Ed Smith-Rowland0 siblings, 0 replies; 3+ messages in thread From: Ed Smith-Rowland @ 2017-03-27 19:04 UTC (permalink / raw) To: gsl-discuss On 03/25/2017 05:24 PM, maxgacode wrote: > Il 24/03/2017 10:19, Ed Smith-Rowland ha scritto: >> Greetings, >> >> I've been looking at the Debye integrals >> >> D_n(x) = \frac{n}{x^n}\int_{0}^{x} \frac{t^n}{e^t - 1}dt >> >> The integrand is everywhere positive. >> >> The definite integral must be zero for x=0. > > But the 1/x factor goes to zero and so you get a 0/0 indeterminate > ratio. Computing the limit to zero returns 1.0! > >> >> The values returned by gsl debye functions start at one for x=0 and >> monotonically decrease. > > > Please note the factor > > \frac{n}{x^n} > > That factor is the responsible of the observed behavior. > > >> >> The definite integral of a positive functions must start at zero and >> monotonically increase. >> >> Is it possible that we have a complementary Debye integral? Perhaps >> scaled? >> >> In any case, the functions can't match the formulas in the manual. >> > > I don't think so. Please try to multiply the result of > gsl_sf_debye_n(x) by n/x^n and see. > > Moreover the Chapter 27 of Abramowitz and Stegun (page 998 of my ninth > edition) is listing the values of the Debye functions, you can easily > verify that GSL implementation is correct. > > > Hope this helps > > Max > Ah. This is just a convention. Wolfram and others lose the n/x^n. So the thins look sigmoid and level off at \Gamma(n+1)\zeta(n+1). Sorry for the noise. Ed ^ permalink raw reply [flat|nested] 3+ messages in thread

end of thread, other threads:[~2017-03-27 19:04 UTC | newest]Thread overview:3+ messages (download: mbox.gz / follow: Atom feed) -- links below jump to the message on this page -- 2017-03-24 9:20 Debye functions Ed Smith-Rowland 2017-03-25 21:24 ` maxgacode 2017-03-27 19:04 ` Ed Smith-Rowland

This is a public inbox, see mirroring instructions for how to clone and mirror all data and code used for this inbox; as well as URLs for read-only IMAP folder(s) and NNTP newsgroup(s).