September  1 2007   9:01:49.586 AM
 
LAGUERRE_RULE
  FORTRAN90 version
 
  Compute a Gauss-Laguerre rule for approximating
 
    Integral ( A <= x < oo ) exp(-x) g(x) dx
 
  of order ORDER.
 
  For now, A is fixed at 0.0.
 
  The user specifies ORDER, OPTION, and OUTPUT.
 
  OPTION is:
 
    0 if the integrand is factored (standard rule):
      Integral ( A to oo ) exp(-x) f(x) dx
 
    1 if the integrand is NOT factored (modified rule):
      Integral ( A to oo )         f(x) dx
 
    For OPTION = 1, the weights of the standard rule
    are multiplied by exp(+x).
 
  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
  A =         0.00000      (Not user input)
  OPTION =        0
  OUTPUT = "C++".
//
//  Weights W, abscissas X and range R
//  for a Gauss-Laguerre quadrature rule
//  ORDER =        4
//  A =        0.00000    
//
//  OPTION = 0, standard rule:
//    Integral ( A <= x < oo ) exp(-x) f(x) dx
//    is to be approximated by
//    sum ( 1 <= I <= ORDER ) w(i) * f(x(i)).
//
  w[ 0] =   0.6031541043416342    ;
  w[ 1] =   0.3574186924377997    ;
  w[ 2] =   0.3888790851500539E-01;
  w[ 3] =   0.5392947055613276E-03;
 
  x[ 0] =   0.3225476896193923    ;
  x[ 1] =    1.745761101158346    ;
  x[ 2] =    4.536620296921128    ;
  x[ 3] =    9.395070912301133    ;
 
  r[ 0] =    0.000000000000000    ;
  r[ 1] =   0.1797693134862000+309;
 
LAGUERRE_RULE:
  Normal end of execution.
 
September  1 2007   9:01:49.587 AM