From 7d5be0d795cfad2b79442fe57b77d6fcefb954e8 Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Rivny=C3=A1k=20T=C3=ADmea?= <rivnyak.timea@hallgato.ppke.hu> Date: Tue, 1 Jun 2021 18:00:04 +0000 Subject: [PATCH] deriv function for SIRC --- deriv.m | 12 ++++++++++++ 1 file changed, 12 insertions(+) create mode 100644 deriv.m diff --git a/deriv.m b/deriv.m new file mode 100644 index 0000000..6e7dad9 --- /dev/null +++ b/deriv.m @@ -0,0 +1,12 @@ +function dydt = deriv(y, nu, mu, epsilon, beta, gamma, Gamma, q) + S = y(1); + I = y(2); + R = y(3); + C = y(4); + dSdt = nu - (beta*I + epsilon*beta*C)*S - mu*S; + dIdt = (beta*I + epsilon*beta*C)*S - gamma*I - mu*I; + dRdt = gamma*(1-q)*I + Gamma*C - mu*R; + dCdt = gamma*q*I - Gamma*C - mu*C; + dydt = [dSdt; dIdt; dRdt; dCdt]; +end + -- GitLab