Cumulative distribution function of 31st January 2019’s daily temperature (left) and expected daily HDD during January 2019 days computed 30 days ahead of the month (right) by the FFT method (smooth and plain curve) and Monte-Carlo (95% confidence interval drawn with dashed lines). The vertical line on the left corresponds to |$T_{b}$|.