How can I incorporate 95% confidence band limits in a recursive forecasting setup?

[dmresearch15]

Specifically, I want to implement 95% confidence intervals for both a single series (e.g., M750) and a panel of series (e.g., M1, M2, M750, and M1000) using the example provided in the Modeltime recursive forecasting guide.

What steps should I follow to achieve this?