Statsmodels VARMAX: confidence / predication intervals with more than one endogenous variable

Question

I am trying to recover confidence/prediction intervals in Python Statsmodels (Version 0.12.1) with two or more endogenous (y) variables, as is common in VARMAX. The followin