Introducing the Stochastic Simulation of Chemical Reactions: Using the Gillespie Algorithm and MATLAB
This paper illustrates the numerical (stochastic and deterministic) simulation of the AB equilibrium process as a suitable tool for introducing the stochastic modeling of chemical reactions. The Gillespie algorithm was selected for the stochastic simulation. Students can develop their own MATLAB programs to simulate the process and then use the results to discover the main differences between the stochastic and deterministic simulation. Two free commercial simulation programs are also proposed as an auxiliary tool for extending the stochastic simulation to other reactions.