Plot of E_min now shows E_min[-1].

plots.py

@@ -32,16 +32,18 @@ def plot_energy(data, fig_path=""):
     fig, axs = plt.subplots(3, tight_layout=True, figsize=(16, 10))
 
     axs[0].plot(beta, E_min)
+    axs[0].plot(beta[-1], E_min[-1], "x", label=f"$E_{{min}}(\u03B2={beta[-1]:.2f}) = {E_min[-1]:.4f}$")
     axs[0].set(
         xlabel="\u03B2", ylabel="E / 1", title="Minimum Energy $E_{{min}}(\u03B2)$"
     )
     axs[0].grid()
+    axs[0].legend()
 
     axs[1].errorbar(beta, E_mean, yerr=2 * np.sqrt(E_var))
     axs[1].set(
         xlabel="\u03B2",
         ylabel="E / 1",
-        title="Mean Energy $E_{{mean}}(\u03B2) \u00B1 2\u03C3_{{mean}}$",
+        title="Mean Energy $E_{{mean}}(\u03B2) \u00B1 2\u03C3$",
     )
     axs[1].grid()