In this paper, an axisymmetric model of the human skin is developed to simulate the steady-state temperature distribution during contact with a hot solid. Simulations are carried out using the boundary element method. This study seeks to investigate the feasibility of using the boundary element method in the studies of burn. A sensitivity analysis is carried out to examine the effects of various parameters on the temperature distribution inside the skin during burn. Furthermore, a statistical analysis based on the Taguchi method is performed to determine the combination of factors that produce the desired outcome (least increase in temperature). In order to validate the accuracy of the numerical scheme, results obtained using the boundary element method are compared with the solutions obtained using the more established finite-element method.