% small metab network
% S -> P -> A -> degrade
% A inhibite S -> P

global Vmax_msp
global Ca
global Kmsp
global Vmaxpa
global Kmap
global B
global V

global Ns   % initial S
global Np   % initial P
global Na   % initial A


Vmax_msp = 1;
Ca = 0.2;
Kmsp = 0.2;
Vmaxpa = 1;
Kmap = 0.2;
B = 0.05;
V = 1;

Ns = 1;
%Ns = 10;
%Ns = 100;
Np = 0;
Na = 0;

tstart = 0;
tend = 10;
%delta_t = (tend-tstart)./1000;  %v2

v_init(1) = Np;
v_init(2) = Na;


%Ca_v = [0.1 1.0 100.0];     %v2


%clear Np_t  %v2
%clear Na_t  %v2
%for i = 1:length(Ca_v)  %v2

%    Ca = Ca_v(i);   %v2
    
    [t,v] = ode45(@met_network,[tstart tend],v_init);
%    [t,v] = ode45(@met_network,[tstart:delta_t:tend],v_init);   %v2

    Np = v(:,1);
    Na = v(:,2);

%    Np_t(:,i) = v(:,1);    %v2
%    Na_t(:,i) = v(:,2);    %v2
    
%end    %v2

plot(t,Np,t,Na)     %v2
%plot(t,Np_t,'b-',t,Na_t,'r--')     %v2