Fast calculation or simulation of the ablation zone induced by radiofrequency ablation (RFA) has a critical role in hepatic RFA planning and therapy. However, it remains challenging to approximate the ablation zone in real time, especially when more than one probe is involved in one ablation session. This paper presents a novel computational technique to calculate the 3D ablation zone of one probe RFA and two-probe switching RFA. The main idea is to get an approximate solution of the temperature distribution from a simplified Pennes bioheat equation, and further fit the solution to the coagulation measurements on ex vivo porcine liver. With a closed-form solution of temperature distribution, the calculation of the ablation zone is as simple as the commonly used ellipsoidal model, but it allows a more realistic prediction of combined ablation zones with different inter-probe spacing. The new approximation technique could potentially replace the original ellipsoidal model in the intervention planning step.