%define constant parameters k1on = 20e-4; k1off = 10e-4; k2on = 10e-4; k2off = 5e-4; k3off = 1e-3; %define initial conditions L0 = 10; %initial concentration of L R0 = 5; %initial concentration of R Rstar0 = 0; %initial concentration of R* C0 = 0; %initial concentration of C