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

%*****************************************************************************80
%
%% hexagon_rule_stroud3 sets the Stroud hexagon quadrature rule #3.
%
%  Licensing:
%
%    This code is distributed under the GNU LGPL license. 
%
%  Modified:
%
%    20 April 2019
%
%  Author:
%
%    John Burkardt
%
%  Reference:
%
%    Arthur Stroud,
%    Approximate Calculation of Multiple Integrals,
%    Prentice Hall, 1971.
%
%  Parameters:
%
%    Output, integer n, the number of points.
%
%    Output, real x(n), y(n), the points.
%
%    Output, real w(n), the weights.
%
%    Output, integer s, the strength.
%
  a = sqrt ( 3.0 ) / 2.0;
  b =  0.5;
  c =  1.0;
  d =  5.0 / 72.0;
  e = 42.0 / 72.0;
  z =  0.0;

  n = 7;
  x = [  z;  c; -c;  b; -b;  b; -b ];
  y = [  z;  z;  z;  a;  a; -a; -a ];
  w = [  e;  d;  d;  d;  d;  d;  d ];
  s = 3;

  return
end
