function [ n, x, y, w, s ] = hexagon_stroud_rule2 ( )

%*****************************************************************************80
%
%% hexagon_stroud_rule2() sets the Stroud hexagon quadrature rule #2.
%
%  Licensing:
%
%    This code is distributed under the MIT license. 
%
%  Modified:
%
%    20 April 2019
%
%  Author:
%
%    John Burkardt
%
%  Reference:
%
%    Arthur Stroud,
%    Approximate Calculation of Multiple Integrals,
%    Prentice Hall, 1971.
%
%  Output:
%
%    integer n, the number of points.
%
%    real x(n), y(n), the points.
%
%    real w(n), the weights.
%
%    integer s, the strength.
%
  a = sqrt ( 5.0 / 12.0 );
  b = 1.0 / 4.0;
  z = 0.0;

  n = 4;
  x = [  a; -a;  z;  z ];
  y = [  z;  z;  a; -a ];
  w = [  b;  b;  b;  b ];
  s = 3;

  return
end
