First function: SolveSimpleSIR