In the field of bio heat transfer, as of now, the main concern of researchers lies in the proper and accurate thermal damage of the diseased tissues without destroying or damaging the neighboring healthy tissues during the tumor treatment. The present work aims to develop a new approach toward solving the bio heat transfer equations for the skin burn and hyperthermia treatments. Both analytical and numerical solutions are proposed. For the analytical study, a differential transform method is used to solve steady and unsteady state heat equations. The finite volume method is adopted to solve these equations numerically, which provides a better scope to solve the highly nonlinear complex equations. To obtain a complete solution, a code is developed in MATLAB and MATHEMATICA. The variation of different parameters, such as perfusion constant, space heating, surface step heating, and thermal conductivity, with time were observed. Apart from the above analysis of temperature distribution during skin burn through the spilling of hot beverage, its numerical solution was also performed for this problem at different boundary conditions. It was observed that with the help of the temperature distribution, depending on the time and the severity of the burn, different ranges of depths of the burn can be determined.