September  2 2007   9:48:33.554 AM
 
JACOBI_RULE
  FORTRAN90 version
 
  Compute a Gauss-Jacobi rule for approximating
 
    Integral ( -1 <= x <= +1 ) (1-x)^ALPHA (1+x)^BETA f(x) dx
 
  of order ORDER.
 
  The user specifies ORDER, ALPHA, BETA, 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.
 
  Input summary:
 
  ORDER =         4
  ALPHA =     0.00000    
  BETA =     0.00000    
  OUTPUT = "C++".
//
//  Weights W, abscissas X and range R
//  for a Gauss-Jacobi quadrature rule
//  ORDER =        4
//  ALPHA =    0.00000    
//  BETA  =    0.00000    
//
//  Standard rule:
//    Integral ( -1 <= x <= +1 ) (1-x)^ALPHA (1+x)^BETA f(x) dx
//    is to be approximated by
//    sum ( 1 <= I <= ORDER ) w(i) * f(x(i)).
//
  w[ 0] =   0.3478548451374422    ;
  w[ 1] =   0.6521451548625246    ;
  w[ 2] =   0.6521451548625246    ;
  w[ 3] =   0.3478548451374422    ;
 
  x[ 0] =  -0.8611363115940526    ;
  x[ 1] =  -0.3399810435848563    ;
  x[ 2] =   0.3399810435848563    ;
  x[ 3] =   0.8611363115940526    ;
 
  r[ 0] =  -1.0000000000000000    ;
  r[ 1] =   1.0000000000000000    ;
 
JACOBI_RULE:
  Normal end of execution.
 
September  2 2007   9:48:33.554 AM