sim_sir.jl 209 Bytes
Newer Older
1
2
3
4

using MarkovProcesses 
using PyPlot

5
load_model("sir")
6
7
8

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