Solver and ODE model function

Revised code and figures. All variables treated as normalized.
Oscillations can be seen but contain negative numbers. S and R1 grow
significantly, with R1 taking negative values.