function twins_test ( )

%*****************************************************************************80
%
%% twins_test() tests twins().
%
%  Licensing:
%
%    This code is distributed under the MIT license.
%
%  Modified:
%
%    10 January 2022
%
%  Author:
%
%    John Burkardt
%
  fprintf ( 1, '\n' );
  fprintf ( 1, 'twins_test():\n' );
  fprintf ( 1, '  twins() estimates the probability that 2 twins in 20 students\n' );
  fprintf ( 1, '  will be randomly assigned the same lab group of 4 students.\n' );

  p = twins ( );
  exact = 3 / 19;

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

  return
end

