Implement heat capacity in free_fermi_proj.f90 via numerical derivative