September  1 2007   9:02:46.227 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 =         8
  A =         0.00000      (Not user input)
  OPTION =        0
  OUTPUT = "MAT".
%
%  Weights W, abscissas X and range R
%  for a Gauss-Laguerre quadrature rule
%  ORDER =        8
%  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( 1) =   0.3691885893416379    ;
    w( 2) =   0.4187867808143433    ;
    w( 3) =   0.1757949866371718    ;
    w( 4) =   0.3334349226121565E-01;
    w( 5) =   0.2794536235225671E-02;
    w( 6) =   0.9076508773358207E-04;
    w( 7) =   0.8485746716272539E-06;
    w( 8) =   0.1048001174871509E-08;
 
    x( 1) =   0.1702796323051010    ;
    x( 2) =   0.9037017767993799    ;
    x( 3) =    2.251086629866131    ;
    x( 4) =    4.266700170287659    ;
    x( 5) =    7.045905402393466    ;
    x( 6) =    10.75851601018100    ;
    x( 7) =    15.74067864127800    ;
    x( 8) =    22.86313173688927    ;
 
    r(1) =    0.000000000000000    ;
    r(2) = Inf;
 
LAGUERRE_RULE:
  Normal end of execution.
 
September  1 2007   9:02:46.228 AM