SIR-Model-ODE.jl

# SIR-Model-ODE.jl
# Juhani Kaukoranta 06.03.2020

using DifferentialEquations
using Plots

function sir!(du,u,p,t)
	du[1] = -p[1]*u[1]*u[2]
	du[2] = p[1]*u[1]*u[2] - p[2]*u[2]
	du[3] = p[2]*u[2]
end
u0 = [1.0, 1.27e-6, 0.0] # alkuarvot sir
tspan = (0.0, 150.0) # aikaväli 150 päivää
p = [1.0/2,1.0/3]  # parametrit b ja k
prob = ODEProblem(sir!,u0,tspan,p)
sol = solve(prob)
plot(sol,title = "SIR-ODE- model",label = ["s" "i" "r"])