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

%*****************************************************************************80
%
%% rule16() returns the rule of degree 16.
%
%  Discussion:
%
%    Order 16 (55 pts)
%    1/6 data for 16-th order quadrature with 13 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.16279607394216944053888791563128, ...
       -0.36133516018585422888297687826246, ...
        0.00000000000000000000000000000000, ...
       -0.80996884917848457646874211083899, ...
       -0.30011609466308557056235128866156, ...
        0.00000000000000000000000000000000, ...
        0.00000000000000000000000000000000, ...
       -0.60911897435540735910318766737015, ...
        0.00000000000000000000000000000000, ...
       -0.56820549744094988938273289446974, ...
        0.00000000000000000000000000000000, ...
        0.00000000000000000000000000000000, ...
        0.00000000000000000000000000000000 ];
  y = [ ... ...
       -0.56061007710964974635124810013984, ...
       -0.52485883552713864174217026883955, ...
        0.92373339140338472305543117791951, ...
       -0.55862165096055508852660603004778, ...
       -0.39262157635321675549452096003032, ...
        0.31873771449384928263441293513356, ...
       -0.27527401254502863816148429232788, ...
       -0.48840198582588110096280723102788, ...
        0.63486450609449659829895244998069, ...
       -0.57094697658781472615514574003317, ...
       -0.47186155418111734634796824747685, ...
        0.00000000000000000000000000000000, ...
        0.10956666024303233991933608013003E+01 ];
  w = [ ... ...
        0.10768394677601291424698896830942E-01, ...
        0.18403461944093881198776232104775E-01, ...
        0.81763898630962949174826141110659E-02, ...
        0.75698865940531301551762664596891E-02, ...
        0.41648559391949793671705497882676E-01, ...
        0.27100623100161070449309097892289E-01, ...
        0.26969791338286996589221138520990E-01, ...
        0.23232761214636807134379750252095E-01, ...
        0.18940606006195797095777111782639E-01, ...
        0.60732744287626264990846070830020E-02, ...
        0.17828637150989055090637459591942E-01, ...
        0.10139891906267619230864844708608E-01, ...
        0.24933912093210466794229664121214E-02 ];

  return
end
