August 30 2007   3:11:59.087 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.
//
//  Weights W, abscissas X and range R
//  for a generalized Gauss-Laguerre quadrature rule
//  of ORDER =        4
//  with ALPHA =   0.500000    
//
  w[ 0] =   0.4530087825780014    ;
  w[ 1] =   0.3816169905147275    ;
  w[ 2] =   0.5079463161099657E-01;
  w[ 3] =   0.8065912142432840E-03;
 
  x[ 0] =   0.5235260767382691    ;
  x[ 1] =    2.156648763269094    ;
  x[ 2] =    5.137387546176711    ;
  x[ 3] =    10.18243761381592    ;
 
  r[ 0] =    0.000000000000000    ;
  r[ 1] =   0.1797693134862000+309;
 
GEN_LAGUERRE_RULE:
  Normal end of execution.
 
August 30 2007   3:11:59.087 PM