Outstanding post, Christopher!
It's intriguing that the shape of the infected curve in the Simple Example section is fairly symmetrical and does not have the long exponential tail.
Generating the phase plane analysis from the simulations is amazing. I've seen (and done) these analyses using the traditional calculus-based approach, but it's wonderful that you can do the same with this model and show that many of the same principles apply.
