```{code-cell} python3

n = 5

P = 0.0125 * np.ones((n, n))

P = np.full((n, n), 0.0125)

P += np.diag(0.95 - 0.0125 * np.ones(5))

# State values of the Markov chain

s = np.array([0.95, 0.975, 1.0, 1.025, 1.05])

```{code-cell} python3

n = 5

P = 0.0125 * np.ones((n, n))

P = np.full((n, n), 0.0125)

P += np.diag(0.95 - 0.0125 * np.ones(5))

s = np.array([0.95, 0.975, 1.0, 1.025, 1.05]) # State values

mc = qe.MarkovChain(P, state_values=s)