begin parameters
  R0   100
  L0   500
  A0   100
  kon  0.01
  koff  0.1
  kAon  0.01
  kAoff  0.1
  kAp 0.01
  kAdp 0.1
end parameters

begin molecule types
   R(l,r,a)
   L(r)
   A(r,Y~U~P)
end molecule types

begin seed species
  R(l,r,a)   R0    
  L(r)   L0
  A(r,Y~U)   A0
end seed species

begin reaction rules
   R(l) + L(r) <-> R(l!1).L(r!1)    kon,koff
   R(a) + A(r) <-> R(a!1).A(r!1)    kAon,kAoff
   R(l!1,r) + R(l!2,r) <-> R(l!1,r!3).R(l!2,r!3)   kon,koff
   R().R().A(Y~U) -> R().R().A(Y~P)  kAp
   A(Y~P) -> A(Y~U)   kAdp
end reaction rules

begin observables
  Molecules R_Dim_M1  R().R()
  Molecules R_Dim_M2  R(r!+)
  Species R_Dim_S R(r!+)
end observables

generate_network();
writeSBML();
simulate_ode({t_end=>50,n_steps=>20});