diff --git a/figures/SchueckerSchmidt2017/Fig9_simulation.eps b/figures/SchueckerSchmidt2017/Fig9_simulation.eps
index 6ae2d4d270fb52999a1b15e133992de26751a64f..eed822c92ba0db94e8d41af4370b4b4e632cd9f4 100644
--- a/figures/SchueckerSchmidt2017/Fig9_simulation.eps
+++ b/figures/SchueckerSchmidt2017/Fig9_simulation.eps
@@ -3,7 +3,7 @@
 %%HiResBoundingBox: 0.000000 0.000000 374.400000 231.391925
 %%Creator: GPL Ghostscript 922 (ps2write)
 %%LanguageLevel: 2
-%%CreationDate: D:20180529165002+09'00'
+%%CreationDate: D:20180529173051+09'00'
 %%EndComments
 %%BeginProlog
 save
diff --git a/figures/SchueckerSchmidt2017/Fig9_simulation.py b/figures/SchueckerSchmidt2017/Fig9_simulation.py
index b907d884b1f9c549ba9fc49bcab28a40ebaf3372..38245aa1d451e36c1a3e694d19787f7a78b37917 100644
--- a/figures/SchueckerSchmidt2017/Fig9_simulation.py
+++ b/figures/SchueckerSchmidt2017/Fig9_simulation.py
@@ -50,12 +50,13 @@ network_params = {'N_scaling': 1.,
 ana_spec = {'load_areas': []}
 M = MultiAreaModel(network_params,
                    analysis=True,
+                   ana_spec=ana_spec,
                    simulation=True,
                    sim_spec=sim_params,
                    theory=True,
                    theory_spec=theory_params)
 p, r_stab = M.theory.integrate_siegert()
-M.analysis.create_pop_rates()
+M.analysis.create_pop_rates(areas=M.area_list)
 
 
 """
@@ -83,8 +84,6 @@ for i, area in enumerate(area_list):
 matrix = np.transpose(matrix)
 rate_matrix_plot(panel_factory.figure, ax, matrix, position='left')
 
-print("Meanfield: Area-averaged population rates: ", np.mean(matrix, axis=1))
-
 
 """
 Panel B
@@ -111,8 +110,6 @@ for i, area in enumerate(M.area_list):
 matrix_sim = np.transpose(matrix_sim)
 rate_matrix_plot(panel_factory.figure, ax, matrix, position='right')
 
-print("Area-averaged population rates: ", np.mean(matrix, axis=0))
-
 
 """
 Panel C