The behavior of the processes taking place in furnaces determines the efficiency of fuel chemical energy utilization, the quality of the final products and the environmental safety of the production. Mathematical models of the processes of gas dynamics and heat transfer in the working space of heating equipment are quite complex, and do not allow the establishment of a direct analytical relationship between the quality indicator of the process (F) and the influencing parameters (xi). To simplify the procedure for obtaining the values of the function F depending on the change of parameters xi, a method of successive approximation is presented in the article. The main idea of the method is that the representation of the function around a point from the domain of the function can be extended to the entire domain for many problems of mechanics. The relative error in the definition of the function acquires its maximum value at the border of the area, and a reasonable narrowing of it allows control of the size of the error. Thus, the advantages of using the method are obvious; it is able to provide approximation of the function in a multiplicative form with a controlled error. The distribution of the method to the field of heat transfer problems is presented in this paper. The successful implementation of this method for solving problems of this kind shows that the solution of practical problems may be generalized for the entire domain of the function, despite the fact that the errors of such a representation increase to 5–7% when approaching its limit, which, however, may be considered acceptable for engineering calculations.