% Chemotherapy ODE functions

function out = chemo_func(t,y)

global I k u b a km kn

out = [-a*y(1)*log(b*y(1)) - km*y(1)*y(2)/(kn+y(2)); I - k*y(1)-u*y(2)];