Update figure comparing Debye and Einstein models. Added experimental data for silver.

src/2_debye_model.md

@@ -123,31 +123,20 @@ temp = np.linspace(1, 215, 100)
 
 fit = curve_fit(c_einstein, T, c, 500)
 T_E = fit[0][0]
-#delta_T_E = np.sqrt(fit[1][0, 0])
-#print(f"T_E = {T_E:.5} ± {delta_T_E:.3} K")
 
 fit = curve_fit(c_debye, T, c, 500)
 T_D = fit[0][0]
-#delta_T_D = np.sqrt(fit[1][0, 0])
-#print(f"T_D = {T_D:.5} ± {delta_T_D:.3} K")
 
 fig, ax = pyplot.subplots()
 ax.scatter(T, c)
 ax.set_title('Heat capacity of silver compared to the Debye and Einstein models')
-ax.plot(temp, c_einstein(temp, T_E), label='Einstein model')
-ax.plot(temp, c_debye(temp, T_D), label='Debye model')
-ax.set_ylim(bottom=0, top=3.4)
+ax.plot(temp, c_debye(temp, T_D), label=f'Debye model, $T_D={T_D:.5}K$')
+ax.plot(temp, c_einstein(temp, T_E), label=f'Einstein model, $T_E={T_E:.5}K$')
+ax.set_ylim(bottom=0, top=3)
 ax.set_xlim(0, 215)
 ax.set_xlabel('$T(K)$')
 ax.set_ylabel(r'$C/k_B$')
-ax.set_xticks([T_E, T_D])
-ax.set_xticklabels(['$T_E$','$T_D$'])
-ax.set_yticks([3])
-ax.set_yticklabels(['$3$'])
-ax.legend(loc='upper left')
-pyplot.hlines([3], 0, 215, linestyles='dashed')
-pyplot.vlines([T_E,T_D], 0, 3.4, linestyles='dashed')
-#draw_classic_axes(ax, xlabeloffset=0.3)
+ax.legend(loc='lower right');
 ```
 
 ## Exercises