% bms925, process modeling 5

global Q
global Cmfa
global V
global kon
global koff
global Ptotal


Q = 1.0;
Cmfa = 1;
V = 1;
kon = 10;
%kon = 1000;
koff = 10;
Ptotal = 1;


tstart = 0;
tend = 10;

Ncx_init = 0;
Nmf_init = 0;

v_init(1) = Ncx_init;
v_init(2) = Nmf_init;

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

Ncx = v(:,1);
Nmf = v(:,2);

Npunoc = Ptotal - Ncx;
Nmtotal = Ncx + Nmf;

plot(t,Ncx,t,Nmf,t,Npunoc,t,Nmtotal)

legend('Ncx','Nmf','Npunouc','Nmtotal');