function [ gr, gr_conv, gr_count ] = gr_init ( gr_num, par_num )

%*****************************************************************************80
%
%% GR_INIT initializes Gelman-Rubin variables.
%
%  Licensing:
%
%    This code is distributed under the GNU LGPL license.
%
%  Modified:
%
%    22 January 2018
%
%  Author:
%
%    Original FORTRAN90 version by Guannan Zhang.
%    MATLAB version by John Burkardt.
%
%  Input parameters:
%
%    Input, integer GR_NUM, the number of times the Gelman-Rubin
%    statistic may be computed.
%
%    Input, integer PAR_NUM, the number of parameters.
%    1 <= PAR_NUM.
%
%  Output parameters:
%
%    Output, real GR(PAR_NUM,GR_NUM), the Gelman-Rubin statistic.
%
%    Output, logical GR_CONV, the convergence flag.
%
%    Output, integer GR_COUNT, counts the number of generations
%    at which the Gelman-Rubin statistic has been computed.
%
  gr = zeros ( par_num, gr_num );
  gr_conv = false;
  gr_count = 0;

  return
end
