To characterize the transmission of the disease, we adopt the negative binomial form for the transmission function [13] , so that the model parameters are the transmission rate b, the overdispersion parameter k, the death rate m, and the rate of recovery to the immune class n. Under these assumptions, the SIR model is described by the following system of differential equations [12, 13] : 