From 40bfac372fb7f18728e3b195d809f7d4755ae72c Mon Sep 17 00:00:00 2001
From: Paul CACHEUX <paulcacheux@gmail.com>
Date: Mon, 11 Jun 2018 08:54:45 +0200
Subject: [PATCH] Add new figure
---
SRI_matlab/main.m | 4 ++--
SRI_matlab/step.m | 6 +++---
2 files changed, 5 insertions(+), 5 deletions(-)
diff --git a/SRI_matlab/main.m b/SRI_matlab/main.m
index f520f4c..fb9df2e 100644
--- a/SRI_matlab/main.m
+++ b/SRI_matlab/main.m
@@ -14,11 +14,11 @@ coeffs.ImmunisationLossRate = 2.7e-4;
coeffs.DeathRate = 1.3e-2;
coeffs.NaturalDeathRate = 2e-4;
coeffs.NatalityRate = 5.3e-4;
-coeffs.QuarantineRate = 0.01;
+coeffs.QuarantineRate = 0;
VaccinationVector = build_vacc_rate(0, 0, 0, N);
% VaccinationVector = build_vacc_rate(0, 0, 50, N);
-m = model(S, I, Q, R, 12, N);
+m = model(S, I, Q, R, 1, N);
t = 1:N;
diff --git a/SRI_matlab/step.m b/SRI_matlab/step.m
index 2353e56..c874baf 100644
--- a/SRI_matlab/step.m
+++ b/SRI_matlab/step.m
@@ -12,9 +12,9 @@ for i=6:length(X)
end
if length(X) > 4
- % infection = S * (1 - 25 * options.InfectionRate / N) ^ I;
- infecter = sum(X(5:length(X)) + I);
- infection = options.InfectionRate * S * infecter/N;
+ infection = S * (1 - 25 * options.InfectionRate / N) ^ I;
+ % infecter = sum(X(5:length(X)) + I);
+ % infection = options.InfectionRate * S * infecter/N;
new_waiting(1) = infection;
dS = -infection;
dI = X(length(X));
--
GitLab