function [ k, t0, y0 ] = reaction_parameters ( )

%*****************************************************************************80
%
%% reaction_parameters returns parameters for the reaction ODE.
%
%  Licensing:
%
%    This code is distributed under the GNU LGPL license.
%
%  Modified:
%
%    22 January 2021
%
%  Author:
%
%    John Burkardt
%
%  Output:
%
%    real K: the reaction rate.
%
%    real T0: the initial time, in s.
%
%    real Y0[3]: the initial condition at time T0, in Mol/L.
%
  k = 0.02;
  t0 = 0.0;
  a0 = 10.0;
  b0 = 12.0;
  c0 = 3.0;

  y0 = [ a0; b0; c0 ];

  return
end

