  We give a prescription for calculating the high-temperature expansion of the
thermal sunset integral to arbitrary order. We derive all terms odd in $T$, and
rederive previous results up to $\mathcal{O}(T^0) $ for both bosonic and
fermionic thermal sunsets in dimensional regularisation. We perform analytical
and numerical cross-checks. Intermediate steps involve integrals over three
Bessel functions.
