August 30 2007   3:18:08.021 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
%  ORDER =        8
%  A     =    0.00000    
%  ALPHA =    1.50000    
%
    w( 1) =   0.1900505753913614    ;
    w( 2) =   0.5604380973781597    ;
    w( 3) =   0.4360070168210791    ;
    w( 4) =   0.1274485793102171    ;
    w( 5) =   0.1476947042302810E-01;
    w( 6) =   0.6195794633122050E-03;
    w( 7) =   0.7164389487568510E-05;
    w( 8) =   0.1070030801415182E-07;
 
    x( 1) =   0.5487420203051165    ;
    x( 2) =    1.638181097086517    ;
    x( 3) =    3.315036043679075    ;
    x( 4) =    5.640313098679649    ;
    x( 5) =    8.715661629230398    ;
    x( 6) =    12.71773901227562    ;
    x( 7) =    17.99242472404181    ;
    x( 8) =    25.43190237470182    ;
 
    r(1) =    0.000000000000000    ;
    r(2) = Inf;
 
GEN_LAGUERRE_RULE:
  Normal end of execution.
 
August 30 2007   3:18:08.021 PM