function dy = cbc(t,y) 
% dy = cbc(t,y)

% Kåre A. Madsen, Nov. 30, 1999.

% The parameters.
mum   = 0.9;
s0    = 0.005;
K     = 0.15;
gamma = 0.05;
D     = 0.3;
alfa  = K + gamma+D;
YP    = 0.5;
SR    = 0.25; 

% The model.
dy = [(mum*y(3)/(y(3) + s0))*y(1) - alfa*y(1);
      gamma*y(1) - D*y(2);
      D*SR - (mum*y(3)/(y(3) + s0))*y(1)/YP - D*y(3)];

--------------------------------------------------------------------------

function dy = rdbd1(t,y) 
% dy = rdbd1(t,y) 

% Kåre A. Madsen, Nov. 30, 1999.

% The parameters.
alfa     = 0.01;
epsilonr = 0;
epsilonf = 0; 

% The model.
dy = [2*y(1) - alfa*y(1)*y(2) - epsilonr;
       -y(2) + alfa*y(1)*y(2) - epsilonf];

--------------------------------------------------------------------------

function dy = rdbd2(t,y) 
% dy = rdbd2(t,y) 

% Kåre A. Madsen, Nov. 30, 1999.

% The parameters.
alfa     = 0.01;
R        = 500;
epsilonr = 0;
epsilonf = 0; 

% The model.
dy = [2*(1-y(1)/R)*y(1) - alfa*y(1)*y(2) - epsilonr;
      -y(2)+alfa*y(1)*y(2) - epsilonf];

--------------------------------------------------------------------------

function dy = rdbd3(t,y) 
% dy = rdbd3(t,y) 

% Kåre A. Madsen, Nov. 30, 1999.

% The parameters.
alfa     = 0.01;
R        = 500;
F        = 500;
epsilonr = 0;
epsilonf = 0; 

% The model.
dy = [2*(1-y(1)/R)*y(1) - alfa*y(1)*y(2) - epsilonr;
      -(1-y(2)/F)*y(2)+alfa*y(1)*y(2) - epsilonf];

--------------------------------------------------------------------------

function dy = tbl(t,y)
% dy = tbl(t,y)

% Kåre A. Madsen, Nov. 30, 1999.
 
% The model (dy/dt = A*y).
dy = [-0.3  0.1  0.3; 0.2 -0.2  0.1; 0.1  0.1 -0.4]*y;