This paper presents a method for determining the temperature distribution in the cross-section of waterwall tubes connected by fins using an in-house numerical algorithm prepared in the MATLAB environment, based on differential equations with separable variables. In order to verify the correctness of the algorithm operation, the temperature values obtained from it, determined for the frontal area of the tubes, are compared with the temperatures found in the Ansys Fluent environment and those measured on the test stand. A system corresponding to a fragment of the combustion chamber wall of a supercritical steam boiler was selected to perform the analysis. The system consists of three tubes connected by fins. The temperature distributions in the cross-sections of the tubes were compared for the case when each of the tubes was heated on one side with the same heat flux and when the heat flux falling on the central tube was by 50% higher than the heat flux incident on the neighbouring tubes. Experimental verification was carried out on a stand equipped with three vertical tubes connected by fins, heated on one side by infrared radiators.