Rename ICsolvers/N803_single.m to N803_single.m

ElsjePienaarGroup · Dec 25, 2022 · b269bba · b269bba
1 parent 682b4a0
commit b269bba
Showing 1 changed file with 19 additions and 19 deletions.
diff --git a/ICsolvers/N803_single.m → N803_single.m b/ICsolvers/N803_single.m → N803_single.m
@@ -8,12 +8,12 @@
 %     |               Pienaar Computational Systems Pharmacology Lab |
 %     \--------------------------------------------------------------/
 %
-% Nomenclature: V  = SIV virions [#/無]
-%               T8 = total CD8+ T cells [#/無]
-%               S0 = resting SIV-specific CD8+ T cells [#/無]
-%               Sa = active  SIV-specific CD8+ T cells [#/無]
-%               N0 = resting non-SIV-specific CD8+ T cells [#/無]
-%               Na = active  non-SIV-specific CD8+ T cells [#/無]
+% Nomenclature: V  = SIV virions [#/無]
+%               T8 = total CD8+ T cells [#/無]
+%               S0 = resting SIV-specific CD8+ T cells [#/無]
+%               Sa = active  SIV-specific CD8+ T cells [#/無]
+%               N0 = resting non-SIV-specific CD8+ T cells [#/無]
+%               Na = active  non-SIV-specific CD8+ T cells [#/無]
 %               X  = N803 at absorption site [pmol/kg]
 %               C  = N803 plasma concentration [pM]
 %               R  = regulation [] (dimensionless quantity)
@@ -59,7 +59,7 @@
 
 % Rename inputed parameters -----------------------------------------------
 Vi   = AllPars(01) ;% V initial value [log(#/mL)]
-SNi  = AllPars(02) ;% S+N initial value [#/無]
+SNi  = AllPars(02) ;% S+N initial value [#/無]
 fS   = AllPars(03) ;% initial frequency: S/(S+N)
 fSA  = AllPars(04) ;% initial frequency: Sa/S
 q    = AllPars(05) ;% V growth rate (if S+N were absent) [/d]
@@ -68,7 +68,7 @@
 V50N = AllPars(08) ;% 50% viral stimulation saturation for N [#/mL]
 mSn  = AllPars(09) ;% normalized Sa reversion rate constant []
 mNn  = AllPars(10) ;% normalized Na reversion rate constant []
-SN50 = AllPars(11) ;% 50% S+N proliferation saturation [#/無]
+SN50 = AllPars(11) ;% 50% S+N proliferation saturation [#/無]
 pS   = AllPars(12) ;% Sa proliferation rate constant [/d]
 pN   = AllPars(13) ;% Na proliferation rate constant [/d]
 d    = AllPars(14) ;% S0/N0 death rate constant [/d]
@@ -92,9 +92,9 @@
 %% ------------------------------------------------------------------------
 % Calculate some initial conditions & parameters --------------------------
 
-Vi   = 10^(Vi-3) ;% V initial value [#/無]
-V50S = V50S/1000 ;% 50% viral stimulation saturation for S [#/無]
-V50N = V50N/1000 ;% 50% viral stimulation saturation for N [#/無]
+Vi   = 10^(Vi-3) ;% V initial value [#/無]
+V50S = V50S/1000 ;% 50% viral stimulation saturation for S [#/無]
+V50N = V50N/1000 ;% 50% viral stimulation saturation for N [#/無]
 
 % restrict mS and mN such that initial activation aS and aN are positive
 US = 2*(2*pS/(pS+dA))^7 ;
@@ -146,8 +146,8 @@
 p0 = pi * (SN50+SNi)/SN50 * (1+p2)  ;% S0/N0 proliferation rate con [/d]
 p1 = pm/p0  ;% S0/N0 proliferation stimulation factor
 
-gS0 = q / ( SA/(1+gS2) + psi*NA/(1+gN2) )  ;% Sa killing rate con [無/#-d]
-gN0 = psi*gS0  ;% Na killing rate constant [無/#-d]
+gS0 = q / ( SA/(1+gS2) + psi*NA/(1+gN2) )  ;% Sa killing rate con [無/#-d]
+gN0 = psi*gS0  ;% Na killing rate constant [無/#-d]
 
 sS  = dR / ( Vi/(V50S+Vi) + sig*Vi/(V50N+Vi) )  ;% regulation gen by S act
 sN = sig*sS  ;% regulation generation by N activation
@@ -156,17 +156,17 @@
 % Prepare parameter and initial value vectors and call 'N803_model_2' -----
 
 Pars(01) = q    ;% V growth rate (if S+N were absent) [/d]
-Pars(02) = gS0  ;% Sa killing rate constant [無/#-d]
-Pars(03) = gN0  ;% Na killing rate constant [無/#-d]
+Pars(02) = gS0  ;% Sa killing rate constant [無/#-d]
+Pars(03) = gN0  ;% Na killing rate constant [無/#-d]
 
-Pars(04) = V50S ;% 50% viral stimulation saturation for S [#/無]
-Pars(05) = V50N ;% 50% viral stimulation saturation for N [#/無]
+Pars(04) = V50S ;% 50% viral stimulation saturation for S [#/無]
+Pars(05) = V50N ;% 50% viral stimulation saturation for N [#/無]
 Pars(06) = aS0  ;% S0 activation rate constant [/d]
 Pars(07) = aN0  ;% N0 activation rate constant [/d]
 Pars(08) = mS   ;% Sa reversion rate constant [/d]
 Pars(09) = mN   ;% Na reversion rate constant [/d]
 
-Pars(10) = SN50 ;% 50% S+N proliferation saturation [#/無]
+Pars(10) = SN50 ;% 50% S+N proliferation saturation [#/無]
 Pars(11) = p0   ;% S0/N0 proliferation rate constant [/d]
 Pars(12) = pS   ;% Sa proliferation rate constant [/d]
 Pars(13) = pN   ;% Na proliferation rate constant [/d]
@@ -214,4 +214,4 @@
 
 Pars = [ Pars Yic ]  ;% adding initial conditions to parameter output
 
-end
+end