diff --git a/examples/tutorials/starting/analysis_2.py b/examples/tutorials/starting/analysis_2.py
index e6e003ddbe..7949e211f3 100644
--- a/examples/tutorials/starting/analysis_2.py
+++ b/examples/tutorials/starting/analysis_2.py
@@ -394,8 +394,10 @@
 #
 
 energy_bounds = [1, 10] * u.TeV
-spec.plot(energy_bounds=energy_bounds, energy_power=2)
-ax = spec.plot_error(energy_bounds=energy_bounds, energy_power=2)
+
+fig, ax = plt.subplots(figsize=(8, 6))
+spec.plot(ax=ax, energy_bounds=energy_bounds, sed_type="e2dnde")
+spec.plot_error(ax=ax, energy_bounds=energy_bounds, sed_type="e2dnde")
 plt.show()
 
 
@@ -417,6 +419,8 @@
 # %%time
 flux_points = fpe.run(datasets=[stacked])
 
-ax = spec.plot_error(energy_bounds=energy_bounds, energy_power=2)
-flux_points.plot(ax=ax, energy_power=2)
+fig, ax = plt.subplots(figsize=(8, 6))
+spec.plot(ax=ax, energy_bounds=energy_bounds, sed_type="e2dnde")
+spec.plot_error(ax=ax, energy_bounds=energy_bounds, sed_type="e2dnde")
+flux_points.plot(ax=ax, sed_type="e2dnde")
 plt.show()