%% Pulse Generator Network model. 
% ODE of the Pulse Generator Network. Parameters in input are:
% - t is the time
% - x is the vector of state variables
% - p is the parameter vector
% - u is the vector of inputs
% - xT is the vector of total concentrations

function dx=PulseGeneratorNetwork(t,x,p,u,xT)

%names of model parameters
k1=p(1); 
K1=p(2); 
lambda2=p(3); 
k12=p(4); 
K2=p(5);  
lambda=p(6);

%input signal
S1=u(1);

%state variables
R2=x(1);
Y=x(2);

%ODEs
dR2 = k1*((S1/K1)^3)/(1+(S1/K1)^3) - lambda2*R2;
dY = k12/(1+(R2/K2)^3)*((S1/K1)^3/(1+(S1/K1)^3)) - lambda*Y;

dx=[dR2;dY];
    
end