function [ a, b, c, d, w ] = rule08 ( )

%*****************************************************************************80
%
%% rule08() returns the rule of precision 8.
%
%  Licensing:
%
%    This code is distributed under the MIT license.
%
%  Modified:
%
%    02 April 2023
%
%  Author:
%
%    John Burkardt
%
%  Reference:
%
%    Jan Jaskowiec, Natarajan Sukumar,
%    High order symmetric cubature rules for tetrahedra and pyramids,
%    International Journal of Numerical Methods in Engineering,
%    Volume 122, Number 1, pages 148-171, 24 August 2020.
%
%  Output:
%
%    a(n), b(n), c(n), d(n): the barycentric coordinates of quadrature points.
%
%    w(n): the quadrature weights.
%
  a = [ ...
      0.18681275827072096885, ...
      0.43956172518783709346, ...
      0.18681275827072096885, ...
      0.18681275827072096885, ...
      0.11446240676123049795, ...
      0.65661277971630849226, ...
      0.11446240676123049795, ...
      0.11446240676123049795, ...
      0.31380659227240159659, ...
      0.05858022318279514085, ...
      0.31380659227240159659, ...
      0.31380659227240159659, ...
      0.04457737944884624520, ...
      0.86626786165346125745, ...
      0.04457737944884624520, ...
      0.04457737944884624520, ...
      0.43451650945615760691, ...
      0.06548349054384239309, ...
      0.43451650945615760691, ...
      0.06548349054384239309, ...
      0.43451650945615760691, ...
      0.06548349054384239309, ...
      0.00507332075042159041, ...
      0.20389317466211026586, ...
      0.20389317466211026586, ...
      0.58714032992535780675, ...
      0.58714032992535780675, ...
      0.58714032992535780675, ...
      0.20389317466211026586, ...
      0.00507332075042159041, ...
      0.20389317466211026586, ...
      0.20389317466211026586, ...
      0.00507332075042159041, ...
      0.20389317466211026586, ...
      0.71467692639330493432, ...
      0.02124777966957167377, ...
      0.02124777966957167377, ...
      0.24282751426755175284, ...
      0.24282751426755175284, ...
      0.24282751426755175284, ...
      0.02124777966957167377, ...
      0.71467692639330493432, ...
      0.02124777966957167377, ...
      0.02124777966957167377, ...
      0.71467692639330493432, ...
      0.02124777966957167377 ];

  b = [ ...
      0.18681275827072096885, ...
      0.18681275827072096885, ...
      0.43956172518783709346, ...
      0.18681275827072096885, ...
      0.11446240676123049795, ...
      0.11446240676123049795, ...
      0.65661277971630849226, ...
      0.11446240676123049795, ...
      0.31380659227240159659, ...
      0.31380659227240159659, ...
      0.05858022318279514085, ...
      0.31380659227240159659, ...
      0.04457737944884624520, ...
      0.04457737944884624520, ...
      0.86626786165346125745, ...
      0.04457737944884624520, ...
      0.06548349054384239309, ...
      0.43451650945615760691, ...
      0.43451650945615760691, ...
      0.06548349054384239309, ...
      0.06548349054384239309, ...
      0.43451650945615760691, ...
      0.20389317466211026586, ...
      0.00507332075042159041, ...
      0.20389317466211026586, ...
      0.20389317466211026586, ...
      0.00507332075042159041, ...
      0.20389317466211026586, ...
      0.58714032992535780675, ...
      0.58714032992535780675, ...
      0.58714032992535780675, ...
      0.20389317466211026586, ...
      0.20389317466211026586, ...
      0.00507332075042159041, ...
      0.02124777966957167377, ...
      0.71467692639330493432, ...
      0.02124777966957167377, ...
      0.02124777966957167377, ...
      0.71467692639330493432, ...
      0.02124777966957167377, ...
      0.24282751426755175284, ...
      0.24282751426755175284, ...
      0.24282751426755175284, ...
      0.02124777966957167377, ...
      0.02124777966957167377, ...
      0.71467692639330493432 ];

  c = [ ...
      0.18681275827072096885, ...
      0.18681275827072096885, ...
      0.18681275827072096885, ...
      0.43956172518783709346, ...
      0.11446240676123049795, ...
      0.11446240676123049795, ...
      0.11446240676123049795, ...
      0.65661277971630849226, ...
      0.31380659227240159659, ...
      0.31380659227240159659, ...
      0.31380659227240159659, ...
      0.05858022318279514085, ...
      0.04457737944884624520, ...
      0.04457737944884624520, ...
      0.04457737944884624520, ...
      0.86626786165346125745, ...
      0.06548349054384239309, ...
      0.06548349054384239309, ...
      0.06548349054384239309, ...
      0.43451650945615760691, ...
      0.43451650945615760691, ...
      0.43451650945615760691, ...
      0.20389317466211026586, ...
      0.20389317466211026586, ...
      0.00507332075042159041, ...
      0.20389317466211026586, ...
      0.20389317466211026586, ...
      0.00507332075042159041, ...
      0.20389317466211026586, ...
      0.20389317466211026586, ...
      0.00507332075042159041, ...
      0.58714032992535780675, ...
      0.58714032992535780675, ...
      0.58714032992535780675, ...
      0.02124777966957167377, ...
      0.02124777966957167377, ...
      0.71467692639330493432, ...
      0.02124777966957167377, ...
      0.02124777966957167377, ...
      0.71467692639330493432, ...
      0.02124777966957167377, ...
      0.02124777966957167377, ...
      0.71467692639330493432, ...
      0.24282751426755175284, ...
      0.24282751426755175284, ...
      0.24282751426755175284 ];

  d = [ ...
      0.43956172518783709346, ...
      0.18681275827072096885, ...
      0.18681275827072096885, ...
      0.18681275827072096885, ...
      0.65661277971630849226, ...
      0.11446240676123049795, ...
      0.11446240676123049795, ...
      0.11446240676123049795, ...
      0.05858022318279514085, ...
      0.31380659227240159659, ...
      0.31380659227240159659, ...
      0.31380659227240159659, ...
      0.86626786165346125745, ...
      0.04457737944884624520, ...
      0.04457737944884624520, ...
      0.04457737944884624520, ...
      0.43451650945615760691, ...
      0.43451650945615760691, ...
      0.06548349054384239309, ...
      0.43451650945615760691, ...
      0.06548349054384239309, ...
      0.06548349054384239309, ...
      0.58714032992535780675, ...
      0.58714032992535780675, ...
      0.58714032992535780675, ...
      0.00507332075042159041, ...
      0.20389317466211026586, ...
      0.20389317466211026586, ...
      0.00507332075042159041, ...
      0.20389317466211026586, ...
      0.20389317466211026586, ...
      0.00507332075042159041, ...
      0.20389317466211026586, ...
      0.20389317466211026586, ...
      0.24282751426755175284, ...
      0.24282751426755175284, ...
      0.24282751426755175284, ...
      0.71467692639330493432, ...
      0.02124777966957167377, ...
      0.02124777966957167377, ...
      0.71467692639330493432, ...
      0.02124777966957167377, ...
      0.02124777966957167377, ...
      0.71467692639330493432, ...
      0.02124777966957167377, ...
      0.02124777966957167377 ];

  w = [ ...
      0.04781298803430665545, ...
      0.04781298803430665545, ...
      0.04781298803430665545, ...
      0.04781298803430665545, ...
      0.02972785408679422259, ...
      0.02972785408679422259, ...
      0.02972785408679422259, ...
      0.02972785408679422259, ...
      0.04352732421897526766, ...
      0.04352732421897526766, ...
      0.04352732421897526766, ...
      0.04352732421897526766, ...
      0.00863273602081497497, ...
      0.00863273602081497497, ...
      0.00863273602081497497, ...
      0.00863273602081497497, ...
      0.03689054126986636428, ...
      0.03689054126986636428, ...
      0.03689054126986636428, ...
      0.03689054126986636428, ...
      0.03689054126986636428, ...
      0.03689054126986636428, ...
      0.01449702671085950460, ...
      0.01449702671085950460, ...
      0.01449702671085950460, ...
      0.01449702671085950460, ...
      0.01449702671085950460, ...
      0.01449702671085950460, ...
      0.01449702671085950460, ...
      0.01449702671085950460, ...
      0.01449702671085950460, ...
      0.01449702671085950460, ...
      0.01449702671085950460, ...
      0.01449702671085950460, ...
      0.00715740186724360752, ...
      0.00715740186724360752, ...
      0.00715740186724360752, ...
      0.00715740186724360752, ...
      0.00715740186724360752, ...
      0.00715740186724360752, ...
      0.00715740186724360752, ...
      0.00715740186724360752, ...
      0.00715740186724360752, ...
      0.00715740186724360752, ...
      0.00715740186724360752, ...
      0.00715740186724360752 ];

  return
end
