Rate law dependent synthesis

Enable the following model:

```python
with reaction_rules():
    ~X > X | (k * X / (K + X))

run_simulation(t, solver='gillespie')
```

Except for `ode`, no solver supports rate law dependent reactions. This proposal suggests accepting the rate law for synthetic reactions (zero-th order).