An analytical solution given by Bessel series to the transient and one-dimensional (1D) bioheat equation in a multilayer region with spatial dependent heat sources is derived. Multilayer regions with 1D Cartesian, cylindrical or spherical geometries and composed of different types of biological tissues characterised by temperature-invariant physiological parameters are considered. Boundary conditions of first, second and third kinds to the temperature at the inner and outer surfaces are also assumed. In this work, the bioheat transfer model is applied to obtain the temperature profiles in a tumor bed and a surrounding healthy tissue using two spatial dependent heat source terms to simulate a magnetic fluid hyperthermia technique in the cancer treatment. The influence of these two heat sources, described by polynomial and exponential functions, on temperature is investigated.