The main objective of this article is to develop a numerical algorithm based on Haar wavelet operational matrix method for solving the dual-phase-lag (DPL) bioheat transfer model under different boundary conditions during hyperthermia treatment. Unlike the existing operational matrix methods based on orthogonal functions, we construct the Haar wavelet operational matrices of integration without invoking the block pulse functions. The substantial importance of the proposed method lies in the fact that it converts the DPL bioheat transfer model into a system of algebraic equations which significantly speeds up the entire computational process. The efficiency and validity of the proposed method are demonstrated by comparing the numerical outcomes with the finite element wavelet Galerkin method and the exact solution of the model.