pub fn sir_petri(th: Rc<ModalDblTheory>) -> ModalDblModel
The SIR model viewed as a reaction network.