%ODE example equations
function dYdt = ODEexample_equations(t,Y,params)

%this function defines the ODEs describing the bimolecular
%reaction L + R <--> C with association and dissociation
%rates kon and koff

%initial conditions
L = Y(1);
R = Y(2);
C = Y(3);

%parameters
kon = params(1);
koff = params(2);

dLdt = -kon*L*R + koff*C;
dRdt = -kon*L*R + koff*C;
dCdt = kon*L*R - koff*C;

dYdt = [dLdt; dRdt; dCdt];