compute_stationary_variance function

Compute the stationary variance matrix