Add conservation test: (s + i + r == population) for all days simulated. Check non-negative s, i, r for all days simulated.