An object-oriented MATLAB implementation of the Gillespie [1] and the 'Brusselator' [2] stochastic simulation algorithms.
This project requires only MATLAB to run the stochastic simulation algorithms, however the Signal Processing Toolbox and the Statistics and Machine Learning Toolbox are used for subsequent analysis of the simulations.
Run the main.m file with the class file ChemicalReaction.m in the same directory, and the ChemicalReaction objects will be instantiated, calculations will be executed, and results of the simulation and analysis will be displayed.
[1] Gillespie, D. T. (1977). “Exact Stochastic Simulation of Coupled Chemical Reactions”. The Journal of Physical Chemistry. 81 (25): 2340-2361. doi: 10.1021/j100540a008
[2] I. Prigogine, From Being to Becoming: Time and Complexity in the Physical Sciences, New York: W.H. Freeman and Company, 1980