cpr_sns: bigger alpha

examples/cpr_sns.py

@@ -78,7 +78,7 @@ def j(T, h, phi, n=15, eta=0.1, alpha_soc=0.1, L=10):
 
 def main():
     T = 0.1
-    alpha_soc = 0.3
+    alpha_soc = 0.9
     h = np.r_[0.0, 0.25, 0.5, 1.0, 1.5]
     phi = np.linspace(-pi, pi, 37)
 
@@ -91,7 +91,7 @@ def main():
     plt.plot(phi / pi, Jx_mean.T)
     plt.xlabel(r"$\varphi / \pi$")
     plt.ylabel(r"$I$")
-    plt.legend()
+    plt.legend(h, title="$h$")
     plt.savefig("cpr_sns.pdf")