Mathematical analysis of swine influenza epidemic model with optimal control

A deterministic model is designed and used to analyze the transmission dynamics and the impact of antiviral drugs in controlling the spread of the 2009 swine influenza pandemic. In particular, the model considers the administration of the antiviral both as a preventive as well as a therapeutic agent. Rigorous analysis of the model reveals that its disease-free equilibrium is globally asymptotically stable under a condition involving the threshold quantity-reproduction number Rc\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathcal {R}}_c$$\end{document}. The disease persists uniformly if Rc>1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathcal {R}}_c>1$$\end{document} and the model has a unique endemic equilibrium under certain condition. The model undergoes backward bifurcation if the antiviral drugs are completely efficient. Uncertainty and sensitivity analysis is presented to identify and study the impact of critical model parameters on the reproduction number. A time dependent optimal treatment strategy is designed using Pontryagin’s maximum principle to minimize the treatment cost and the infected population. Finally the reproduction number is estimated for the influenza outbreak and model provides a reasonable fit to the observed swine (H1N1) pandemic data in Manitoba, Canada, in 2009.