The phenomenon of thermal ablation is described by Pennes' bioheat equation. This model is based on Newton's law of cooling. Many approximate methods have been considered because of the importance of this issue. We propose an implicit numerical scheme which has better stability properties than other approaches.