Rename HIV_model_2.m to Legacy Code/HIV_model_2.m

HIV_model_2.m → Legacy Code/HIV_model_2.m

@@ -34,8 +34,8 @@
 %   -> NOTE: parameter pairs are for E and K, respectively
 
 %     P(01) = V i.c. [log(/mL)]
-%     P(02) = E i.c. [/�L]
-%     P(03) = K i.c. [/�L]
+%     P(02) = E i.c. [/無]
+%     P(03) = K i.c. [/無]
 
 %     P(04) = N i.c. [pmol/kg]
 %     P(05) = C absorption rate const. [/d]
@@ -44,10 +44,10 @@
 
 %     P(08) = death rate constant [/d]
 %     P(09)
-%     P(10) = maximum concentration [/�L]
+%     P(10) = maximum concentration [/無]
 %     P(11)
 
-%     P(12) = killing rate constant [�L/d]
+%     P(12) = killing rate constant [無/d]
 %     P(13)
 %     P(14) = Ce killing modifier []
 %     P(15)
@@ -56,7 +56,7 @@
 %     P(17)
 %     P(18) = constant IL15 factor [pM]
 %     P(19)
-%     P(20) = V-dep IL15 factor [pM*�L]
+%     P(20) = V-dep IL15 factor [pM*無]
 %     P(21)
 
 %     P(22) = Ct max level []
@@ -71,13 +71,13 @@
 
 %     P(30) = escape mutation rate []
 %     P(31) = initial escape proportion []
-%     P(32) = response activation threshold [/�L]
+%     P(32) = response activation threshold [/無]
 %     P(33) = response prolif. rate const. [/d]
 
-%     P(34) = T i.c. [/�L]
+%     P(34) = T i.c. [/無]
 %     P(35) = T death rate constant [/d]
-%     P(36) = T maximum concentration [/�L]
-%     P(37) = T infection rate constant [/�L-d]
+%     P(36) = T maximum concentration [/無]
+%     P(37) = T infection rate constant [/無-d]
 %     P(38) = infected T decay rate constant [*(P(35)]
 
 % EquatedPars = vector of indices in 'P' used for equating parameter values
@@ -91,9 +91,9 @@
 
 % Y_MAIN = model output at each time point in 'SolTimes'
 % Y_MAIN(:,1) = V [log(#/mL)]
-% Y_MAIN(:,2) = E [/�L]
-% Y_MAIN(:,3) = K [/�L]
-% Y_MAIN(:,4) = T [/�L]
+% Y_MAIN(:,2) = E [/無]
+% Y_MAIN(:,3) = K [/無]
+% Y_MAIN(:,4) = T [/無]
 
 % Y_EXTRA = cell array of extra outputs (see 'output' section)
 
@@ -131,7 +131,7 @@
 % -------------------------------------------------------------------------
 % Declare initial conditions & calculated parameters ----------------------
 
-V = 10^(P(01)-3) ;% converting V i.c. to [/�L]
+V = 10^(P(01)-3) ;% converting V i.c. to [/無]
 
 V        = V       *10^(3*mutanum) ;% scaling up V to avoid rounding errors
 P(20:21) = P(20:21)/10^(3*mutanum) ;% scaling parameters to match
@@ -280,7 +280,7 @@
 Y = abs(Y)      ;% removing negatives resulting from numerical error
 
 Y_MAIN(:,1) = log10(sum(Y(:,10:10+mutanum),2))+3 ;% V [log(#/mL)]
-Y_MAIN(:,2:4) = Y(:,7:9)                         ;% E:K:T [/�L]
+Y_MAIN(:,2:4) = Y(:,7:9)                         ;% E:K:T [/無]
 
 % -------------------------------------------------------------------------
 % Prepare extra outputs 'Y_EXTRA' -----------------------------------------
@@ -385,4 +385,4 @@
     Y_EXTRA{15}{3,1}(n) =  funOn(2) ;
 end
 
-end
+end