Hi Sebastian,
I'm not an expert on the Gillespie algorithm, but from what I understand, this approach should be compatible.
The main issue with using constant rates is the risk of numerical instability during integration. You need to ensure that there's always a sufficient amount of substrate A available to support the reaction rate over a given time step (Δt); otherwise, you may end up with negative concentrations.
Moreover, systems without saturation mechanisms can easily become unstable, diverging towards ±∞ under certain conditions. To avoid this, I often model reactions with saturation kinetics using very low values. This effectively yields zero-order kinetics—since the system is nearly always saturated—while maintaining stability even when A is depleted.
Best regards,
Matthias
--
You received this message because you are subscribed to the Google Groups "COPASI User Forum" group.
To unsubscribe from this group and stop receiving emails from it, send an email to copasi-user-fo...@googlegroups.com.
To view this discussion visit https://groups.google.com/d/msgid/copasi-user-forum/8ee131bb-26ee-45c8-a6a3-2eeeb779f456n%40googlegroups.com.