@@ -227,7 +227,7 @@ def consumption_incomplete(cp, s_path):
227227 db = ((1 - β) * v - y) / β
228228229229 for i, s in enumerate(s_path):
230- c_path[i] = (1 - β) * (v - b_path[i] * np.ones((n, 1)))[s, 0]
230+ c_path[i] = (1 - β) * (v - np.full((n, 1), b_path[i]))[s, 0]
231231 b_path[i + 1] = b_path[i] + db[s, 0]
232232233233 return c_path, b_path[:-1], y[s_path]
@@ -255,7 +255,7 @@ fig, ax = plt.subplots(1, 2, figsize=(14, 4))
255255256256ax[0].set_title('Consumption paths')
257257ax[0].plot(np.arange(N_simul), c_path, label='incomplete market')
258-ax[0].plot(np.arange(N_simul), c_bar * np.ones(N_simul), label='complete market')
258+ax[0].plot(np.arange(N_simul), np.full(N_simul, c_bar), label='complete market')
259259ax[0].plot(np.arange(N_simul), y_path, label='income', alpha=.6, ls='--')
260260ax[0].legend()
261261ax[0].set_xlabel('Periods')
@@ -289,7 +289,7 @@ fig, ax = plt.subplots(1, 2, figsize=(14, 4))
289289290290ax[0].set_title('Tax collection paths')
291291ax[0].plot(np.arange(N_simul), c_path, label='incomplete market')
292-ax[0].plot(np.arange(N_simul), c_bar * np.ones(N_simul), label='complete market')
292+ax[0].plot(np.arange(N_simul), np.full(N_simul, c_bar), label='complete market')
293293ax[0].plot(np.arange(N_simul), y_path, label='govt expenditures', alpha=.6, ls='--')
294294ax[0].legend()
295295ax[0].set_xlabel('Periods')
@@ -584,7 +584,7 @@ class TaxSmoothingExample:
584584 plt.figure()
585585 plt.title('Tax collection paths')
586586 plt.plot(np.arange(N), self.T_path, label='incomplete market')
587- plt.plot(np.arange(N), self.T_bar * np.ones(N), label='complete market')
587+ plt.plot(np.arange(N), np.full(N, self.T_bar), label='complete market')
588588 plt.plot(np.arange(N), self.g_path, label='govt expenditures', alpha=.6, ls='--')
589589 plt.legend()
590590 plt.xlabel('Periods')