How to fit multivariate normal distribution with autocorrelation to data in Python?

Question

I have several data series. They both covary with each other and are autocorrelated with themselves. I\'d like to fit a multivariate normal distribution to these data, while als