August 30 2007   3:12:15.281 PM
 
GEN_LAGUERRE_RULE
  FORTRAN90 version
 
  Compute a generalized Gauss-Laguerre rule for approximating
 
    Integral ( 0 <= x < oo ) x^ALPHA exp(-x) g(x) dx
 
  of order ORDER and parameter ALPHA.
 
  The user specifies ORDER, ALPHA, 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.
c
c  Weights W, abscissas X and range R
c  for a generalized Gauss-Laguerre quadrature rule
c  of ORDER =        4
c  with ALPHA =  -0.500000    
c
      w( 1) =    1.322294130254031    
      w( 2) =   0.4156046846751229    
      w( 3) =   0.3415596873061822E-01
      w( 4) =   0.3992081761643634E-03
 
      x( 1) =   0.1453035215033171    
      x( 2) =    1.339097288126361    
      x( 3) =    3.926963501358287    
      x( 4) =    8.588635689012035    
 
      r( 1) =    0.000000000000000    
      r( 2) =   0.1797693134862000+309
 
GEN_LAGUERRE_RULE:
  Normal end of execution.
 
August 30 2007   3:12:15.282 PM