The technique of integral Hankel transform to find the solution of heat conduction in half- -space coated by a multilayered package of homogenous laminae is applied. The half-space is heated by the given heat flux on the boundary surface. The temperature and heat flux distribution in the radial direction is analyzed for two types of coatings: 1) when the heat conductivity coefficient is described by a power or exponential function of the distance to the boundary surface; 2) multilayered coating has a periodic structure.