function [ x, y, w ] = rule21 ( )

%*****************************************************************************80
%
%% rule21() returns the rule of degree 21.
%
%  Discussion:
%
%    Order 21 (87 pts)
%    1/6 data for 21-th order quadrature with 19 nodes.
%
%  Licensing:
%
%    This code is distributed under the MIT license.
%
%  Modified:
%
%    26 June 2014
%
%  Author:
%
%    Original FORTRAN77 version by Hong Xiao, Zydrunas Gimbutas.
%    This version by John Burkardt.
%
%  Output:
%
%    real X(*), Y(*), the coordinates of the nodes.
%
%    real W(*), the weights.
%
  x = [ ...
       0.00000000000000000000000000000000, ...
      -0.21632544850764896966019269484948, ...
       0.00000000000000000000000000000000, ...
      -0.51732142577254394519688365999304, ...
       0.00000000000000000000000000000000, ...
      -0.23848997800485075950680771153334, ...
      -0.49726138663597772473398444337837, ...
      -0.72412950658754887280604342699955, ...
      -0.25333184867896101624130032052433, ...
       0.00000000000000000000000000000000, ...
       0.00000000000000000000000000000000, ...
      -0.46634273354592047897176529761299, ...
       0.00000000000000000000000000000000, ...
       0.00000000000000000000000000000000, ...
      -0.69288931850735639864308358232905, ...
       0.00000000000000000000000000000000, ...
      -0.26783854808923110943600296777461, ...
       0.00000000000000000000000000000000, ...
      -0.87825896621183219776432484250545 ];
  y = [ ... ...
       0.11915504280142053470339944811982, ...
      -0.22176792985285491406799950307586, ...
      -0.56698524576800618707986473265874, ...
      -0.56534402375875828330832096864561, ...
      -0.24347179507408359903990857698094, ...
      -0.36295805875690062325727639767329, ...
      -0.50981524806986605698152573307478, ...
      -0.56083300395461387753780504376756, ...
      -0.48546117771023431955108992405636, ...
       0.74251201445344514748786620220395, ...
       0.49552108269676082383881281439860, ...
      -0.40335187439348537560379590829971, ...
      -0.51360341626925635313747432894532, ...
      -0.40349668011940382629268981404988, ...
      -0.49169851113649364618692601028995, ...
       0.96892916731392889311013905422251, ...
      -0.55958871884460214682886011537049, ...
       0.11174875776997437235819388085953E+01, ...
      -0.55950684866453675407246048332841 ];
  w = [ ... ...
       0.14115632059803205472519666814894E-01, ...
       0.23025262548389042223555118865658E-01, ...
       0.29202561691086176915750601020522E-02, ...
       0.55355656065861078854118874771594E-02, ...
       0.15135314836891179901866111095885E-01, ...
       0.24278255412699264374036359934954E-01, ...
       0.13779168816408564476674209048439E-01, ...
       0.58970817066310603066620906137064E-02, ...
       0.19083475799818969202994305068534E-01, ...
       0.89861747152080313144225770269461E-02, ...
       0.12802269114319752046734087041489E-01, ...
       0.20930889409056329487663117084041E-01, ...
       0.80375338997864963118800061347869E-02, ...
       0.12907050562519806085002873135871E-01, ...
       0.12913281131480700366519063232496E-01, ...
       0.47063366701120542614104893201110E-02, ...
       0.90017947131145907457116462987061E-02, ...
       0.99277995286896342401230585292059E-03, ...
       0.42975457006126745578855094841825E-02 ];

  return
end
