Hi,
You correctly understand the current constraints of SimBiology's stochastic solver: reactions must use linear, mass action kinetics. I can't think of any easy workaround for this constraint. If I needed to do it, I would also proabably implement my own solver and profile it to figure out where the performance bottlenecks are.
The idea of supporting other kinds of kinetics has come up occasionally before. If I remember correctly, one reason we have not done this is because we were concerned about correctly generalizing our mathematical approach. Do you have references that discuss the SSA/Gillespie algorithm for other kinetic laws?
-Arthur
Best Answer