sim_sir.jl

using MarkovProcesses 
using PyPlot

load_model("sir")

σ = simulate(SIR)
plt.figure()
plt.step(σ["times"], σ["I"], "ro--", marker="x", where="post", linewidth=1.0)
plt.savefig("sim_sir.png")
plt.close()