The equation of state p(ǫ) and speed of sound squared c 2 s are studied in grand canonical ensemble of all hadron resonances having masses ≤ 2 GeV. This large ensemble is divided into strange and non-strange hadron resonances and furthermore to pionic, bosonic and femionic sectors. It is found that the pions represent the main contributors to c 2 s and other thermodynamic quantities including the equation of state p(ǫ) at low temperatures. At high temperatures, the main contributions are added in by the massive hadron resonances. The speed of sound squared can be calculated from the derivative of pressure with respect to the energy density, ∂p/∂ǫ, or from the entropy-specific heat ratio, s/cv. It is concluded that the physics of these two expressions is not necessarily identical. They are distinguishable below and above the critical temperature Tc. This behavior is observed at vanishing and finite chemical potential. At high temperatures, both expressions get very close to each other and both of them approach the asymptotic value, 1/3. In the HRG results, which are only valid below Tc, the difference decreases with increasing the temperature and almost vanishes near Tc. It is concluded that the HRG model can very well reproduce the results of the lattice quantum chromodynamics (QCD) of ∂p/∂ǫ and s/cv, especially at finite chemical potential. In light of this, energy fluctuations and other collective phenomena associated with the specific heat might be present in the HRG model. At fixed temperatures, it is found that c 2 s is not sensitive to the chemical potential.