How to use Gillespie's algorithm in MATLAB to display chemical reactions

Hi everyone, I am trying to use MATLAB to do Gillespie's algorithm for chemical reactions:

A + B --> C with rate f,

and C --> decay with rate d,

in which f = d = 1/s.

At time 0, A(0) = 20, B(0) = 10 and C(0) =10.

Could you please help me with this, probably using random number generator to represent probalities of each reactions?

Thank you all very much!

Answer 1

[Dinesh Yadav]

Hi Yuhong

Kindly go through the following links for existing implementations.

https://www.mathworks.com/matlabcentral/answers/224604-implementing-gillespie-s-algorithm

https://www.mathworks.com/matlabcentral/fileexchange/34707-gillespie-stochastic-simulation-algorithm

https://www.mathworks.com/matlabcentral/fileexchange/34151-gillespie-s-direct-method-stochastic-simulation-algorithm

Hope it helps.