The higher heating value (HHV) is the most important indicator of a coal's potential energy yield. It is commonly used in the efficiency and optimal design calculations pertaining to the coal combustion and gasification processes. Since the experimental determination of coal's HHV is tedious and time-consuming, a number of proximate and/or ultimate analyses based correlations—which are mostly linear—have been proposed for its estimation. Owing to the fact that relationships between some of the constituents of the proximate/ultimate analyses and the HHV are nonlinear, the linear models make suboptimal predictions. Also, a majority of the currently available HHV models are restricted to the coals of specific ranks or particular geographical regions. Accordingly, in this study three proximate and ultimate analysis based nonlinear correlations have been developed for the prediction of HHV of coals by utilizing the computational intelligence (CI) based genetic programming (GP) formalism. Each of these correlations possesses following noteworthy characteristics: (i) the highest HHV prediction accuracy and generalization capability as compared to the existing models, (ii) wider applicability for coals of different ranks and from diverse geographies, and (iii) structurally lower complex than the other CI-based existing HHV models. It may also be noted that in this study, the GP technique has been used for the first time for developing coal-specific HHV models. Owing to the stated attractive features, the GP-based models proposed here possess a significant potential to replace the existing models for predicting the HHV of coals.