Mathematical modelling is the theoretically established tool for developing solar heating systems, e.g. with using transfer functions. If we know the transfer functions of the system, the outlet temperature can be predicted as a function of the input variables (solar irradiance, inlet temperature, environment temperatures), dynamic analysis can be carried out, and stable system control can be effectively designed based on the well-tried methods of control engineering. For these purposes, new, validated transfer functions for solar heating systems are worked out in this study based on a mathematical model, which can be found in the literature and has been applied successfully in the field. The transfer functions are used for dynamic analysis and control design of solar heating systems. The dynamic analysis is presented and the efficiency of the proposed stable control is demonstrated with respect to a real solar heating system. (C) 2014 Elsevier Ltd. All rights reserved.