def METHOD_NAME(self):
    """
    Test the get_heat_capacity() method.
    """
    t_list = np.array([300, 500, 1000, 1500, 2000])
    cv_exp_list = np.array([1, 1, 1, 1, 1])
    for temperature, cv_exp in zip(t_list, cv_exp_list):
        cv_act = get_heat_capacity(temperature, ENERGY_F, DEGENERACY_F, N0)
        assert round(abs(cv_exp / cv_act - 1.0), 4) == 0