Compute the stationary values for W

weather.py

"""Licensed under GPLv3 <http://www.gnu.org/licenses/>.
Authors:
Eric Lopes
"""
import scipy.linalg
import numpy

W = numpy.matrix(
    [[0.4 ,0.4 ,0.2 ,0.0],
     [0.2 ,0.2 ,0.3 ,0.3],
     [0.25,0.25,0.25,0.25],
     [0.1 ,0.4 ,0.2 ,0.3]])

evals, lvecs = scipy.linalg.eig(W, right=False, left=True)
tolerance = 1e-15
mask = abs(evals - 1) < tolerance
evals = evals[mask]
lvecs = lvecs[:, mask]
equilibrium = lvecs/lvecs.sum(axis=0, keepdims=True)
print(equilibrium)