begin parameters kp1 0.00001 km1 0.01 kp2 0.00001 km2 0.01 kp3 0.00001 km3 0.01 R0 100 L0 100 end parameters begin molecule types R(r,r) L(l,l) end molecule types begin seed species R(r,r) R0 L(l,l) L0 end seed species begin reaction rules R(r) + L(l,l) <-> R(r!1).L(l!1,l) kp1,km1 # Ligand addition R(r) + L(l,l!2) <-> R(r!1).L(l!1,l!2) kp2,km2 # Chain elongation # R(r).L(l) <-> R(r!1).L(l!1) kp3,km3 # Ring closure end reaction rules generate_network({max_stoich=>{R=>5,L=>5}});