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

%*****************************************************************************80
%
%% hexagon_stroud_rule4() sets the Stroud hexagon quadrature rule #4.
%
%  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 ( 14.0 ) / 5.0;
  b = sqrt ( 14.0 ) / 10.0;
  c = sqrt ( 42.0 ) / 10.0;
  d = 125.0 / 1008.0;
  e = 258.0 / 1008.0;
  z = 0.0;

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

  return
end
