We provide an overview of the mathematical modeling of deterministic and stochastic infectious disease models. These models enable one to understand the outbreak, spread, and extinction of disease. We then focus on stochastic models with a disease reservoir to understand outbreak vulnerability for zoonotic diseases such as Ebola Virus Disease (EVD). Numerical results from a more complicated EVD model are compared with the theoretical results of a simplified stochastic SISk model. We also demonstrate the effect that vaccine has on outbreak vulnerability in a population that is connected to a disease reservoir