using MarkovProcesses

load_model("SIR")
SIR.time_bound = 100.0
σ1, σ2 = simulate(SIR), simulate(SIR)

dist_lp([σ1], [σ2])

return true