reaction_twoway_ode

reaction_twoway_ode, a MATLAB 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.

Licensing:

The computer code and data files described and made available on this web page are distributed under the MIT license

Languages:

reaction_twoway_ode is available in a MATLAB version and an Octave version and a Python version.

Related Data and codes:

reaction_twoway_ode_test

matlab_ode, MATLAB codes which set up various ordinary differential equations (ODE).

References:

1. Willem Hundsdorfer, Jan Verwer,
   Numerical solution of time-dependent advection-diffusion-reaction equations,
   Springer, 2003
   ISBN: 978-3-662-09017-6

Source Code:

• reaction_twoway_conserved.m, returns the value of a quantity that should be conserved.
• reaction_twoway_deriv.m, returns the right hand side of the ODE.
• reaction_twoway_exact.m, returns the exact solution of the ODE.
• reaction_twoway_parameters.m, returns parameters of the ODE.