Modified Sir Model

I am making a modified SIR model with an added vaccination parameter V. InitIALLY all the nodes in the graph are susceptible and there are a few initially infected people. The init

Solution 1:

Your problem is that V is only measuring the total number of individuals becoming vaccinated in that time step, rather than the cumulative number of vaccinations.

So S+I+V+R is not constant. If you want V to be the cumulative number of vaccinations, then do V.append(V[-1]+len(vaccinated))

I'm not sure that this test if len(vaccinated)< (Vl*N) is doing what you think as well. It is checking whether the number of people vaccinated in a given time step is less than some fraction of the entire population. I suspect you want to use the cumulative number of vaccinations.