function [ quad_w, quad_x ] = quad ( quad_num )

%% QUAD returns the abscissas and weights for gaussian quadrature on [-1,1].
%
%  Modified:
%
%    03 November 2006
%
%  Author:
%
%    Max Gunzburger
%    Teresa Hodge
%
%    MATLAB translation by John Burkardt
%
%  Parameters:
%
%    Input, integer QUAD_NUM, the order of the quadrature rule.
%
%    Output, real QUAD_W(QUAD_NUM), the quadrature weights.
%
%    Output, real QUAD_X(QUAD_NUM), the quadrature abscissas.
%

%
%  Quadrature points on [-1,1]
%
  quad_x(1) = -0.973906528517172;
  quad_x(2) = -0.865063366688985;
  quad_x(3) = -0.679409568299024;
  quad_x(4) = -0.433395394129247;
  quad_x(5) = -0.148874338981631;
  quad_x(6) =  0.148874338981631;
  quad_x(7) =  0.433395394129247;
  quad_x(8) =  0.679409568299024;
  quad_x(9) =  0.865063366688985;
  quad_x(10) = 0.973906528517172;
%
%  Weight factors
%
  quad_w(1) =  0.066671344308688;
  quad_w(2) =  0.149451349150581;
  quad_w(3) =  0.219086362515982;
  quad_w(4) =  0.269266719309996;
  quad_w(5) =  0.295524224714753;
  quad_w(6) =  0.295524224714753;
  quad_w(7) =  0.269266719309996;
  quad_w(8) =  0.219086362515982;
  quad_w(9) =  0.149451349150581;
  quad_w(10) = 0.066671344308688;