01 March 2008 10:49:54 AM

CHEBYSHEV1_RULE
  C++ version

  Compiled on Mar  1 2008 at 10:49:35.

  Compute a Gauss-Chebyshev type 1 rule for approximating

    Integral ( -1 <= x <= +1 ) f(x) / sqrt ( 1 - x^2 ) 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 = 4

  OUTPUT option is "C++".
//
//  Weights W, abscissas X and range R
//  for a Gauss-Chebyshev type 1 quadrature rule.
//  ORDER = 4
//
//  Standard rule:
//    Integral ( -1 <= x <= +1 ) f(x) / sqrt ( 1 - x^2 ) dx
//    is to be approximated by
//    sum ( 1 <= I <= ORDER ) w(i) * f(x(i)).
//
  w[0] = 0.7853981633974483;
  w[1] = 0.7853981633974483;
  w[2] = 0.7853981633974483;
  w[3] = 0.7853981633974483;

  x[0] = -0.9238795325112867;
  x[1] = -0.3826834323650897;
  x[2] = 0.3826834323650898;
  x[3] = 0.9238795325112867;

  r[0] = -1;
  r[1] = 1;

CHEBYSHEV1_RULE:
  Normal end of execution.

01 March 2008 10:49:54 AM