The high-temperature expansion of the thermal sunset

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.