--- a
+++ b/Sequential/duringrad.m
@@ -0,0 +1,65 @@
+% 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
+% 
+%