specific_heat function

Calculate specific heat