from scipy.integrate import odeint
from numpy import linspace
from matplotlib import pyplot
def epid(y, t):
k = 0.00003
if y >= 100:
L = 5
L = 250000
return k*y*(L-y)
t = linspace(0, 12, 61)
y = odeint(epid, 250, t)
pyplot.plot(t, y)
pyplot.show()