From 3f90693e524e5a6330af45831085e2ead3f91fb0 Mon Sep 17 00:00:00 2001
From: Laura Mieczkowska <s180881@student.pg.edu.pl>
Date: Sun, 19 Jun 2022 12:24:47 +0000
Subject: [PATCH] Program- FINALNA WERSJA
---
...zprzestrzeniania_sie_infekcji-_program.cpp | 95 +++++++++++++++++++
1 file changed, 95 insertions(+)
create mode 100644 Modele_rozprzestrzeniania_sie_infekcji-_program.cpp
diff --git a/Modele_rozprzestrzeniania_sie_infekcji-_program.cpp b/Modele_rozprzestrzeniania_sie_infekcji-_program.cpp
new file mode 100644
index 0000000..af31925
--- /dev/null
+++ b/Modele_rozprzestrzeniania_sie_infekcji-_program.cpp
@@ -0,0 +1,95 @@
+/*
+Projekt specjalnosciowy
+
+"Modele rozprzestrzeniania sie infekcji"
+Mirella Narewska, Laura Mieczkowska, Konrad Korejwo
+
+*/
+
+
+
+#include<iostream>
+#include<fstream>
+
+using namespace std;
+
+float dhdt(float N, float bHN, float adH, float t, float C, float H) {
+
+
+ return (H * (bHN * (N - C - H) - adH));
+}
+
+float dcdt(float N, float bCN, float adC, float t, float C, float H) {
+ return (C * (bCN * (N - C - H) - adC));
+}
+
+void RK4(float N, float bHN, float bCN, float adH, float adC,
+ float t0, float H0, float C0, float t, float h) {
+ int n = (int)((t - t0) / h);
+ float k1, k2, k3, k4, k5, l1, l2, l3, l4;
+ float H = H0;
+ float C = C0;
+
+ cout << "t\tH\tC\n" << endl;
+ for (int i = 1; i <= n; i++) {
+
+ k1 = h * dhdt(N, bHN, adH, t0, C, H);
+ l1 = h * dcdt(N, bCN, adC, t0, C, H);
+
+ k2 = h * dhdt(N, bHN, adH, t0 + 0.5 * h, C + 0.5 * l1, H + 0.5 * k1);
+ l2 = h * dcdt(N, bCN, adC, t0 + 0.5 * h, C + 0.5 * l1, H + 0.5 * k1);
+
+ k3 = h * dhdt(N, bHN, adH, t0 + 0.5 * h, C + 0.5 * l2, H + 0.5 * k2);
+ l3 = h * dcdt(N, bCN, adC, t0 + 0.5 * h, C + 0.5 * l2, H + 0.5 * k2);
+
+ k4 = h * dhdt(N, bHN, adH, t0 + h, C + l3, H + k3);
+ l4 = h * dcdt(N, bCN, adC, t0 + h, C + l3, H + k3);
+
+ H = H + (1.0 / 6.0) * (k1 + 2 * k2 + 2 * k3 + k4);
+ C = C + (1.0 / 6.0) * (l1 + 2 * l2 + 2 * l3 + l4);
+
+ t0 = t0 + h;
+
+ cout << t0 << "\t" << H << "\t" << C << endl;
+
+ fstream plik1;
+ fstream plik2;
+
+ plik1.open("dane_HA.txt", ios::out|ios::app);
+ plik1 << t0 << "\t" << H << endl;
+ plik1.close();
+
+ plik2.open("dane_CA.txt", ios::out | ios::app);
+ plik2 << t0 << "\t" << C << endl;
+ plik2.close();
+ }
+}
+
+int main() {
+ float t0 = 0, t = 170, H0 = 90, C0 = 90, h = 0.5; //beth israel deconess- c=1, h=1
+ float beta_h, N, alfa_h, delta_h, beta_c, alfa_c, delta_c;
+
+ cout << "Wielkosc populacji: ";
+ cin >> N;
+ cout << "\n---SZCZEP HA-MRSA---" << endl;
+ cout << "Czestotliwosc rozprzestrzeniania sie: ";
+ cin >> beta_h;
+ cout << "Czestotliwosc opuszczania szpitala przez zarazonych pacjentow: ";
+ cin >> delta_h;
+ cout << "Czestotliwosc skutecznego przechodzenia leczenia przez zarazonych pacjentow: ";
+ cin >> alfa_h;
+
+ float bHN = beta_h / N;
+ float adH = alfa_h + delta_h;
+ cout << "\n---SZCZEP CA-MRSA---" << endl;
+ cout << "Czestotliwosc rozprzestrzeniania sie: ";
+ cin >> beta_c;
+ cout << "Czestotliwosc opuszczania szpitala przez zarazonych pacjentow: ";
+ cin >> delta_c;
+ cout << "Czestotliwosc skutecznego przechodzenia leczenia przez zarazonych pacjentow: ";
+ cin >> alfa_c;
+
+ float bCN = beta_c / N;
+ float adC = alfa_c + delta_c;
+ RK4(N, bHN, bCN, adH, adC, t0, H0, C0, t, h);
+}
\ No newline at end of file
--
GitLab