This paper introduces a mathematical model that includes four categories of hosts and five categories of vectors. The study investigates three control function strategies: personal protection, larval reduction, and adult mosquito vector reduction. To analyze the stability of the model, disease-free equilibria (DFE) were obtained by setting the adult population to zero, while endemic equilibria (EE) were identified at a fixed value. The basic reproduction number \(({\mathcal {R}}_0)\) was specified, and the necessary equations for its computation were derived to evaluate the stability of the model. The analysis utilizing the Jacobian matrix (JM) revealed that both DFE and EE exhibit stability when the value of \({\mathcal {R}}_0<1\) and instability when the value of \({\mathcal {R}}_0\ge 1\) . Furthermore, a Lyapunov function (LF) was formulated to demonstrate the overall stability of these equilibrium locations. The results suggest that the DFE is globally stable for \({\mathcal {R}}_0\le 1\) . Furthermore, when \({\mathcal {R}}_0>1\) , there is a single unstable equilibrium (EE), indicating that the disease persists. The EE is stable under specific parameter conditions. Numerical simulations demonstrate that effectively implementing control functions \(u_1,~ u_2\) , and \(u_3\) substantially reduces the infected human and adult mosquito populations and maintains the basic reproduction number \({\mathcal {R}}_0\) below 1. Additional simulations reveal that the parameter \(c\) plays a critical role in infection growth. Infections increase exponentially when \(c\le 0.0006\) , whereas the development of infected humans is regulated when \(c>0.0006\) . The results highlight how crucial it is to combine vector control and treatment tactics with meticulous parameter management to reduce and possibly completely eradicate the disease transmission.
