Performs a chemical reaction from given initial reactant concentrations, with randomness, to compute product concentrations at final time, t_f. The reaction is A → B → C, with kinetic mechanism: a'(t) = −k₁a, b'(t) = −k₂b + k₁a, c'(t) = k₂b ([Winkel, 2011]). Integration is through Python 'scipy' module solver solve_ivp.

^* Caution: Trials > ~1.e+5 can lead to excessive run time and time out failure.

Files involved: P-randKinet.php (this one), RandKinet.php, randKinet.py.

Draws plots of: the kinetics (compositions vs. time; and the random behavior of the resulting composition of the incumbent component.

Python file (download): randKinet.py

Computation structure — A PHP file (this page: 'P-randKinet.php') calls (via 'action=...') an intermediate PHP file ('RandKinet.php'), which (through $_POST) sends the problem data as command line arguments to a Python script. This last does the computing, and makes 'system' call(s) to 'gnuplot', for the final (temporary) webpage.

Keywords: PHP, command line arguments, Python script, gnuplot.

• (Wikipedia) PHP (from Personal Home Page) • Python (programming language)

• Python module 'scipy' solver scipy.integrate.solve_ivp, accessed 02-Oct-2023.

• Winkel, B. (pdf), 2011, International Journal of Mathematical Education in Science and Engineering, 42(1), 37–51. DOI  10.1080/0020739X.2010.500806.

• 1895-09-29: Hotelling, Harold (†1973-12-26, 78 yrs.).