function [ y_equi ] = glycolysis_equilibrium ( )

%*****************************************************************************80
%
%% glycolysis_equilibrium() returns the equilibrium solution of glycolysis_ode().
%
%  Licensing:
%
%    This code is distributed under the MIT license.
%
%  Modified:
%
%    01 January 2022
%
%  Author:
%
%    John Burkardt
%
%  Output:
%
%    real Y_EQUI(2): the equilibrium solution.
%
  [ a, b, t0, y0, tstop ] = glycolysis_parameters ( );

  y_equi = [ b, b/(a+b^2) ];

  return
end

