In recent years, the availability of highly pure stable isotopes has made possible the investigation of the dependence of the physical properties of crystals, in particular semiconductors, on their isotopic composition. Following the investigation of the specific heat ($C_p$, $C_v$) of monatomic crystals such as diamond, silicon, and germanium, similar investigations have been undertaken for the tetrahedral diatomic systems ZnO and GaN (wurtzite structure), for which the effect of the mass of the cation differs from that of the anion. In this article we present measurements for a semiconductor with rock salt structure, namely lead sulfide. Because of the large difference in the atomic mass of both constituents ($M_{rm Pb}$= 207.21 and ($M_{rm S}$=32.06 a.m.u., for the natural isotopic abundance) the effects of varying the cation and that of the anion mass are very different for this canonical semiconductor. We compare the measured temperature dependence of $C_p approx C_v$, and the corresponding derivatives with respect to ($M_{rm Pb}$ and $M_{rm S}$), with textit{textit{ab initio}} calculations based on the lattice dynamics obtained from the local density approximation (LDA) electronic band structure. Quantitative deviations between theory and experiment are attributed to the absence of spin-orbit interaction in the ABINIT program used for the electronic band structure calculations.