function [ a, b, c, w ] = rule17 ( )

%*****************************************************************************80
%
%% rule17() returns the rule of precision 17.
%
%  Licensing:
%
%    This code is distributed under the MIT license.
%
%  Modified:
%
%    09 June 2023
%
%  Author:
%
%    John Burkardt
%
%  Reference:
%
%    Hong Xiao, Zydrunas Gimbutas,
%    A numerical algorithm for the construction of efficient
%    quadrature rules in two and higher dimensions, 
%    Computers and Mathematics with Applications,
%    Volume 59, 2010, pages 663-676.
%
%  Output:
%
%    real a(n), b(n), c(n): the barycentric coordinates of quadrature points.
%
%    real w(n): the quadrature weights.
%
  a = [ ...
          0.4171034443615992, ...
          0.4171034443615992, ...
          0.1657931112768016, ...
          0.0725054707990024, ...
          0.9159193532978169, ...
          0.0115751759031809, ...
          0.9159193532978169, ...
          0.0725054707990024, ...
          0.0115751759031804, ...
          0.1803581162663705, ...
          0.6392837674672587, ...
          0.1803581162663708, ...
          0.4154754592952291, ...
          0.5712948679446841, ...
          0.0132296727600869, ...
          0.5712948679446841, ...
          0.4154754592952290, ...
          0.0132296727600868, ...
          0.2717918700553548, ...
          0.7150722591106424, ...
          0.0131358708340029, ...
          0.7150722591106424, ...
          0.2717918700553550, ...
          0.0131358708340026, ...
          0.2992189424769703, ...
          0.5432755795961597, ...
          0.1575054779268699, ...
          0.5432755795961597, ...
          0.2992189424769702, ...
          0.1575054779268698, ...
          0.2857065024365866, ...
          0.4285869951268267, ...
          0.2857065024365867, ...
          0.3062815917461865, ...
          0.6263690303864522, ...
          0.0673493778673613, ...
          0.6263690303864522, ...
          0.3062815917461865, ...
          0.0673493778673612, ...
          0.1687225134952594, ...
          0.7532351459364581, ...
          0.0780423405682825, ...
          0.7532351459364581, ...
          0.1687225134952594, ...
          0.0780423405682823, ...
          0.0666540634795968, ...
          0.8666918730408060, ...
          0.0666540634795971, ...
          0.1591922874727927, ...
          0.8247900701650880, ...
          0.0160176423621194, ...
          0.8247900701650880, ...
          0.1591922874727927, ...
          0.0160176423621192, ...
          0.0147554916607538, ...
          0.9704890166784920, ...
          0.0147554916607542, ...
          0.4655978716188903, ...
          0.4655978716188902, ...
          0.0688042567622194 ];

  b = [ ...
          0.1657931112768017, ...
          0.4171034443615993, ...
          0.4171034443615994, ...
          0.0115751759031807, ...
          0.0725054707990024, ...
          0.9159193532978170, ...
          0.0115751759031807, ...
          0.9159193532978170, ...
          0.0725054707990027, ...
          0.1803581162663707, ...
          0.1803581162663706, ...
          0.6392837674672588, ...
          0.0132296727600870, ...
          0.4154754592952291, ...
          0.5712948679446843, ...
          0.0132296727600870, ...
          0.5712948679446841, ...
          0.4154754592952293, ...
          0.0131358708340028, ...
          0.2717918700553549, ...
          0.7150722591106425, ...
          0.0131358708340028, ...
          0.7150722591106424, ...
          0.2717918700553552, ...
          0.1575054779268700, ...
          0.2992189424769703, ...
          0.5432755795961599, ...
          0.1575054779268700, ...
          0.5432755795961599, ...
          0.2992189424769705, ...
          0.2857065024365867, ...
          0.2857065024365867, ...
          0.4285869951268269, ...
          0.0673493778673613, ...
          0.3062815917461865, ...
          0.6263690303864524, ...
          0.0673493778673613, ...
          0.6263690303864524, ...
          0.3062815917461867, ...
          0.0780423405682825, ...
          0.1687225134952594, ...
          0.7532351459364584, ...
          0.0780423405682825, ...
          0.7532351459364582, ...
          0.1687225134952597, ...
          0.0666540634795971, ...
          0.0666540634795968, ...
          0.8666918730408062, ...
          0.0160176423621193, ...
          0.1591922874727927, ...
          0.8247900701650882, ...
          0.0160176423621193, ...
          0.8247900701650882, ...
          0.1591922874727930, ...
          0.0147554916607541, ...
          0.0147554916607538, ...
          0.9704890166784922, ...
          0.0688042567622194, ...
          0.4655978716188904, ...
          0.4655978716188905 ];

  c = [ ...
          0.4171034443615991, ...
          0.1657931112768015, ...
          0.4171034443615990, ...
          0.9159193532978169, ...
          0.0115751759031807, ...
          0.0725054707990022, ...
          0.0725054707990024, ...
          0.0115751759031805, ...
          0.9159193532978168, ...
          0.6392837674672587, ...
          0.1803581162663707, ...
          0.1803581162663704, ...
          0.5712948679446841, ...
          0.0132296727600869, ...
          0.4154754592952288, ...
          0.4154754592952289, ...
          0.0132296727600869, ...
          0.5712948679446839, ...
          0.7150722591106424, ...
          0.0131358708340027, ...
          0.2717918700553547, ...
          0.2717918700553548, ...
          0.0131358708340026, ...
          0.7150722591106422, ...
          0.5432755795961598, ...
          0.1575054779268700, ...
          0.2992189424769702, ...
          0.2992189424769703, ...
          0.1575054779268699, ...
          0.5432755795961597, ...
          0.4285869951268267, ...
          0.2857065024365867, ...
          0.2857065024365865, ...
          0.6263690303864522, ...
          0.0673493778673613, ...
          0.3062815917461863, ...
          0.3062815917461865, ...
          0.0673493778673611, ...
          0.6263690303864521, ...
          0.7532351459364580, ...
          0.0780423405682825, ...
          0.1687225134952591, ...
          0.1687225134952594, ...
          0.0780423405682823, ...
          0.7532351459364579, ...
          0.8666918730408061, ...
          0.0666540634795972, ...
          0.0666540634795967, ...
          0.8247900701650880, ...
          0.0160176423621193, ...
          0.1591922874727923, ...
          0.1591922874727927, ...
          0.0160176423621191, ...
          0.8247900701650878, ...
          0.9704890166784921, ...
          0.0147554916607542, ...
          0.0147554916607536, ...
          0.4655978716188903, ...
          0.0688042567622194, ...
          0.4655978716188902 ];

  w = [ ...
          0.0273109265281021, ...
          0.0273109265281021, ...
          0.0273109265281021, ...
          0.0045843484017359, ...
          0.0045843484017359, ...
          0.0045843484017359, ...
          0.0045843484017359, ...
          0.0045843484017359, ...
          0.0045843484017359, ...
          0.0263126305880180, ...
          0.0263126305880180, ...
          0.0263126305880180, ...
          0.0103984399558395, ...
          0.0103984399558395, ...
          0.0103984399558395, ...
          0.0103984399558395, ...
          0.0103984399558395, ...
          0.0103984399558395, ...
          0.0086922145010012, ...
          0.0086922145010012, ...
          0.0086922145010012, ...
          0.0086922145010012, ...
          0.0086922145010012, ...
          0.0086922145010012, ...
          0.0261716259353370, ...
          0.0261716259353370, ...
          0.0261716259353370, ...
          0.0261716259353370, ...
          0.0261716259353370, ...
          0.0261716259353370, ...
          0.0377162371527953, ...
          0.0377162371527953, ...
          0.0377162371527953, ...
          0.0224877725466911, ...
          0.0224877725466911, ...
          0.0224877725466911, ...
          0.0224877725466911, ...
          0.0224877725466911, ...
          0.0224877725466911, ...
          0.0205578983204545, ...
          0.0205578983204545, ...
          0.0205578983204545, ...
          0.0205578983204545, ...
          0.0205578983204545, ...
          0.0205578983204545, ...
          0.0124590008023054, ...
          0.0124590008023054, ...
          0.0124590008023054, ...
          0.0079783002059296, ...
          0.0079783002059296, ...
          0.0079783002059296, ...
          0.0079783002059296, ...
          0.0079783002059296, ...
          0.0079783002059296, ...
          0.0027738875776376, ...
          0.0027738875776376, ...
          0.0027738875776376, ...
          0.0250194509504974, ...
          0.0250194509504974, ...
          0.0250194509504974 ];

  return
end
