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

%*****************************************************************************80
%
%% rule12() returns the rule of degree 12.
%
%  Discussion:
%
%    Order 12 (33 pts)
%    1/6 data for   12-th order quadrature with   8 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.00000000000000000000000000000000, ...
       0.00000000000000000000000000000000, ...
      -0.37279552304503872030654041088416, ...
      -0.72405807527665067283016173568468, ...
      -0.39365448416805093945053087420319, ...
       0.00000000000000000000000000000000, ...
       0.00000000000000000000000000000000 ];
  y = [ ... ...
       0.21432682937950204373207766849052, ...
       0.77622032111803989421823022245158, ...
      -0.36989093457992084889029420576697, ...
      -0.37591965438942697512621254516818, ...
      -0.54031470967359184488712033281188, ...
      -0.53745394007281752834331765730183, ...
      -0.53648686378750904365565925207548, ...
       0.10693230309921692958873012567969E+01 ];
  w = [ ... ...
       0.41154432712824560996161389102816E-01, ...
       0.18744876429730660409284009808805E-01, ...
       0.32848111684339866575373457398258E-01, ...
       0.56890409960601999633893263927781E-01, ...
       0.19851236078200935116627045158968E-01, ...
       0.28668810178252124180808712582749E-01, ...
       0.15968477487762470187363091359494E-01, ...
       0.52193142937027930370255142939621E-02 ];

  return
end
