% This code is for implementing RT 
% Radiation therapy first and then Chemotherapy 

close all;
clear all;
clc;

% This code is for Radiation Therapy

t_final = 40;
%D = 0:1:45;
D = 40;
t = 1:0.01:t_final;

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

for i = 1:length(D)
    i
    x(:,2*i-1:2*i) = ode4(@sys_ndRad,t,[1;0],D(i)); 
    
    PMCells_after_Saturation(i) = x(end,2*i)
    ERR(i) = x(end,2*i)*1.2;  
end

%ERR_lung

Ncells_EndofRad = x(end,1*i);
PMcells_EndofRad = x(end,2*i);

figure(1)
%hold on
plot(D,ERR)

% %------- Thyroid ERR data due to HL from Shuryak--------% 
% 
% A = dlmread('LungDataHL.txt','\t');
% 
% x = A(:,1);
% y = A(:,2);
% 
% M = dlmread('LungErrorBars.txt','\t');
% 
% temp1 = M(:,1);
% temp2 = M(:,2);
% 
% for j = 1:length(temp2)    
%     if (mod(j,2) == 0)        
%         ly(j) = temp2(j);
%     end
%     if (mod(j,2) == 1) 
%         uy(j) = temp2(j);
%     end
% end
% 
% l = ly(find(ly~=0));
% u = uy(find(uy~=0));
% for i = 1:length(y)    
%     L(i) = y(i) - l(i);
%     U(i) = -y(i) + u(i);
% end
% errorbar(x,y,U,L)
% %scatter(x,y,'bo');
% hold off
% 
%