function [ x, y, z, w ] = rule08 ( )

%*****************************************************************************80
%
%% rule08() returns the prism rule of precision 8.
%
%  Discussion:
%
%    We suppose we are given a triangular prism P with vertices
%    (1,0,0), (0,1,0), (0,0,0),
%    (1,0,1), (0,1,1), (0,0,1),
%    We call a rule with n points, returning coordinates
%    x, y, z, and weights w.  Then the integral I of f(x,y,z) over P is 
%    approximated by Q as follows:
%
%    Q = volume(P) * sum ( 1 <= i <= n ) w(i) * f(x(i),y(i),z(i)) 
%
%  Licensing:
%
%    This code is distributed under the MIT license.
%
%  Modified:
%
%    02 May 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:
%
%    real x(n), y(n), z(n): the coordinates of quadrature points.
%
%    real w(n): the quadrature weights, which sum to 1.
%
  x = [ ...
          0.3333333333333333, ...
          0.3333333333333333, ...
          0.3333333333333333, ...
          0.3333333333333333, ...
          0.4600889628137106, ...
          0.4600889628137106, ...
          0.0798220743725788, ...
          0.0534123509369071, ...
          0.0534123509369071, ...
          0.8931752981261857, ...
          0.0472387858397694, ...
          0.0472387858397694, ...
          0.9055224283204613, ...
          0.0472387858397694, ...
          0.0472387858397694, ...
          0.9055224283204613, ...
          0.1740616079243704, ...
          0.1740616079243704, ...
          0.6518767841512593, ...
          0.1740616079243704, ...
          0.1740616079243704, ...
          0.6518767841512593, ...
          0.1597492639425890, ...
          0.1597492639425890, ...
          0.6805014721148220, ...
          0.1597492639425890, ...
          0.1597492639425890, ...
          0.6805014721148220, ...
          0.4585690687909513, ...
          0.4585690687909513, ...
          0.0828618624180973, ...
          0.4585690687909513, ...
          0.4585690687909513, ...
          0.0828618624180973, ...
          0.0085881275077590, ...
          0.7285980718010000, ...
          0.0085881275077590, ...
          0.2628138006912410, ...
          0.7285980718010000, ...
          0.2628138006912410, ...
          0.0085881275077590, ...
          0.7285980718010000, ...
          0.0085881275077590, ...
          0.2628138006912410, ...
          0.7285980718010000, ...
          0.2628138006912410 ];

  y = [ ...
          0.3333333333333333, ...
          0.3333333333333333, ...
          0.3333333333333333, ...
          0.3333333333333333, ...
          0.4600889628137106, ...
          0.0798220743725788, ...
          0.4600889628137106, ...
          0.0534123509369071, ...
          0.8931752981261857, ...
          0.0534123509369071, ...
          0.0472387858397694, ...
          0.9055224283204613, ...
          0.0472387858397694, ...
          0.0472387858397694, ...
          0.9055224283204613, ...
          0.0472387858397694, ...
          0.1740616079243704, ...
          0.6518767841512593, ...
          0.1740616079243704, ...
          0.1740616079243704, ...
          0.6518767841512593, ...
          0.1740616079243704, ...
          0.1597492639425890, ...
          0.6805014721148220, ...
          0.1597492639425890, ...
          0.1597492639425890, ...
          0.6805014721148220, ...
          0.1597492639425890, ...
          0.4585690687909513, ...
          0.0828618624180973, ...
          0.4585690687909513, ...
          0.4585690687909513, ...
          0.0828618624180973, ...
          0.4585690687909513, ...
          0.7285980718010000, ...
          0.0085881275077590, ...
          0.2628138006912410, ...
          0.0085881275077590, ...
          0.2628138006912410, ...
          0.7285980718010000, ...
          0.7285980718010000, ...
          0.0085881275077590, ...
          0.2628138006912410, ...
          0.0085881275077590, ...
          0.2628138006912410, ...
          0.7285980718010000 ];

  z = [ ...
          0.6122169330529953, ...
          0.3877830669470047, ...
          0.1590872792145671, ...
          0.8409127207854329, ...
          0.5000000000000000, ...
          0.5000000000000000, ...
          0.5000000000000000, ...
          0.5000000000000000, ...
          0.5000000000000000, ...
          0.5000000000000000, ...
          0.9144077215293400, ...
          0.9144077215293400, ...
          0.9144077215293400, ...
          0.0855922784706599, ...
          0.0855922784706599, ...
          0.0855922784706599, ...
          0.7030421146772952, ...
          0.7030421146772952, ...
          0.7030421146772952, ...
          0.2969578853227048, ...
          0.2969578853227048, ...
          0.2969578853227048, ...
          0.9829325832703272, ...
          0.9829325832703272, ...
          0.9829325832703272, ...
          0.0170674167296728, ...
          0.0170674167296728, ...
          0.0170674167296728, ...
          0.9380979550001272, ...
          0.9380979550001272, ...
          0.9380979550001272, ...
          0.0619020449998729, ...
          0.0619020449998729, ...
          0.0619020449998729, ...
          0.7886751345948129, ...
          0.7886751345948129, ...
          0.7886751345948129, ...
          0.7886751345948129, ...
          0.7886751345948129, ...
          0.7886751345948129, ...
          0.2113248654051871, ...
          0.2113248654051871, ...
          0.2113248654051871, ...
          0.2113248654051871, ...
          0.2113248654051871, ...
          0.2113248654051871 ];

  w = [ ...
          0.0207349366428555, ...
          0.0207349366428555, ...
          0.0510894271121815, ...
          0.0510894271121815, ...
          0.0526997449065072, ...
          0.0526997449065072, ...
          0.0526997449065072, ...
          0.0181522487841497, ...
          0.0181522487841497, ...
          0.0181522487841497, ...
          0.0071453577535386, ...
          0.0071453577535386, ...
          0.0071453577535386, ...
          0.0071453577535386, ...
          0.0071453577535386, ...
          0.0071453577535386, ...
          0.0406292626589893, ...
          0.0406292626589893, ...
          0.0406292626589893, ...
          0.0406292626589893, ...
          0.0406292626589893, ...
          0.0406292626589893, ...
          0.0111430273484653, ...
          0.0111430273484653, ...
          0.0111430273484653, ...
          0.0111430273484653, ...
          0.0111430273484653, ...
          0.0111430273484653, ...
          0.0210865092935865, ...
          0.0210865092935865, ...
          0.0210865092935865, ...
          0.0210865092935865, ...
          0.0210865092935865, ...
          0.0210865092935865, ...
          0.0136475290908731, ...
          0.0136475290908731, ...
          0.0136475290908731, ...
          0.0136475290908731, ...
          0.0136475290908731, ...
          0.0136475290908731, ...
          0.0136475290908731, ...
          0.0136475290908731, ...
          0.0136475290908731, ...
          0.0136475290908731, ...
          0.0136475290908731, ...
          0.0136475290908731 ];

  return
end
