<p>This study presents a numerical investigation of a fractional-order malaria disease model using the Boubaker wavelet collocation method. The malaria model, formulated with fractional-order derivatives, is analyzed to capture the complex transmission dynamics between humans and vectors. The Boubaker wavelet approach, renowned for its efficiency and accuracy in solving FODEs, approximates the model's solutions. We constructed the integration operational matrices using the Boubaker wavelets. The Boubaker wavelet collocation method (BWCM) is employed to efficiently solve fractional-order systems. The effectiveness of the suggested approach for resolving fractional-order disease models is demonstrated by comparing numerical results obtained using NDSolve and Runge–Kutta methods. Compared to conventional numerical techniques, the BWCM requires fewer computer resources and delivers superior accuracy and efficiency. Sensitivity studies assess the impact of fractional-order parameters and other key variables on the disease dynamics. The study provides new insights into malaria transmission dynamics and highlights the importance of fractional-order models and advanced numerical techniques in epidemiological research.</p>