# Best way to use/implement Debye function in Julia?

**URL:** <https://discourse.julialang.org/t/best-way-to-use-implement-debye-function-in-julia/111090>\
**Category:** General Usage\
**Tags:** numerics\
**Created:** [March 3, 2024, 12:36pm UTC](https://discourse.julialang.org/t/best-way-to-use-implement-debye-function-in-julia/111090 "2024-03-03T12:36:06Z")\
**Posts on this page:** 1\
**Showing post:** 16

<div class="post-metadata">

**Author:** ![stevengj](https://sea2.discourse-cdn.com/julialang/user_avatar/discourse.julialang.org/stevengj/32/71_2.png) [@stevengj](https://discourse.julialang.org/u/stevengj)\
**Post date:** [March 4, 2024, 9:16pm UTC](https://discourse.julialang.org/t/best-way-to-use-implement-debye-function-in-julia/111090/16 "2024-03-04T21:16:35Z")

</div>

> [@trg818](#):
>
> ```julia-auto
> for l in 0:k
> ll=l
> is=is+(-1)^l *(big(ll)^n)*binomial(k,l)
> end
> 
> ```

This is a slow way to implement a power series. You always want to compute each term as a recurrence from the next. See also [And Julia keeps amazing me: high precision computation - #6 by stevengj](https://discourse.julialang.org/t/and-julia-keeps-amazing-me-high-precision-computation/107740/6)

---

_[View the full topic](https://discourse.julialang.org/t/best-way-to-use-implement-debye-function-in-julia/111090)._
