Modeling of a natural circulation boiler for a coal-fired thermal power station is presented here. The boiler system is divided into seven subcomponents, and for each section, models based on conservation of mass, momentum, and energy are formulated. The pressure drop at various sections and the heat transfer coefficients are computed using empirical correlations. Solutions are obtained by using SIMULINK. The model is validated by comparing its steady state and dynamic responses with the actual plant data. Open loop responses of the model to the step changes in the operating parameters, such as pressure, temperature, steam flow, feed water flow, are also analyzed. The present model can be used for the development and design of effective boiler control systems.