Abstract: The mathematical model is established for multiple layer cylindrical walls. The iteration and inversion methods are used for Laplace transformation by application of the Laplace transf ormation and Stehfest numerical inversion method. The temperature profile is obtained. The room measurem ent results verify that the solution method is at the high precision. The solution method is also used under different initial and boundary conditions with different heat source systems. The solution method can be applied to several engineering fields.