clear all

global cH cT cV cE cVT aET sH sT sV sE sVT sET dH dT dV dE ... 
       mH mT mV mE D P0 r dsgC 
   
cH = 4; 
cT = 6; 
cV = 0.05; 
cE = 2;
cVT = 4; 
aET = .25;

sH = .05; 
sT = .002; 
sV = .12; 
sE = .0001;
sVT = 1; 
sET = .0012;

dH = .008;
dT = .48;  
dV = 0.01; 
dE = 0.1; 

mH = .1;
mT = 2; 
mV = .1; 
mE=.1;

cT = 5.85;

dE = .10108; %periodic

dE = .1010775;



D = 1008;
P0 = .00035;

r = .7;

init=[0.2,0.001,0,.00001,.0003];
%init=[0.0363    0.1523    0.0019    0.0222    0.0003];
%init=[0.0526    0.0724    0.0010    0.0565    0.0003];
Tran = 0; %
%Evol = input('Evolution Time:' ); % Length of evolution after transient. 
Evol = 3000; %
blah=odeset('RelTol', 1e-16, 'AbsTol', [1e-8, 1e-16, 1e-8, 1e-16, 1e-8], 'BDF', 'on');

[t,y]=ode15s('tumormodel',[0 Evol],init,blah);

figure(1)
plot(t,y(:,1),'r-',t,y(:,2),'k-',t,y(:,3),'m-',t,y(:,4),'b-',t,y(:,5),'g-')

figure(2) 
plot3(y(:,1),y(:,2),y(:,4))
xlabel('H')
ylabel('T')
zlabel('E')
