diff --git a/examples/slds_example_figure.py b/examples/slds_example_figure.py
index a276805..933b436 100644
--- a/examples/slds_example_figure.py
+++ b/examples/slds_example_figure.py
@@ -97,6 +97,7 @@ def random_dynamics(n):
     z = np.repeat(z, D)
 
     statesobj = slds._states_class(model=slds, T=z.size, stateseq=z, inputs=np.ones((z.size, 1)))
+    statesobj.generate_gaussian_states()
     y = statesobj.data = statesobj.generate_obs()
     x = statesobj.gaussian_states
     slds.states_list.append(statesobj)