diff --git a/HIV_model_4.m b/Legacy Code/HIV_model_4.m
similarity index 95%
rename from HIV_model_4.m
rename to Legacy Code/HIV_model_4.m
index 5d95d84..744574d 100644
--- a/HIV_model_4.m
+++ b/Legacy Code/HIV_model_4.m	
@@ -42,9 +42,9 @@
 % normout     = scalar used to change units of Y_MAIN (see OUTPUTS)
 %               (if empty, normout is ignored)
 
-%     P(01) = V i.c. [log(/mL)] (will convert to [/µL])
-%     P(02) = E i.c. [/µL]
-%     P(03) = K i.c. [/µL]
+%     P(01) = V i.c. [log(/mL)] (will convert to [/ç„¡])
+%     P(02) = E i.c. [/ç„¡]
+%     P(03) = K i.c. [/ç„¡]
 
 %     P(04) = Q i.c. [pmol/kg]
 %     P(05) = C absorption rate constant [/d]
@@ -53,11 +53,11 @@
 
 %     P(08) = death rate constant [/d]
 %     P(09)
-%     P(10) = max proliferating concentration [/µL]
+%     P(10) = max proliferating concentration [/ç„¡]
 %     P(11)
 
-%     P(12) = total initial killing rate [/d]     (will convert to [/µL-d])
-%     P(13) =                            [*P(12)] (will convert to [/µL-d])
+%     P(12) = total initial killing rate [/d]     (will convert to [/ç„¡-d])
+%     P(13) =                            [*P(12)] (will convert to [/ç„¡-d])
 %     P(14) = drug killing modifier []
 %     P(15)
 
@@ -82,10 +82,10 @@
 %     P(32) = REG killing modifier []
 %     P(33)
 
-%     P(34) = T i.c. [/µL]
+%     P(34) = T i.c. [/ç„¡]
 %     P(35) = T death rate constant [/d]
-%     P(36) = T max proliferating concentration [/µL]
-%     P(37) = T infection by V1 rate constant [/µL-d]
+%     P(36) = T max proliferating concentration [/ç„¡]
+%     P(37) = T infection by V1 rate constant [/ç„¡-d]
 %     P(38) = infected T death rate [*P(35)] (will convert to [/d])
 
 %     P(39) = V2 initial frequency []
@@ -99,9 +99,9 @@
 
 %                                                         if normout = 1
 % Y_MAIN(:,1) = V at points in 'SolTimes' [log(#/mL)]     [log fold change]
-% Y_MAIN(:,2) = E at points in 'SolTimes' [/µL]           [fold change]
-% Y_MAIN(:,3) = K at points in 'SolTimes' [/µL]           [fold change]
-% Y_MAIN(:,4) = T at points in 'SolTimes' [/µL]           [fold change]
+% Y_MAIN(:,2) = E at points in 'SolTimes' [/ç„¡]           [fold change]
+% Y_MAIN(:,3) = K at points in 'SolTimes' [/ç„¡]           [fold change]
+% Y_MAIN(:,4) = T at points in 'SolTimes' [/ç„¡]           [fold change]
 %                                                         (see INPUTS)
 
 % Y_ALL = cell array of extra outputs (see end of script)
@@ -136,7 +136,7 @@
 % -------------------------------------------------------------------------
 % Declare initial conditions & calculated parameters ----------------------
 
-P(01) = 10^(P(01)-3) ;% convert from [log(/mL)] to [/µL]
+P(01) = 10^(P(01)-3) ;% convert from [log(/mL)] to [/ç„¡]
 
 Yic       = zeros(1,In(5))                          ;% initial conditions
 Yic(7:11) = [P([02,03,34])' P(01)*[1-P(39) P(39)] ] ;% E,K,T,V1,V2 i.c.
@@ -147,7 +147,7 @@
 if P(41) == 0 ; P(41) = 1 ; end % if V2 targeting factor is zero (off)...
 if P(42) == 0 ; P(42) = 1 ; end % then set to 1 (no reduction in targeting)
 
-% calculate beta [µL/d] and b [µL/d]
+% calculate beta [ç„¡/d] and b [ç„¡/d]
 if P(39) == 0 % for single strain
     beta = [ (sum(P(12:13)) + P(38)) / P(34) ; 0 ] ;
     b    = [ P(37)                           ; 0 ] ;
@@ -161,7 +161,7 @@
 end
 
 P(30:31) = P(30:31).*P(18:19) ;% convert from [*P(18:19)] to [/d]
-P(12:13) = P(12:13)./P(02:03) ;% convert from [/d] to [/µL-d]
+P(12:13) = P(12:13)./P(02:03) ;% convert from [/d] to [/ç„¡-d]
 
 % if a maximum proliferating concentration is zero, set it as very large
 % (this avoids division by zero and approximately removes its effect)
@@ -296,7 +296,7 @@
 if normout == 1   % units are [log fold change] and [fold change]
     Y_MAIN(:,1)   = log10( (Y(:,10)+Y(:,11)) / P(01)        ) ;% V
     Y_MAIN(:,2:4) =         Y(:,7:9)        ./ P([02,03,34])' ;% E,K,T
-else              % units are [log(#/mL)] and [/µL]
+else              % units are [log(#/mL)] and [/ç„¡]
     Y_MAIN(:,1)   = log10( Y(:,10)+Y(:,11) ) + 3 ;% V
     Y_MAIN(:,2:4) =        Y(:,7:9)              ;% E,K,T
 end
@@ -415,4 +415,4 @@
     Y_ALL{15}{2,1}(n) = Y(n,6)  ;
 end
 
-end
\ No newline at end of file
+end