Mathematical study for Zika virus transmission with general incidence rate

Abstract

An appropriate mathematical model for describing the Zika virus transmission with nonlinear general incidence rate was proposed. The basic reproduction number $ \mathcal{R}_0 $ was calculated using the next generation matrix method. Analysis of the local and the global stability of the equilibrium points was detailed using Jacobian linearisation method and Lyapunov theory, respectively. We proved that the disease-free equilibrium is locally and globally asymptotically stable when $ \mathcal{R}_0 $ is small than 1, and the infected equilibrium point is locally and globally asymptotically stable when $ \mathcal{R}_0 $ is greater than 1. The overall sensitivity analysis is based on statistical tools. This method consists of varying the parameters of the model to study one by one and then observe the effect of this variation on the model output. Sensitivity indices quantifying the influence of parameters on the output, always depend on the sample list of parameters. Later, we used optimal control to examine the effect of treatment where the purpose is to minimize the number of infected individuals with optimal treatment cost by applying Pontryagin's maximum principle. Therefore, we formulated an optimal control problem using the most parameter that influences the model output as a control parameter. The existence of the solution was proved and characterized using adjointt variables. Finally, a numerical scheme was applied to solve the coupled systems. Obtained results are validated numerically.

Keywords

• zika virus
• global stability
• lyapunov function
• lasalle's invariance principle
• sensitivity analysis
• nonlinear incidence rate
• optimal control.