diff --git a/examples/pincell_depletion/restart_depletion.py b/examples/pincell_depletion/restart_depletion.py
index 6a38b4bd9..fd0a4a8dd 100644
--- a/examples/pincell_depletion/restart_depletion.py
+++ b/examples/pincell_depletion/restart_depletion.py
@@ -72,21 +72,20 @@ time, Xe_capture = results.get_reaction_rate('1', 'Xe135', '(n,gamma)')
 ###############################################################################
 
 days = 24*60*60
-plt.figure()
-plt.plot(time/days, keff[:, 0], label="K-effective")
-plt.xlabel("Time (days)")
-plt.ylabel("Keff")
+fig, ax = plt.subplots()
+ax.errorbar(time/days, keff[:, 0], keff[:, 1], label="K-effective")
+ax.set_xlabel("Time [d]")
+ax.set_ylabel("Keff")
 plt.show()
 
-plt.figure()
-plt.plot(time/days, n_U235, label="U 235")
-plt.xlabel("Time (days)")
-plt.ylabel("n U5 (-)")
+fig, ax = plt.subplots()
+ax.plot(time/days, n_U235, label="U235")
+ax.set_xlabel("Time [d]")
+ax.set_ylabel("U235 atoms")
 plt.show()
 
-plt.figure()
-plt.plot(time/days, Xe_capture, label="Xe135 capture")
-plt.xlabel("Time (days)")
-plt.ylabel("RR (-)")
+fig, ax = plt.subplots()
+ax.plot(time/days, Xe_capture, label="Xe135 capture")
+ax.set_xlabel("Time [d]")
+ax.set_ylabel("Xe135 capture rate")
 plt.show()
-plt.close('all')
diff --git a/examples/pincell_depletion/run_depletion.py b/examples/pincell_depletion/run_depletion.py
index 1343e8b59..45dfd2a87 100644
--- a/examples/pincell_depletion/run_depletion.py
+++ b/examples/pincell_depletion/run_depletion.py
@@ -117,21 +117,20 @@ time, Xe_capture = results.get_reaction_rate('1', 'Xe135', '(n,gamma)')
 ###############################################################################
 
 days = 24*60*60
-plt.figure()
-plt.plot(time/days, keff[:, 0], label="K-effective")
-plt.xlabel("Time (days)")
-plt.ylabel("Keff")
+fig, ax = plt.subplots()
+ax.errorbar(time/days, keff[:, 0], keff[:, 1], label="K-effective")
+ax.set_xlabel("Time [d]")
+ax.set_ylabel("Keff")
 plt.show()
 
-plt.figure()
-plt.plot(time/days, n_U235, label="U235")
-plt.xlabel("Time (days)")
-plt.ylabel("n U5 (-)")
+fig, ax = plt.subplots()
+ax.plot(time/days, n_U235, label="U235")
+ax.set_xlabel("Time [d]")
+ax.set_ylabel("U235 atoms")
 plt.show()
 
-plt.figure()
-plt.plot(time/days, Xe_capture, label="Xe135 capture")
-plt.xlabel("Time (days)")
-plt.ylabel("RR (-)")
+fig, ax = plt.subplots()
+ax.plot(time/days, Xe_capture, label="Xe135 capture")
+ax.set_xlabel("Time [d]")
+ax.set_ylabel("Xe135 capture rate")
 plt.show()
-plt.close('all')