27 June 2009 09:39:13 PM

CLENSHAW_CURTIS_RULE
  C++ version

  Compiled on Jun 27 2009 at 21:35:17.

  Compute a Clenshaw Curtis rule for approximating

    Integral ( -1 <= x <= +1 ) f(x) dx

  of order ORDER.

  The user specifies ORDER and OUTPUT.

  OUTPUT is:

  "C++" for printed C++ output;
  "F77" for printed Fortran77 output;
  "F90" for printed Fortran90 output;
  "MAT" for printed MATLAB output;

  or:

  "filename" to generate 3 files:

    filename_w.txt - the weight file
    filename_x.txt - the abscissa file.
    filename_r.txt - the region file.

  The requested order of the rule is = 8

  OUTPUT option is "MAT".
%
%  Weights W, abscissas X and range R
%  for a Clenshaw Curtis quadrature rule
%  ORDER = 8
%
%  Standard rule:
%    Integral ( -1 <= x <= +1 ) f(x) dx
%    is to be approximated by
%    sum ( 1 <= I <= ORDER ) w(i) * f(x(i)).
%
  w(1) = 0.02040816326530613;
  w(2) = 0.1901410072182084;
  w(3) = 0.3522424237181591;
  w(4) = 0.4372084057983264;
  w(5) = 0.4372084057983264;
  w(6) = 0.3522424237181591;
  w(7) = 0.1901410072182084;
  w(8) = 0.02040816326530613;

  x(1) = -1;
  x(2) = -0.900968867902419;
  x(3) = -0.6234898018587335;
  x(4) = -0.2225209339563143;
  x(5) = 0.2225209339563145;
  x(6) = 0.6234898018587336;
  x(7) = 0.9009688679024191;
  x(8) = 1;

  r(1) = -1;
  r(2) = 1;

CLENSHAW_CURTIS_RULE:
  Normal end of execution.

27 June 2009 09:39:13 PM