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

%*****************************************************************************80
%
%% rule13() returns the pyramid quadrature rule of precision 13.
%
%  Discussion:
%
%    The unit pyramid with square base is the region
%
%      -1 <= X <= 1
%      -1 <= Y <= 1
%       0 <= Z <= 1.
%
%  Licensing:
%
%    This code is distributed under the MIT license.
%
%  Modified:
%
%    15 April 2023
%
%  Author:
%
%    John Burkardt
%
%  Reference:
%
%    Jan Jaskowiec, Natarajan Sukumar,
%    High order symmetric cubature rules for tetrahedra and pyramids,
%    International Journal for Numerical Methods in Engineering,
%    Volume 122, Number 1, pages 148-171, 24 August 2020.
%
%  Output:
%
%    x(n), y(n), z(n): the coordinates of quadrature points.
%
%    w(n): the quadrature weights.
%
  x = [ ...
      0.00000000000000000000, ...
      0.00000000000000000000, ...
      0.00000000000000000000, ...
      0.00000000000000000000, ...
      0.41223892088328983885, ...
      0.00000000000000000000, ...
     -0.41223892088328983885, ...
      0.00000000000000000000, ...
      0.94482408797958883362, ...
      0.00000000000000000000, ...
     -0.94482408797958883362, ...
      0.00000000000000000000, ...
      0.15012306550295120844, ...
      0.00000000000000000000, ...
     -0.15012306550295120844, ...
      0.00000000000000000000, ...
      0.36338283577562219273, ...
      0.00000000000000000000, ...
     -0.36338283577562219273, ...
      0.00000000000000000000, ...
      0.32897122373450704558, ...
      0.00000000000000000000, ...
     -0.32897122373450704558, ...
      0.00000000000000000000, ...
      0.55704935891608620135, ...
      0.00000000000000000000, ...
     -0.55704935891608620135, ...
      0.00000000000000000000, ...
      0.71060354437421247020, ...
      0.00000000000000000000, ...
     -0.71060354437421247020, ...
      0.00000000000000000000, ...
      0.88531748708681257121, ...
      0.00000000000000000000, ...
     -0.88531748708681257121, ...
      0.00000000000000000000, ...
      0.48884754503841615358, ...
     -0.48884754503841615358, ...
      0.48884754503841615358, ...
     -0.48884754503841615358, ...
      0.45031772456672747307, ...
     -0.45031772456672747307, ...
      0.45031772456672747307, ...
     -0.45031772456672747307, ...
      0.65872956915547864476, ...
     -0.65872956915547864476, ...
      0.65872956915547864476, ...
     -0.65872956915547864476, ...
      0.62588271911798842861, ...
     -0.62588271911798842861, ...
      0.62588271911798842861, ...
     -0.62588271911798842861, ...
      0.25178960033193620305, ...
     -0.25178960033193620305, ...
      0.25178960033193620305, ...
     -0.25178960033193620305, ...
      0.91338191280449132492, ...
     -0.91338191280449132492, ...
      0.91338191280449132492, ...
     -0.91338191280449132492, ...
      0.11117861641720541699, ...
     -0.11117861641720541699, ...
      0.11117861641720541699, ...
     -0.11117861641720541699, ...
      0.37814924987375336807, ...
     -0.37814924987375336807, ...
      0.37814924987375336807, ...
     -0.37814924987375336807, ...
      0.23146195691577758913, ...
     -0.23146195691577758913, ...
      0.23146195691577758913, ...
     -0.23146195691577758913, ...
      0.18598667282576661353, ...
     -0.18598667282576661353, ...
      0.18598667282576661353, ...
     -0.18598667282576661353, ...
      0.72941490542892961635, ...
     -0.72941490542892961635, ...
      0.72941490542892961635, ...
     -0.72941490542892961635, ...
      0.80419241777926819825, ...
     -0.80419241777926819825, ...
      0.80419241777926819825, ...
     -0.80419241777926819825, ...
      0.47036082838417470064, ...
     -0.47036082838417470064, ...
      0.47036082838417470064, ...
     -0.47036082838417470064, ...
      0.27535657520011547206, ...
     -0.27535657520011547206, ...
      0.27535657520011547206, ...
     -0.27535657520011547206, ...
      0.25840443662578793660, ...
     -0.25840443662578793660, ...
      0.25840443662578793660, ...
     -0.25840443662578793660, ...
      0.77456301686419115615, ...
      0.43978938971942055369, ...
     -0.77456301686419115615, ...
      0.43978938971942055369, ...
      0.77456301686419115615, ...
     -0.43978938971942055369, ...
     -0.77456301686419115615, ...
     -0.43978938971942055369, ...
      0.72014146551340840752, ...
      0.22174790443082656455, ...
     -0.72014146551340840752, ...
      0.22174790443082656455, ...
      0.72014146551340840752, ...
     -0.22174790443082656455, ...
     -0.72014146551340840752, ...
     -0.22174790443082656455, ...
      0.17509849556381359981, ...
      0.58560797514860796209, ...
     -0.17509849556381359981, ...
      0.58560797514860796209, ...
      0.17509849556381359981, ...
     -0.58560797514860796209, ...
     -0.17509849556381359981, ...
     -0.58560797514860796209, ...
      0.93762547318055500245, ...
      0.73982675325131319610, ...
     -0.93762547318055500245, ...
      0.73982675325131319610, ...
      0.93762547318055500245, ...
     -0.73982675325131319610, ...
     -0.93762547318055500245, ...
     -0.73982675325131319610, ...
      0.92588822421370353677, ...
      0.50002093145521131490, ...
     -0.92588822421370353677, ...
      0.50002093145521131490, ...
      0.92588822421370353677, ...
     -0.50002093145521131490, ...
     -0.92588822421370353677, ...
     -0.50002093145521131490, ...
      0.26820739210416322251, ...
      0.56185198689199644662, ...
     -0.26820739210416322251, ...
      0.56185198689199644662, ...
      0.26820739210416322251, ...
     -0.56185198689199644662, ...
     -0.26820739210416322251, ...
     -0.56185198689199644662, ...
      0.36207959439737902319, ...
      0.79245628070986395830, ...
     -0.36207959439737902319, ...
      0.79245628070986395830, ...
      0.36207959439737902319, ...
     -0.79245628070986395830, ...
     -0.36207959439737902319, ...
     -0.79245628070986395830 ];

  y = [ ...
      0.00000000000000000000, ...
      0.00000000000000000000, ...
      0.00000000000000000000, ...
      0.00000000000000000000, ...
      0.00000000000000000000, ...
      0.41223892088328983885, ...
      0.00000000000000000000, ...
     -0.41223892088328983885, ...
      0.00000000000000000000, ...
      0.94482408797958883362, ...
      0.00000000000000000000, ...
     -0.94482408797958883362, ...
      0.00000000000000000000, ...
      0.15012306550295120844, ...
      0.00000000000000000000, ...
     -0.15012306550295120844, ...
      0.00000000000000000000, ...
      0.36338283577562219273, ...
      0.00000000000000000000, ...
     -0.36338283577562219273, ...
      0.00000000000000000000, ...
      0.32897122373450704558, ...
      0.00000000000000000000, ...
     -0.32897122373450704558, ...
      0.00000000000000000000, ...
      0.55704935891608620135, ...
      0.00000000000000000000, ...
     -0.55704935891608620135, ...
      0.00000000000000000000, ...
      0.71060354437421247020, ...
      0.00000000000000000000, ...
     -0.71060354437421247020, ...
      0.00000000000000000000, ...
      0.88531748708681257121, ...
      0.00000000000000000000, ...
     -0.88531748708681257121, ...
      0.48884754503841615358, ...
      0.48884754503841615358, ...
     -0.48884754503841615358, ...
     -0.48884754503841615358, ...
      0.45031772456672747307, ...
      0.45031772456672747307, ...
     -0.45031772456672747307, ...
     -0.45031772456672747307, ...
      0.65872956915547864476, ...
      0.65872956915547864476, ...
     -0.65872956915547864476, ...
     -0.65872956915547864476, ...
      0.62588271911798842861, ...
      0.62588271911798842861, ...
     -0.62588271911798842861, ...
     -0.62588271911798842861, ...
      0.25178960033193620305, ...
      0.25178960033193620305, ...
     -0.25178960033193620305, ...
     -0.25178960033193620305, ...
      0.91338191280449132492, ...
      0.91338191280449132492, ...
     -0.91338191280449132492, ...
     -0.91338191280449132492, ...
      0.11117861641720541699, ...
      0.11117861641720541699, ...
     -0.11117861641720541699, ...
     -0.11117861641720541699, ...
      0.37814924987375336807, ...
      0.37814924987375336807, ...
     -0.37814924987375336807, ...
     -0.37814924987375336807, ...
      0.23146195691577758913, ...
      0.23146195691577758913, ...
     -0.23146195691577758913, ...
     -0.23146195691577758913, ...
      0.18598667282576661353, ...
      0.18598667282576661353, ...
     -0.18598667282576661353, ...
     -0.18598667282576661353, ...
      0.72941490542892961635, ...
      0.72941490542892961635, ...
     -0.72941490542892961635, ...
     -0.72941490542892961635, ...
      0.80419241777926819825, ...
      0.80419241777926819825, ...
     -0.80419241777926819825, ...
     -0.80419241777926819825, ...
      0.47036082838417470064, ...
      0.47036082838417470064, ...
     -0.47036082838417470064, ...
     -0.47036082838417470064, ...
      0.27535657520011547206, ...
      0.27535657520011547206, ...
     -0.27535657520011547206, ...
     -0.27535657520011547206, ...
      0.25840443662578793660, ...
      0.25840443662578793660, ...
     -0.25840443662578793660, ...
     -0.25840443662578793660, ...
      0.43978938971942055369, ...
      0.77456301686419115615, ...
      0.43978938971942055369, ...
     -0.77456301686419115615, ...
     -0.43978938971942055369, ...
      0.77456301686419115615, ...
     -0.43978938971942055369, ...
     -0.77456301686419115615, ...
      0.22174790443082656455, ...
      0.72014146551340840752, ...
      0.22174790443082656455, ...
     -0.72014146551340840752, ...
     -0.22174790443082656455, ...
      0.72014146551340840752, ...
     -0.22174790443082656455, ...
     -0.72014146551340840752, ...
      0.58560797514860796209, ...
      0.17509849556381359981, ...
      0.58560797514860796209, ...
     -0.17509849556381359981, ...
     -0.58560797514860796209, ...
      0.17509849556381359981, ...
     -0.58560797514860796209, ...
     -0.17509849556381359981, ...
      0.73982675325131319610, ...
      0.93762547318055500245, ...
      0.73982675325131319610, ...
     -0.93762547318055500245, ...
     -0.73982675325131319610, ...
      0.93762547318055500245, ...
     -0.73982675325131319610, ...
     -0.93762547318055500245, ...
      0.50002093145521131490, ...
      0.92588822421370353677, ...
      0.50002093145521131490, ...
     -0.92588822421370353677, ...
     -0.50002093145521131490, ...
      0.92588822421370353677, ...
     -0.50002093145521131490, ...
     -0.92588822421370353677, ...
      0.56185198689199644662, ...
      0.26820739210416322251, ...
      0.56185198689199644662, ...
     -0.26820739210416322251, ...
     -0.56185198689199644662, ...
      0.26820739210416322251, ...
     -0.56185198689199644662, ...
     -0.26820739210416322251, ...
      0.79245628070986395830, ...
      0.36207959439737902319, ...
      0.79245628070986395830, ...
     -0.36207959439737902319, ...
     -0.79245628070986395830, ...
      0.36207959439737902319, ...
     -0.79245628070986395830, ...
     -0.36207959439737902319 ];

  z = [ ...
      0.92974005864256004106, ...
      0.58387036044479012631, ...
      0.10942458150673148309, ...
      0.30043598129646631456, ...
      0.40351765144169587929, ...
      0.40351765144169587929, ...
      0.40351765144169587929, ...
      0.40351765144169587929, ...
      0.00554139647030791119, ...
      0.00554139647030791119, ...
      0.00554139647030791119, ...
      0.00554139647030791119, ...
      0.84979294403286831372, ...
      0.84979294403286831372, ...
      0.84979294403286831372, ...
      0.84979294403286831372, ...
      0.60719504362173426504, ...
      0.60719504362173426504, ...
      0.60719504362173426504, ...
      0.60719504362173426504, ...
      0.02033792783786766978, ...
      0.02033792783786766978, ...
      0.02033792783786766978, ...
      0.02033792783786766978, ...
      0.10601199873730246526, ...
      0.10601199873730246526, ...
      0.10601199873730246526, ...
      0.10601199873730246526, ...
      0.24910419915320258788, ...
      0.24910419915320258788, ...
      0.24910419915320258788, ...
      0.24910419915320258788, ...
      0.08107956240970674855, ...
      0.08107956240970674855, ...
      0.08107956240970674855, ...
      0.08107956240970674855, ...
      0.00950815168106079051, ...
      0.00950815168106079051, ...
      0.00950815168106079051, ...
      0.00950815168106079051, ...
      0.48819360509214093646, ...
      0.48819360509214093646, ...
      0.48819360509214093646, ...
      0.48819360509214093646, ...
      0.27690831156608641805, ...
      0.27690831156608641805, ...
      0.27690831156608641805, ...
      0.27690831156608641805, ...
      0.12248093820428122835, ...
      0.12248093820428122835, ...
      0.12248093820428122835, ...
      0.12248093820428122835, ...
      0.22323057279864513824, ...
      0.22323057279864513824, ...
      0.22323057279864513824, ...
      0.22323057279864513824, ...
      0.01207165002118426068, ...
      0.01207165002118426068, ...
      0.01207165002118426068, ...
      0.01207165002118426068, ...
      0.75518724716690510679, ...
      0.75518724716690510679, ...
      0.75518724716690510679, ...
      0.75518724716690510679, ...
      0.08424825204838617965, ...
      0.08424825204838617965, ...
      0.08424825204838617965, ...
      0.08424825204838617965, ...
      0.63188577231776799081, ...
      0.63188577231776799081, ...
      0.63188577231776799081, ...
      0.63188577231776799081, ...
      0.39892977659234107879, ...
      0.39892977659234107879, ...
      0.39892977659234107879, ...
      0.39892977659234107879, ...
      0.02502635713295603415, ...
      0.02502635713295603415, ...
      0.02502635713295603415, ...
      0.02502635713295603415, ...
      0.12038000765135539738, ...
      0.12038000765135539738, ...
      0.12038000765135539738, ...
      0.12038000765135539738, ...
      0.29270359357039776871, ...
      0.29270359357039776871, ...
      0.29270359357039776871, ...
      0.29270359357039776871, ...
      0.70742748921558806785, ...
      0.70742748921558806785, ...
      0.70742748921558806785, ...
      0.70742748921558806785, ...
      0.52832746509410000169, ...
      0.52832746509410000169, ...
      0.52832746509410000169, ...
      0.52832746509410000169, ...
      0.19538409012904675577, ...
      0.19538409012904675577, ...
      0.19538409012904675577, ...
      0.19538409012904675577, ...
      0.19538409012904675577, ...
      0.19538409012904675577, ...
      0.19538409012904675577, ...
      0.19538409012904675577, ...
      0.02260035188498128386, ...
      0.02260035188498128386, ...
      0.02260035188498128386, ...
      0.02260035188498128386, ...
      0.02260035188498128386, ...
      0.02260035188498128386, ...
      0.02260035188498128386, ...
      0.02260035188498128386, ...
      0.21608520533637087802, ...
      0.21608520533637087802, ...
      0.21608520533637087802, ...
      0.21608520533637087802, ...
      0.21608520533637087802, ...
      0.21608520533637087802, ...
      0.21608520533637087802, ...
      0.21608520533637087802, ...
      0.06232437306564607427, ...
      0.06232437306564607427, ...
      0.06232437306564607427, ...
      0.06232437306564607427, ...
      0.06232437306564607427, ...
      0.06232437306564607427, ...
      0.06232437306564607427, ...
      0.06232437306564607427, ...
      0.01929314718254383776, ...
      0.01929314718254383776, ...
      0.01929314718254383776, ...
      0.01929314718254383776, ...
      0.01929314718254383776, ...
      0.01929314718254383776, ...
      0.01929314718254383776, ...
      0.01929314718254383776, ...
      0.40102466550487736452, ...
      0.40102466550487736452, ...
      0.40102466550487736452, ...
      0.40102466550487736452, ...
      0.40102466550487736452, ...
      0.40102466550487736452, ...
      0.40102466550487736452, ...
      0.40102466550487736452, ...
      0.09200529105614789482, ...
      0.09200529105614789482, ...
      0.09200529105614789482, ...
      0.09200529105614789482, ...
      0.09200529105614789482, ...
      0.09200529105614789482, ...
      0.09200529105614789482, ...
      0.09200529105614789482 ];

  w = [ ...
      0.00120366956557754735, ...
      0.01890734566258640142, ...
      0.02412883548353303431, ...
      0.01430958560741040810, ...
      0.01462233467606358257, ...
      0.01462233467606358257, ...
      0.01462233467606358257, ...
      0.01462233467606358257, ...
      0.00131047133506746131, ...
      0.00131047133506746131, ...
      0.00131047133506746131, ...
      0.00131047133506746131, ...
      0.00131534859459711507, ...
      0.00131534859459711507, ...
      0.00131534859459711507, ...
      0.00131534859459711507, ...
      0.00604172938779333306, ...
      0.00604172938779333306, ...
      0.00604172938779333306, ...
      0.00604172938779333306, ...
      0.00828934369195859742, ...
      0.00828934369195859742, ...
      0.00828934369195859742, ...
      0.00828934369195859742, ...
      0.01293341373858935758, ...
      0.01293341373858935758, ...
      0.01293341373858935758, ...
      0.01293341373858935758, ...
      0.00563885419154383245, ...
      0.00563885419154383245, ...
      0.00563885419154383245, ...
      0.00563885419154383245, ...
      0.00450499139881360439, ...
      0.00450499139881360439, ...
      0.00450499139881360439, ...
      0.00450499139881360439, ...
      0.00363618420610894565, ...
      0.00363618420610894565, ...
      0.00363618420610894565, ...
      0.00363618420610894565, ...
      0.00388842116802182820, ...
      0.00388842116802182820, ...
      0.00388842116802182820, ...
      0.00388842116802182820, ...
      0.00402365320826360111, ...
      0.00402365320826360111, ...
      0.00402365320826360111, ...
      0.00402365320826360111, ...
      0.00955901949244871221, ...
      0.00955901949244871221, ...
      0.00955901949244871221, ...
      0.00955901949244871221, ...
      0.01938656716858040696, ...
      0.01938656716858040696, ...
      0.01938656716858040696, ...
      0.01938656716858040696, ...
      0.00096875488373819618, ...
      0.00096875488373819618, ...
      0.00096875488373819618, ...
      0.00096875488373819618, ...
      0.00502214350286199571, ...
      0.00502214350286199571, ...
      0.00502214350286199571, ...
      0.00502214350286199571, ...
      0.01448769412031481602, ...
      0.01448769412031481602, ...
      0.01448769412031481602, ...
      0.01448769412031481602, ...
      0.00459976431383042356, ...
      0.00459976431383042356, ...
      0.00459976431383042356, ...
      0.00459976431383042356, ...
      0.01107636061349879207, ...
      0.01107636061349879207, ...
      0.01107636061349879207, ...
      0.01107636061349879207, ...
      0.00469260200432727610, ...
      0.00469260200432727610, ...
      0.00469260200432727610, ...
      0.00469260200432727610, ...
      0.00285194982449880072, ...
      0.00285194982449880072, ...
      0.00285194982449880072, ...
      0.00285194982449880072, ...
      0.01323033835935080553, ...
      0.01323033835935080553, ...
      0.01323033835935080553, ...
      0.01323033835935080553, ...
      0.00111972343972882474, ...
      0.00111972343972882474, ...
      0.00111972343972882474, ...
      0.00111972343972882474, ...
      0.00918246462138119605, ...
      0.00918246462138119605, ...
      0.00918246462138119605, ...
      0.00918246462138119605, ...
      0.00456490653701050169, ...
      0.00456490653701050169, ...
      0.00456490653701050169, ...
      0.00456490653701050169, ...
      0.00456490653701050169, ...
      0.00456490653701050169, ...
      0.00456490653701050169, ...
      0.00456490653701050169, ...
      0.00547503483442879964, ...
      0.00547503483442879964, ...
      0.00547503483442879964, ...
      0.00547503483442879964, ...
      0.00547503483442879964, ...
      0.00547503483442879964, ...
      0.00547503483442879964, ...
      0.00547503483442879964, ...
      0.00992261207645742799, ...
      0.00992261207645742799, ...
      0.00992261207645742799, ...
      0.00992261207645742799, ...
      0.00992261207645742799, ...
      0.00992261207645742799, ...
      0.00992261207645742799, ...
      0.00992261207645742799, ...
      0.00110382912156457261, ...
      0.00110382912156457261, ...
      0.00110382912156457261, ...
      0.00110382912156457261, ...
      0.00110382912156457261, ...
      0.00110382912156457261, ...
      0.00110382912156457261, ...
      0.00110382912156457261, ...
      0.00236228293958670472, ...
      0.00236228293958670472, ...
      0.00236228293958670472, ...
      0.00236228293958670472, ...
      0.00236228293958670472, ...
      0.00236228293958670472, ...
      0.00236228293958670472, ...
      0.00236228293958670472, ...
      0.00630666584934871672, ...
      0.00630666584934871672, ...
      0.00630666584934871672, ...
      0.00630666584934871672, ...
      0.00630666584934871672, ...
      0.00630666584934871672, ...
      0.00630666584934871672, ...
      0.00630666584934871672, ...
      0.00675492513102409885, ...
      0.00675492513102409885, ...
      0.00675492513102409885, ...
      0.00675492513102409885, ...
      0.00675492513102409885, ...
      0.00675492513102409885, ...
      0.00675492513102409885, ...
      0.00675492513102409885 ];

  return
end
