Update lotkavolterra_antibiotics.m

lotkavolterra_antibiotics.m

@@ -78,8 +78,13 @@ function dxdt = lotka_volterra_with_antibiotics(t, x, r, K, A, antibiotics)
 
     % Növekedési ráták csökkentése az időfüggő antibiotikum hatással
     r_mod = r_mod .* (1 - antibiotic_effect * sensitivity_vector');
-    resistance_factor = 1 + 0.01 * t; % Rezisztencia fokozatos növekedése
+    
+    scaled_time = t - 30; % Átskálázott idő: antibiotikum hatás kezdete legyen a 0
+    resistance_growth_rate = 0.01;
+    resistance_factor = 1 + resistance_growth_rate * scaled_time;
     r_mod = r_mod .* (1 ./ resistance_factor);
+    max_increase = 1.5; % Maximális eltartóképesség növekedési faktor (pl. 1.5 = 50%-os növekedés)
+    K_mod = K_mod .* (1 + (max_increase - 1) * (1 - exp(-resistance_growth_rate * scaled_time))); 
 
     % Az antibiotikumok antagonisztikus hatásának figyelembevétele
         interaction_modifier = ones(size(r_mod)); % Kezdetben nincs hatás