Rename HIV_model_3.m to Legacy Code/HIV_model_3.m

HIV_model_3.m → Legacy Code/HIV_model_3.m

@@ -33,8 +33,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) = Q i.c. [pmol/kg]
 %     P(05) = C absorption rate constant [/d]
@@ -43,10 +43,10 @@
 
 %     P(08) = death rate constant [/d]
 %     P(09)
-%     P(10) = max proliferating concentration [/�L]
+%     P(10) = max proliferating concentration [/無]
 %     P(11)
 
-%     P(12) = killing rate constant [�L/d]
+%     P(12) = killing rate constant [無/d]
 %     P(13)
 %     P(14) = drug killing modifier
 %     P(15)
@@ -68,13 +68,13 @@
 
 %     P(28) = escape mutation rate []
 %     P(29) = initial escape proportion []
-%     P(30) = response activation threshold [/�L]
+%     P(30) = response activation threshold [/無]
 %     P(31) = response proliferation rate constant [/d]
 
-%     P(32) = T i.c. [/�L]
+%     P(32) = T i.c. [/無]
 %     P(33) = T death rate constant [/d]
-%     P(34) = T max proliferating concentration [/�L]
-%     P(35) = T infection rate constant [/�L-d]
+%     P(34) = T max proliferating concentration [/無]
+%     P(35) = T infection rate constant [/無-d]
 %     P(36) = infected T decay rate constant [*(P(33)]
 
 % Delays      = number of delay chambers for:
@@ -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)
 
@@ -146,7 +146,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(30)    = P(30)   *10^(3*mutanum) ;
@@ -308,7 +308,7 @@
 Y = abs(Y)      ;% removing negatives resulting from numerical error
 
 Y_MAIN(:,1) = log10(sum(Y(:,10:mutend),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' -----------------------------------------
@@ -418,4 +418,4 @@
     Y_EXTRA{15}{3,1}(n) = Y(n,6) ;
 end
 
-end
+end