reaction_twoway_ode, a Python code which sets up and solves the ordinary differential equations (ODE) which model a two-way chemical reaction between species W1 and W2.
The reference suggests an experiment in which the Forward Euler method is used, with a timestep of 1/50. The problem is to be solved three times, with the second reaction coefficient, k2, set to 10 (the default), 100, and 1000. The value k2=1000 makes the solution very difficult to compute accurately.
The computer code and data files described and made available on this web page are distributed under the GNU LGPL license.
reaction_twoway_ode is available in a MATLAB version and an Octave version and a Python version.
python_ode, Python codes which sets up various systems of ordinary differential equations (ODE).