|
a |
|
b/ChemoOnly/sys_ndChemo.m |
|
|
1 |
|
|
|
2 |
function xprime = sys_ndChemo(t,x,D,OneCycle,TotalDays) |
|
|
3 |
|
|
|
4 |
xprime = zeros(2,1); |
|
|
5 |
|
|
|
6 |
N = 10^9; |
|
|
7 |
lambda = 0.40; |
|
|
8 |
r = 0.8175; |
|
|
9 |
|
|
|
10 |
gamma = 10^(-10); |
|
|
11 |
|
|
|
12 |
alphad = 0.1333; |
|
|
13 |
alphaC = 0.2; |
|
|
14 |
|
|
|
15 |
if(t<=TotalDays) |
|
|
16 |
|
|
|
17 |
chemo = exp(-alphad*mod(t,OneCycle)); |
|
|
18 |
|
|
|
19 |
xprime(1) = lambda*x(1)*(1-x(1)) - alphaC*chemo*(D/OneCycle)*x(1) - gamma*(D/OneCycle)*x(1); |
|
|
20 |
xprime(2) = r*lambda*x(2)*(1-x(1)) - alphaC*chemo*(D/OneCycle)*x(2) + gamma*(D/OneCycle)*N*x(1); |
|
|
21 |
else |
|
|
22 |
|
|
|
23 |
xprime(1) = lambda*x(1)*(1-x(1)); |
|
|
24 |
xprime(2) = r*lambda*x(2)*(1-x(1)); |
|
|
25 |
|
|
|
26 |
end |
|
|
27 |
|
|
|
28 |
|
|
|
29 |
|
|
|
30 |
|
|
|
31 |
|
|
|
32 |
|