The goal of this thesis is to study the dynamics behind the evolution of virulence. We examine first the underlying mechanics of linear systems of ordinary differential equations by investigating the classification of fixed points in these systems, then applying these techniques to nonlinear systems. We then seek to establish the validity of a system that models the population dynamics of uninfected and infected hosts---first with one parasite strain, then n strains. We define the basic reproductive ratio of a parasite, and study its relationship to the evolution of virulence. Lastly, we investigate the mathematics behind superinfection.