Die Skizze zeigt den sogenannten Apfelkrieg. Von beiden Seiten werden Äpfel über den Zaun geworfen. Da von links mehr Äpfel nach rechts geworfen werden als umgekehrt, stellt sich ein Gleichgewicht ein, bei dem sich mehr Äpfel rechts befinden. 
Hierbei ist NA0 die Gesamtzahl der Äpfel (Teilchen) und x ist der Umsatz, bzw. die Anzahl der Produkteilchen. Die Geschwindigkeitskonstante für die Hinreaktion ist khin, für die Rückreaktion krueck und dt ist das Zeitintervall, in dem eine bestimmte Anzahl an Äpfeln geworfen wird.

https://doi.org/10.1002/ckon.201900053
https://doi.org/10.1021/acs.jchemed.2c00203
https://doi.org/10.1021/acs.jchemed.0c00081
 

Mit dem TigerJython-Code kann das System simuliert werden. Es handelt sich um die numerische Lösung einer Differentialgleichung. Allerdings werden nicht nur ganze Äpfel, sondern auch Bruchteile von Äpfeln zwischen den beiden Seiten verschoben. Deshalb wird zusätzlich der gerundete Wert für x ausgegeben.

Kommentierter Code für TigerJython

Alternativ kann man eine stochatische Simulation durchführen. Dabei wird auf einem 6x6 Spielbrett ein Feld mit zwei Würfeln bestimmt und z.B. ein Geldstück auf diesem Feld umgedreht. Kopf und Zahl stehen für Edukt- und Produktteilchen. Dieses Spiel lässt sich ebenfalls in einen übersichtlichen TigerJython-Code übersetzen, wie unten abgebildet.

Kommentierter Code für TigerJython

Bei diesem Spiel bzw. Simulation ist die Reaktionsgeschwindigkeit in beide Richtungen gleich, so dass nach einiger Zeit gleich viel Edunkt- und Produktteilchen auf dem Brett liegen. Simuliert man weiter, dann fluktuiert die Zahl der Teilchen um den 50%-Wert. Die Gleichgewichtskonstante ist also immer K=1. Will man andere Gleichgewichtskonstanten berücksichtigen, dann müssen sich die Wahrscheinlichkeiten für Hin- und Rückreaktion unterschieden, wie in dem Code unten umgesetzt. In den Zeilen 17/18 ist die Wahrscheinlichkeit für die Hinreaktion festgelegt, in den Zeilen 22/23 für die Rückreaktion.

Kommentierter Code für TigerJython