February 26 2008   7:41:58.204 PM
 
CHEBYSHEV2_RULE
  FORTRAN90 version
 
  Compute a Gauss-Chebyshev type 1 rule for
 
    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 = "F77".
c
c  Weights W, abscissas X and range R
c  for a Gauss-Chebyshev type 2 quadrature rule
c  ORDER =        4
c
c  Standard rule:
c    Integral ( -1 <= x <= +1 ) f(x)  sqrt ( 1 - x^2 ) dx
c    is to be approximated by
c    sum ( 1 <= I <= ORDER ) w(i) * f(x(i)).
c
      w( 1) =   0.2170787134227061    
      w( 2) =   0.5683194499747424    
      w( 3) =   0.5683194499747423    
      w( 4) =   0.2170787134227060    
 
      x( 1) =  -0.8090169943749473    
      x( 2) =  -0.3090169943749473    
      x( 3) =   0.3090169943749475    
      x( 4) =   0.8090169943749475    
 
      r( 1) =  -1.0000000000000000    
      r( 2) =   1.0000000000000000    
 
CHEBYSHEV2_RULE:
  Normal end of execution.
 
February 26 2008   7:41:58.209 PM