DOI: 10.36347/sjpms.2023.v10i01.002

Pages: 11-42

Background: In this paper, an expanded 10-Dimensional deterministic mathematical dynamic model was formulated that accounted for the role of global stability analysis in the methodological application of dual-bilinear controls with vaccination and impeccable role of adaptive immune response in the control of COVID-19 in Nigeria. In reality, following the introduction of both nonpharmaceutical and pharmacotherapy and the recent availability of vaccines for the control and treatment of the deadly aerosol viral load known as COVID-19, a number of notable scientific investigations on the transmission and treatment dynamics have been conducted but without thoughtful contributions on the combination of these multi-facet control functions that could lead to feasible eradication of the deadly virus. Methods: The model was formulated based expanded 10-Dimensional deterministic dynamic mathematical subpopulations with compartmental interactions investigated using triple-bilinear control functions: bilinear nonpharmaceutical (face-masking and social distancing - , ), bilinear pharmacotherapies (hydroxylchloroquine and azithromycin - , ) and bilinear immunity controls (adaptive immune effectors and BNT162b2 vaccine - , ). Experimental Data was collected from University of Calabar Teaching Hospital from the period July, 2022 through September, 2022, as the initial and final time intervals. Apart from fundamental theory of differential equations explored for system mathematical properties, analytical predictions explored classical method of Lyapunov functions with the incorporation of the theory of Volterra-Lyapunov stable matrices for the analysis of the system global stability conditions. Results: System mass actions and the reproduction numbers for both off-treatment and onset-treatment scenarios was for the first time computed with explicit results obtained ( and ). Moreso, off-treatment scenario showed that population extinction was eminent following the unabated ....