examples: fixup

examples/bulk_j_an.py

@@ -109,16 +109,20 @@ j_anom = np.vectorize(
 
 
 def main():
-    T = 0.05
-    h = np.linspace(0, 1.75, 21)
+    T = 0.02
+    h = np.linspace(0, 1.75, 29)
 
     j, j_an = j_anom(T, h, perturbative=False)
 
     t3 = np.diag([1, 1, -1, -1])
     j = tr(j @ t3)[:, 2, 2, UP]
 
-    plt.plot(h, j.real, "-", label=r"numerics")
-    plt.plot(h, j_an.real, "--", label=r"analytic")
+    plt.plot(h, -j.real, "-", label=r"numerics")
+    plt.plot(h, -j_an.real, "--", label=r"analytic")
+    plt.title(rf"$T = {T} \Delta$")
+    plt.xlabel(r"$h/\Delta$")
+    plt.ylabel(r"$j_{\mathrm{an}}$")
+    plt.legend()
     plt.savefig("bulk_j_an.pdf")