parameter tuning

code/heat_capacity.py

@@ -33,7 +33,7 @@ def c_debye(T, T_D=1):
 
 
 def plot_einstein():
-    T = np.linspace(0.01, 3, 500)
+    T = np.linspace(0.01, 1.5, 500)
     fig, ax = pyplot.subplots()
 
     ax.plot(T, c_einstein(T)/3)
@@ -46,7 +46,7 @@ def plot_einstein():
     ax.set_xticklabels([r'$\hbar \omega/k_B$'])
     ax.set_yticks([1])
     ax.set_yticklabels(['$k_B$'])
-    pyplot.hlines([1], 0, 3, linestyles='dashed')
+    pyplot.hlines([1], 0, 1.5, linestyles='dashed')
     draw_classic_axes(ax)
     pyplot.savefig('zeropoint.svg')