function [ k1, k2, t0, y0 ] = reaction_twoway_parameters ( )

%*****************************************************************************80
%
%% reaction_twoway_parameters returns parameters for the twoway reaction ODE.
%
%  Licensing:
%
%    This code is distributed under the GNU LGPL license.
%
%  Modified:
%
%    31 January 2021
%
%  Author:
%
%    John Burkardt
%
%  Reference:
%
%    Willem Hundsdorfer, Jan Verwer,
%    Numerical solution of time-dependent advection-diffusion-reaction equations,
%    Springer, 2003
%    ISBN: 978-3-662-09017-6
%
%  Output:
%
%    real K1, K2: the reaction rates.
%
%    real T0: the initial time.
%
%    real Y0[2]: the initial conditions.
%
  k1 = 1.0;
  k2 = 10.0;
  t0 = 0.0;
  w1 = 0.1;
  w2 = 0.9;

  y0 = [ w1; w2 ];

  return
end

