The current manuscript investigates the dynamics of a fractal–fractional malaria disease model caused by the Plasmodium parasite, in the framework of Caputo operator. The model under consideration has six compartments, namely, four categories of humans and two of vector populations. By utilizing the fixed point theory, we establish the existence results and uniqueness of solution for the aforementioned model. The Ulam–Hyers (UH) approach is carried out to evaluate the stability of the considered system. For the approximate solution of the model, the fractional Adams–Bashforth iterative techniques have been employed to obtain the numerical solution for the specified system. To simulate the model behavior, we use various values for the fractional order