--- a
+++ b/Sequential/sys_ndChemo.m
@@ -0,0 +1,23 @@
+
+function xprime = sys_ndChemo(t,x,D,OneCycle,TotalDays)
+
+xprime = zeros(2,1);
+
+N = 10^9; %
+lambda = 0.40; % Sachs, radiation
+r = 0.8175; % Sachs, radiation
+
+gamma = 10^(-10); % fitting parameter
+
+alphad = 0.1333; 
+alphaC = 0.;
+
+if(t<=TotalDays)    
+        chemo = exp(-alphad*mod(t,OneCycle));
+        xprime(1) = lambda*x(1)*(1-x(1)) - alphaC*chemo*(D/OneCycle)*x(1) - gamma*(D/OneCycle)*x(1);
+        xprime(2) = r*lambda*x(2)*(1-x(1)) - alphaC*chemo*(D/OneCycle)*x(2) + gamma*(D/OneCycle)*N*x(1);
+else
+    xprime(1) = lambda*x(1)*(1-x(1));
+    xprime(2) = r*lambda*x(2)*(1-x(1));
+end
+