% Chemotherapy model

global I k u b a km kn

I = 0.005;	% infusion rate

k = 0.002;	% uptake rate of drug by cells

u = 0.005; 	% Removal rate of drug by blood

b = 0.1;		% 1/maximum tumor size

a = 0.01;		% Gompertz growth constant

km = 0.6;		% Michaelis-Menton constant

kn = 0.4;		% Michaelis-Menton constant

[t,y] = ode45('chemo_func',[0 200],[1; 0]);

plot(t,y(:,1),'-',t,y(:,2),'-.')

xlabel('Time')

title('Chemotherapy model')

text(100,1.2,'I=0.005, k=0.002, u=0.005')

text(100,1.0,'a=0.01,b=0.1,km=0.6, kn=0.4')

axis([0 200 0 1.3])