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

%*****************************************************************************80
%
%% rule13() returns the prism rule of precision 13.
%
%  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.4898439392870481, ...
          0.4898439392870481, ...
          0.0203121214259039, ...
          0.4898439392870481, ...
          0.4898439392870481, ...
          0.0203121214259039, ...
          0.4652402431593082, ...
          0.4652402431593082, ...
          0.0695195136813836, ...
          0.4652402431593082, ...
          0.4652402431593082, ...
          0.0695195136813836, ...
          0.1189161740974999, ...
          0.1189161740974999, ...
          0.7621676518050002, ...
          0.1189161740974999, ...
          0.1189161740974999, ...
          0.7621676518050002, ...
          0.4065656818638698, ...
          0.4065656818638698, ...
          0.1868686362722604, ...
          0.4065656818638698, ...
          0.4065656818638698, ...
          0.1868686362722604, ...
          0.3731243598486834, ...
          0.3731243598486834, ...
          0.2537512803026333, ...
          0.3731243598486834, ...
          0.3731243598486834, ...
          0.2537512803026333, ...
          0.0222257711836755, ...
          0.0222257711836755, ...
          0.9555484576326490, ...
          0.0222257711836755, ...
          0.0222257711836755, ...
          0.9555484576326490, ...
          0.4252710267490136, ...
          0.4252710267490136, ...
          0.1494579465019728, ...
          0.4252710267490136, ...
          0.4252710267490136, ...
          0.1494579465019728, ...
          0.2896929731593648, ...
          0.2896929731593648, ...
          0.4206140536812703, ...
          0.2896929731593648, ...
          0.2896929731593648, ...
          0.4206140536812703, ...
          0.1072868932263340, ...
          0.1072868932263340, ...
          0.7854262135473320, ...
          0.1072868932263340, ...
          0.1072868932263340, ...
          0.7854262135473320, ...
          0.2080125480772502, ...
          0.2080125480772502, ...
          0.5839749038454997, ...
          0.2080125480772502, ...
          0.2080125480772502, ...
          0.5839749038454997, ...
          0.2230030950549982, ...
          0.2230030950549982, ...
          0.5539938098900037, ...
          0.2230030950549982, ...
          0.2230030950549982, ...
          0.5539938098900037, ...
          0.0261743916084927, ...
          0.0261743916084927, ...
          0.9476512167830147, ...
          0.0261743916084927, ...
          0.0261743916084927, ...
          0.9476512167830147, ...
          0.0069141701592508, ...
          0.2742659465495736, ...
          0.0069141701592508, ...
          0.7188198832911755, ...
          0.2742659465495736, ...
          0.7188198832911755, ...
          0.1210652398684976, ...
          0.0234237758449634, ...
          0.1210652398684976, ...
          0.8555109842865389, ...
          0.0234237758449634, ...
          0.8555109842865389, ...
          0.0935360177546358, ...
          0.5415308571572887, ...
          0.0935360177546358, ...
          0.3649331250880756, ...
          0.5415308571572887, ...
          0.3649331250880756, ...
          0.0190809537819619, ...
          0.4425451813256520, ...
          0.0190809537819619, ...
          0.5383738648923861, ...
          0.4425451813256520, ...
          0.5383738648923861, ...
          0.0190809537819619, ...
          0.4425451813256520, ...
          0.0190809537819619, ...
          0.5383738648923861, ...
          0.4425451813256520, ...
          0.5383738648923861, ...
          0.2923396969545124, ...
          0.0181776522602859, ...
          0.2923396969545124, ...
          0.6894826507852017, ...
          0.0181776522602859, ...
          0.6894826507852017, ...
          0.2923396969545124, ...
          0.0181776522602859, ...
          0.2923396969545124, ...
          0.6894826507852017, ...
          0.0181776522602859, ...
          0.6894826507852017, ...
          0.0212309381671597, ...
          0.8615902483468726, ...
          0.0212309381671597, ...
          0.1171788134859677, ...
          0.8615902483468726, ...
          0.1171788134859677, ...
          0.0212309381671597, ...
          0.8615902483468726, ...
          0.0212309381671597, ...
          0.1171788134859677, ...
          0.8615902483468726, ...
          0.1171788134859677, ...
          0.2541452627735283, ...
          0.0752550719401233, ...
          0.2541452627735283, ...
          0.6705996652863484, ...
          0.0752550719401233, ...
          0.6705996652863484, ...
          0.2541452627735283, ...
          0.0752550719401233, ...
          0.2541452627735283, ...
          0.6705996652863484, ...
          0.0752550719401233, ...
          0.6705996652863484, ...
          0.1414034499801619, ...
          0.0248530932165772, ...
          0.1414034499801619, ...
          0.8337434568032609, ...
          0.0248530932165772, ...
          0.8337434568032609, ...
          0.1414034499801619, ...
          0.0248530932165772, ...
          0.1414034499801619, ...
          0.8337434568032609, ...
          0.0248530932165772, ...
          0.8337434568032609, ...
          0.5920519581168684, ...
          0.1079442302776815, ...
          0.5920519581168684, ...
          0.3000038116054501, ...
          0.1079442302776815, ...
          0.3000038116054501, ...
          0.5920519581168684, ...
          0.1079442302776815, ...
          0.5920519581168684, ...
          0.3000038116054501, ...
          0.1079442302776815, ...
          0.3000038116054501 ];

  y = [ ...
          0.4898439392870481, ...
          0.0203121214259039, ...
          0.4898439392870481, ...
          0.4898439392870481, ...
          0.0203121214259039, ...
          0.4898439392870481, ...
          0.4652402431593082, ...
          0.0695195136813836, ...
          0.4652402431593082, ...
          0.4652402431593082, ...
          0.0695195136813836, ...
          0.4652402431593082, ...
          0.1189161740974999, ...
          0.7621676518050002, ...
          0.1189161740974999, ...
          0.1189161740974999, ...
          0.7621676518050002, ...
          0.1189161740974999, ...
          0.4065656818638698, ...
          0.1868686362722604, ...
          0.4065656818638698, ...
          0.4065656818638698, ...
          0.1868686362722604, ...
          0.4065656818638698, ...
          0.3731243598486834, ...
          0.2537512803026333, ...
          0.3731243598486834, ...
          0.3731243598486834, ...
          0.2537512803026333, ...
          0.3731243598486834, ...
          0.0222257711836755, ...
          0.9555484576326490, ...
          0.0222257711836755, ...
          0.0222257711836755, ...
          0.9555484576326490, ...
          0.0222257711836755, ...
          0.4252710267490136, ...
          0.1494579465019728, ...
          0.4252710267490136, ...
          0.4252710267490136, ...
          0.1494579465019728, ...
          0.4252710267490136, ...
          0.2896929731593648, ...
          0.4206140536812703, ...
          0.2896929731593648, ...
          0.2896929731593648, ...
          0.4206140536812703, ...
          0.2896929731593648, ...
          0.1072868932263340, ...
          0.7854262135473320, ...
          0.1072868932263340, ...
          0.1072868932263340, ...
          0.7854262135473320, ...
          0.1072868932263340, ...
          0.2080125480772502, ...
          0.5839749038454997, ...
          0.2080125480772502, ...
          0.2080125480772502, ...
          0.5839749038454997, ...
          0.2080125480772502, ...
          0.2230030950549982, ...
          0.5539938098900037, ...
          0.2230030950549982, ...
          0.2230030950549982, ...
          0.5539938098900037, ...
          0.2230030950549982, ...
          0.0261743916084927, ...
          0.9476512167830147, ...
          0.0261743916084927, ...
          0.0261743916084927, ...
          0.9476512167830147, ...
          0.0261743916084927, ...
          0.2742659465495736, ...
          0.0069141701592508, ...
          0.7188198832911755, ...
          0.0069141701592508, ...
          0.7188198832911755, ...
          0.2742659465495736, ...
          0.0234237758449634, ...
          0.1210652398684976, ...
          0.8555109842865389, ...
          0.1210652398684976, ...
          0.8555109842865389, ...
          0.0234237758449634, ...
          0.5415308571572887, ...
          0.0935360177546358, ...
          0.3649331250880756, ...
          0.0935360177546358, ...
          0.3649331250880756, ...
          0.5415308571572887, ...
          0.4425451813256520, ...
          0.0190809537819619, ...
          0.5383738648923861, ...
          0.0190809537819619, ...
          0.5383738648923861, ...
          0.4425451813256520, ...
          0.4425451813256520, ...
          0.0190809537819619, ...
          0.5383738648923861, ...
          0.0190809537819619, ...
          0.5383738648923861, ...
          0.4425451813256520, ...
          0.0181776522602859, ...
          0.2923396969545124, ...
          0.6894826507852017, ...
          0.2923396969545124, ...
          0.6894826507852017, ...
          0.0181776522602859, ...
          0.0181776522602859, ...
          0.2923396969545124, ...
          0.6894826507852017, ...
          0.2923396969545124, ...
          0.6894826507852017, ...
          0.0181776522602859, ...
          0.8615902483468726, ...
          0.0212309381671597, ...
          0.1171788134859677, ...
          0.0212309381671597, ...
          0.1171788134859677, ...
          0.8615902483468726, ...
          0.8615902483468726, ...
          0.0212309381671597, ...
          0.1171788134859677, ...
          0.0212309381671597, ...
          0.1171788134859677, ...
          0.8615902483468726, ...
          0.0752550719401233, ...
          0.2541452627735283, ...
          0.6705996652863484, ...
          0.2541452627735283, ...
          0.6705996652863484, ...
          0.0752550719401233, ...
          0.0752550719401233, ...
          0.2541452627735283, ...
          0.6705996652863484, ...
          0.2541452627735283, ...
          0.6705996652863484, ...
          0.0752550719401233, ...
          0.0248530932165772, ...
          0.1414034499801619, ...
          0.8337434568032609, ...
          0.1414034499801619, ...
          0.8337434568032609, ...
          0.0248530932165772, ...
          0.0248530932165772, ...
          0.1414034499801619, ...
          0.8337434568032609, ...
          0.1414034499801619, ...
          0.8337434568032609, ...
          0.0248530932165772, ...
          0.1079442302776815, ...
          0.5920519581168684, ...
          0.3000038116054501, ...
          0.5920519581168684, ...
          0.3000038116054501, ...
          0.1079442302776815, ...
          0.1079442302776815, ...
          0.5920519581168684, ...
          0.3000038116054501, ...
          0.5920519581168684, ...
          0.3000038116054501, ...
          0.1079442302776815 ];

  z = [ ...
          0.9628749136752521, ...
          0.9628749136752521, ...
          0.9628749136752521, ...
          0.0371250863247479, ...
          0.0371250863247479, ...
          0.0371250863247479, ...
          0.1596095438845206, ...
          0.1596095438845206, ...
          0.1596095438845206, ...
          0.8403904561154794, ...
          0.8403904561154794, ...
          0.8403904561154794, ...
          0.9056527042389668, ...
          0.9056527042389668, ...
          0.9056527042389668, ...
          0.0943472957610332, ...
          0.0943472957610332, ...
          0.0943472957610332, ...
          0.5096252745686303, ...
          0.5096252745686303, ...
          0.5096252745686303, ...
          0.4903747254313698, ...
          0.4903747254313698, ...
          0.4903747254313698, ...
          0.9622105990610692, ...
          0.9622105990610692, ...
          0.9622105990610692, ...
          0.0377894009389309, ...
          0.0377894009389309, ...
          0.0377894009389309, ...
          0.6451848165106220, ...
          0.6451848165106220, ...
          0.6451848165106220, ...
          0.3548151834893781, ...
          0.3548151834893781, ...
          0.3548151834893781, ...
          0.8358323923694781, ...
          0.8358323923694781, ...
          0.8358323923694781, ...
          0.1641676076305219, ...
          0.1641676076305219, ...
          0.1641676076305219, ...
          0.6970082735594627, ...
          0.6970082735594627, ...
          0.6970082735594627, ...
          0.3029917264405373, ...
          0.3029917264405373, ...
          0.3029917264405373, ...
          0.5806209783551364, ...
          0.5806209783551364, ...
          0.5806209783551364, ...
          0.4193790216448636, ...
          0.4193790216448636, ...
          0.4193790216448636, ...
          0.6310517191712730, ...
          0.6310517191712730, ...
          0.6310517191712730, ...
          0.3689482808287270, ...
          0.3689482808287270, ...
          0.3689482808287270, ...
          0.8758646517044397, ...
          0.8758646517044397, ...
          0.8758646517044397, ...
          0.1241353482955603, ...
          0.1241353482955603, ...
          0.1241353482955603, ...
          0.9265272967069804, ...
          0.9265272967069804, ...
          0.9265272967069804, ...
          0.0734727032930196, ...
          0.0734727032930196, ...
          0.0734727032930196, ...
          0.5000000000000000, ...
          0.5000000000000000, ...
          0.5000000000000000, ...
          0.5000000000000000, ...
          0.5000000000000000, ...
          0.5000000000000000, ...
          0.5000000000000000, ...
          0.5000000000000000, ...
          0.5000000000000000, ...
          0.5000000000000000, ...
          0.5000000000000000, ...
          0.5000000000000000, ...
          0.5000000000000000, ...
          0.5000000000000000, ...
          0.5000000000000000, ...
          0.5000000000000000, ...
          0.5000000000000000, ...
          0.5000000000000000, ...
          0.6831055612129073, ...
          0.6831055612129073, ...
          0.6831055612129073, ...
          0.6831055612129073, ...
          0.6831055612129073, ...
          0.6831055612129073, ...
          0.3168944387870927, ...
          0.3168944387870927, ...
          0.3168944387870927, ...
          0.3168944387870927, ...
          0.3168944387870927, ...
          0.3168944387870927, ...
          0.8946991133344783, ...
          0.8946991133344783, ...
          0.8946991133344783, ...
          0.8946991133344783, ...
          0.8946991133344783, ...
          0.8946991133344783, ...
          0.1053008866655217, ...
          0.1053008866655217, ...
          0.1053008866655217, ...
          0.1053008866655217, ...
          0.1053008866655217, ...
          0.1053008866655217, ...
          0.7783161730088111, ...
          0.7783161730088111, ...
          0.7783161730088111, ...
          0.7783161730088111, ...
          0.7783161730088111, ...
          0.7783161730088111, ...
          0.2216838269911889, ...
          0.2216838269911889, ...
          0.2216838269911889, ...
          0.2216838269911889, ...
          0.2216838269911889, ...
          0.2216838269911889, ...
          0.2838393063545257, ...
          0.2838393063545257, ...
          0.2838393063545257, ...
          0.2838393063545257, ...
          0.2838393063545257, ...
          0.2838393063545257, ...
          0.7161606936454743, ...
          0.7161606936454743, ...
          0.7161606936454743, ...
          0.7161606936454743, ...
          0.7161606936454743, ...
          0.7161606936454743, ...
          0.9877489005025584, ...
          0.9877489005025584, ...
          0.9877489005025584, ...
          0.9877489005025584, ...
          0.9877489005025584, ...
          0.9877489005025584, ...
          0.0122510994974416, ...
          0.0122510994974416, ...
          0.0122510994974416, ...
          0.0122510994974416, ...
          0.0122510994974416, ...
          0.0122510994974416, ...
          0.9801022671238835, ...
          0.9801022671238835, ...
          0.9801022671238835, ...
          0.9801022671238835, ...
          0.9801022671238835, ...
          0.9801022671238835, ...
          0.0198977328761166, ...
          0.0198977328761166, ...
          0.0198977328761166, ...
          0.0198977328761166, ...
          0.0198977328761166, ...
          0.0198977328761166 ];

  w = [ ...
          0.0026290597067936, ...
          0.0026290597067936, ...
          0.0026290597067936, ...
          0.0026290597067936, ...
          0.0026290597067936, ...
          0.0026290597067936, ...
          0.0059734237493126, ...
          0.0059734237493126, ...
          0.0059734237493126, ...
          0.0059734237493126, ...
          0.0059734237493126, ...
          0.0059734237493126, ...
          0.0080935494241557, ...
          0.0080935494241557, ...
          0.0080935494241557, ...
          0.0080935494241557, ...
          0.0080935494241557, ...
          0.0080935494241557, ...
          0.0068191162590357, ...
          0.0068191162590357, ...
          0.0068191162590357, ...
          0.0068191162590357, ...
          0.0068191162590357, ...
          0.0068191162590357, ...
          0.0049334513687730, ...
          0.0049334513687730, ...
          0.0049334513687730, ...
          0.0049334513687730, ...
          0.0049334513687730, ...
          0.0049334513687730, ...
          0.0020147080132884, ...
          0.0020147080132884, ...
          0.0020147080132884, ...
          0.0020147080132884, ...
          0.0020147080132884, ...
          0.0020147080132884, ...
          0.0123357707370753, ...
          0.0123357707370753, ...
          0.0123357707370753, ...
          0.0123357707370753, ...
          0.0123357707370753, ...
          0.0123357707370753, ...
          0.0113647943376711, ...
          0.0113647943376711, ...
          0.0113647943376711, ...
          0.0113647943376711, ...
          0.0113647943376711, ...
          0.0113647943376711, ...
          0.0066668649606825, ...
          0.0066668649606825, ...
          0.0066668649606825, ...
          0.0066668649606825, ...
          0.0066668649606825, ...
          0.0066668649606825, ...
          0.0119452201363976, ...
          0.0119452201363976, ...
          0.0119452201363976, ...
          0.0119452201363976, ...
          0.0119452201363976, ...
          0.0119452201363976, ...
          0.0106878809406112, ...
          0.0106878809406112, ...
          0.0106878809406112, ...
          0.0106878809406112, ...
          0.0106878809406112, ...
          0.0106878809406112, ...
          0.0016663443274632, ...
          0.0016663443274632, ...
          0.0016663443274632, ...
          0.0016663443274632, ...
          0.0016663443274632, ...
          0.0016663443274632, ...
          0.0045120516170652, ...
          0.0045120516170652, ...
          0.0045120516170652, ...
          0.0045120516170652, ...
          0.0045120516170652, ...
          0.0045120516170652, ...
          0.0034887068394433, ...
          0.0034887068394433, ...
          0.0034887068394433, ...
          0.0034887068394433, ...
          0.0034887068394433, ...
          0.0034887068394433, ...
          0.0111803506944158, ...
          0.0111803506944158, ...
          0.0111803506944158, ...
          0.0111803506944158, ...
          0.0111803506944158, ...
          0.0111803506944158, ...
          0.0042052602324505, ...
          0.0042052602324505, ...
          0.0042052602324505, ...
          0.0042052602324505, ...
          0.0042052602324505, ...
          0.0042052602324505, ...
          0.0042052602324505, ...
          0.0042052602324505, ...
          0.0042052602324505, ...
          0.0042052602324505, ...
          0.0042052602324505, ...
          0.0042052602324505, ...
          0.0043767917497427, ...
          0.0043767917497427, ...
          0.0043767917497427, ...
          0.0043767917497427, ...
          0.0043767917497427, ...
          0.0043767917497427, ...
          0.0043767917497427, ...
          0.0043767917497427, ...
          0.0043767917497427, ...
          0.0043767917497427, ...
          0.0043767917497427, ...
          0.0043767917497427, ...
          0.0044067550419370, ...
          0.0044067550419370, ...
          0.0044067550419370, ...
          0.0044067550419370, ...
          0.0044067550419370, ...
          0.0044067550419370, ...
          0.0044067550419370, ...
          0.0044067550419370, ...
          0.0044067550419370, ...
          0.0044067550419370, ...
          0.0044067550419370, ...
          0.0044067550419370, ...
          0.0120350326772681, ...
          0.0120350326772681, ...
          0.0120350326772681, ...
          0.0120350326772681, ...
          0.0120350326772681, ...
          0.0120350326772681, ...
          0.0120350326772681, ...
          0.0120350326772681, ...
          0.0120350326772681, ...
          0.0120350326772681, ...
          0.0120350326772681, ...
          0.0120350326772681, ...
          0.0013986921884649, ...
          0.0013986921884649, ...
          0.0013986921884649, ...
          0.0013986921884649, ...
          0.0013986921884649, ...
          0.0013986921884649, ...
          0.0013986921884649, ...
          0.0013986921884649, ...
          0.0013986921884649, ...
          0.0013986921884649, ...
          0.0013986921884649, ...
          0.0013986921884649, ...
          0.0047551548873781, ...
          0.0047551548873781, ...
          0.0047551548873781, ...
          0.0047551548873781, ...
          0.0047551548873781, ...
          0.0047551548873781, ...
          0.0047551548873781, ...
          0.0047551548873781, ...
          0.0047551548873781, ...
          0.0047551548873781, ...
          0.0047551548873781, ...
          0.0047551548873781 ];

  return
end
