This paper presents a Finite Element Method (FEM) application to thermal study of natural three layers of human dermal parts of varying properties. This paper carries out investigation of temperature distributions in these layers namely epidermis, dermis and under lying tissue layer. It is assumed that the outer skin is exposed to atmosphere and the loss of heat due to convection, radiation and evaporation of water have also been taken into account. The computations are carried out for one dimensional unsteady state case and the shape functions in dermal parts have been considered to be quadratic. A Finite Element scheme that uses the Crank-Nicolson method is used to solve the problem and the results computed have been exhibited graphically.