I am trying to simulate the transmission of viruses in a population using the function ode from the deSolve package. The basic of my model is a SIR
ode
deSolve