% bms925, process modeling 6

global Q
global Ca
global Ve
global Vi
global B
global Rmax
global Kp5
global tp
global frac_on

Q = 1.0;
%Q = 0.01;
Ca = 1;
Ve = 0.1;
Vi = 1;
B = 1;
Rmax = 1;
Kp5 = 0.4;
tp = 1;
frac_on = 0.20;


tstart = 0;
tend = 10;
%tend = 100;

Ne_init = 0;
Ni_init = 0;

v_init(1) = Ne_init;
v_init(2) = Ni_init;

[t,v] = ode45(@simp_cell2,[tstart tend],v_init);

Ne = v(:,1);
Ni = v(:,2);

plot(t,Ne,t,Ni)