diff --git a/HIV_model_5.m b/Legacy Code/HIV_model_5.m
similarity index 96%
rename from HIV_model_5.m
rename to Legacy Code/HIV_model_5.m
index 2c77898..aba9672 100644
--- a/HIV_model_5.m
+++ b/Legacy Code/HIV_model_5.m	
@@ -36,8 +36,8 @@
 
 %     Initial Conditions
 %     P(01) = virion i.c. [log(/mL)]
-%     P(02) = CD8+ T cell (R+A+E) i.c. [/µL]
-%     P(03) = NK cell i.c. [/µL]
+%     P(02) = CD8+ T cell (R+A+E) i.c. [/ç„¡]
+%     P(03) = NK cell i.c. [/ç„¡]
 
 %     Pharmacokinetic Parameters
 %     P(04) = N803 injection site i.c. [pmol/kg]
@@ -58,7 +58,7 @@
 %     Parameters governing R,K proliferation and death
 %     P(16) = R death rate constant [/d]
 %     P(17) = K "
-%     P(18) = R max proliferating concentration [/µL]
+%     P(18) = R max proliferating concentration [/ç„¡]
 %     P(19) = K "
 %     P(20) = R drug-induced proliferation factor [/d]
 %     P(21) = K "
@@ -70,12 +70,12 @@
 %     P(25) = E initial killing rate (a_E * E_0) [/d]
 %     P(26) = K "
 %     P(27) = R drug-induced activation factor [/d]
-%     P(28) = E drug-induced killing factor [µL/d]
+%     P(28) = E drug-induced killing factor [ç„¡/d]
 %     P(29) = K "
 %     P(30) = R activation suppression factor []
 %     P(31) = E killing suppression factor []
 %     P(32) = K "
-%     P(33) = R half-activation virus concentration [/µL]
+%     P(33) = R half-activation virus concentration [/ç„¡]
 
 %     Parameters governing CD8+ T cell clonal expansion and retraction
 %     P(34) = A death rate constant [/d]
@@ -85,10 +85,10 @@
 %     P(38) = E deactivation rate constant [/d]
 
 %     Parameters governing CD4+ T cells
-%     P(39) = CD4+ T cell i.c. [/µL]
+%     P(39) = CD4+ T cell i.c. [/ç„¡]
 %     P(40) = T death rate constant [/d]
-%     P(41) = T max proliferating concentration [/µL]
-%     P(42) = T infection rate constant [/µL-d]
+%     P(41) = T max proliferating concentration [/ç„¡]
+%     P(42) = T infection rate constant [/ç„¡-d]
 %     P(43) = infected T decay rate constant [*(P(40)]
 
 %% ========================================================================
@@ -97,9 +97,9 @@
 
 % Y_MAIN = model output at each time point in 'SolTimes'
 % Y_MAIN(:,1) = virion [log(#/mL)]
-% Y_MAIN(:,2) = CD8+ T cell (R+A+E) [/µL]
-% Y_MAIN(:,3) = NK cell [/µL]
-% Y_MAIN(:,4) = CD4+ T cell [/µL]
+% Y_MAIN(:,2) = CD8+ T cell (R+A+E) [/ç„¡]
+% Y_MAIN(:,3) = NK cell [/ç„¡]
+% Y_MAIN(:,4) = CD4+ T cell [/ç„¡]
 
 % Y_ALL = cell array of extra outputs (see 'output' section)
 
@@ -133,7 +133,7 @@
 % -------------------------------------------------------------------------
 % Declare initial conditions & calculated parameters ----------------------
 
-V0 = 10^(P(01)-3) ;% convert V i.c. to [/µL]
+V0 = 10^(P(01)-3) ;% convert V i.c. to [/ç„¡]
 
 % calculate i.c. for CD8+ T cell subsets (R,A,E)
 CD8TCOEFS          = [ ones(1,P(36)+2) ; zeros(P(36)+1,P(36)+2) ] ;
@@ -156,7 +156,7 @@
 % redefine P(25:26) as E,K per-capita killing rate constants
 P(25:26) = P(25:26)./ Yic(6:7) ;
 
-% (infection rate constant) * (infected cell to virus conversion) [/µL-d]
+% (infection rate constant) * (infected cell to virus conversion) [/ç„¡-d]
 beta = ( P(40)*P(43) + sum(P(25:26).*Yic(6:7)) ) / Yic(8) ;
 
 % calculate [R;K;T] proliferation rate [/d]
@@ -305,9 +305,9 @@
 Y = abs(Y)	;% removing negatives resulting from numerical error
 
 Y_MAIN(:,1) = log10(Y(:,9))+3               ;% virion [log(#/mL)]
-Y_MAIN(:,2) = sum(Y(:,[5,6,10:in(1)]),2)  	;% CD8+ T cell (R+A+E) [/µL]
-Y_MAIN(:,3) = Y(:,7)                        ;% NK cell [/µL]
-Y_MAIN(:,4) = Y(:,8)                        ;% CD4+ T cell [/µL]
+Y_MAIN(:,2) = sum(Y(:,[5,6,10:in(1)]),2)  	;% CD8+ T cell (R+A+E) [/ç„¡]
+Y_MAIN(:,3) = Y(:,7)                        ;% NK cell [/ç„¡]
+Y_MAIN(:,4) = Y(:,8)                        ;% CD4+ T cell [/ç„¡]
 
 % -------------------------------------------------------------------------
 % Prepare extra outputs 'Y_ALL' -----------------------------------------
@@ -467,4 +467,4 @@
     Y_ALL{20}{3,1}(n) = Yn(4) ;
 end
 
-end
\ No newline at end of file
+end