function double_six_test ( )

%*****************************************************************************80
%
%% double_six_test() tests double_six().
%
%  Licensing:
%
%    This code is distributed under the MIT license.
%
%  Modified:
%
%    10 January 2022
%
%  Author:
%
%    John Burkardt
%
  fprintf ( 1, '\n' );
  fprintf ( 1, 'double_six_test():\n' );
  fprintf ( 1, '  double_six() considers Newton''s double six dice problem.\n' );
  fprintf ( 1, '  What is the expected number of times you must roll a fair die\n' );
  fprintf ( 1, '  in order to observe two consecutive values of six?\n' );

  p = double_six ( );
  p2 = 1.0 / 6.0;
  exact = 1.0 / p2^2 + 1.0 / p;

  fprintf ( 1, '\n' );
  fprintf ( 1, '  Estimate = %g, exact = %g\n', p, exact );

  return
end

