Parabolic heat equation based on Fourier's theory (FHE), and hyperbolic heat equation (HHE), has been used to mathematically model the temperature distributions of biological tissue during thermal ablation. However, both equations have certain theoretical limitations. The FHE assumes an infinite thermal energy propagation speed, whereas the HHE might possibly be in breach of the second law of thermodynamics. The relativistic heat equation (RHE) is a hyperbolic-like equation, whose theoretical model is based on the theory of relativity and which was designed to overcome these theoretical impediments. In this study, the three heat equations for modelling of thermal ablation of biological tissues (FHE, HHE and RHE) were solved analytically and the temperature distributions compared. We found that RHE temperature values were always lower than those of the FHE, while the HHE values were higher than the FHE, except for the early stages of heating and at points away from the electrode. Although both HHE and RHE are mathematically hyperbolic, peaks were only found in the HHE temperature profiles. The three solutions converged for infinite time or infinite distance from the electrode. The percentage differences between the FHE and the other equations were larger for higher values of thermal relaxation time in HHE.