% This code is for implementing the Chemotherapy

% Travis et al. 2002: Table 2A - Relative Risk is 4.2 without the dose

% Table 4A = average of all doses categories is 5.35 and average dose=50

close all;

clear all;

clc;

OneCycle = 30; % number of days per cycle, Travis 2002

NumbCycles = 3; % number of cycles, Historical studies

TotalDays = OneCycle * NumbCycles; 

t_final = TotalDays + 50; % End of treatment

t = 1:0.01:t_final;

D = 20;

x = zeros(length(t),length(D)*2);

for i = 1:length(D)

    i

    x(:,2*i-1:2*i) = ode4(@sys_ndChemo,t,[1;0],D(i),OneCycle,TotalDays); 

    PMCells_after_Saturation(i) = x(end,2*i);

    ERR_Breast(i) = x(end,2*i)*1;  

end

Ncells_EndofChemo = x(end,1*i);

PMcells_EndofChemo = x(end,2*i);

ERR_Breast

RR = 1+ERR_Breast;

%figure(1)

%plot(t,x(:,2))