# representation

S1, K1 = kmuth.stationary_values()

# Extract scalars from the nested array

S1, K1 = S1.item(), K1.item()

# Form innovation representation state-space

Ak, Ck, Gk, Hk = A, K1, G, 1

simulate it

```{code-cell} python3

# Extract the scalar K1 from the nested array

K1 = K1[0][0]

# Create grand state-space for y_t, a_t as observed vars -- Use

# stacking trick above

Af = np.array([[ 1, 0, 0],