August 30 2007 3:12:27.118 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.00000 ! w( 1) = 0.6031541043416342 w( 2) = 0.3574186924377997 w( 3) = 0.3888790851500539E-01 w( 4) = 0.5392947055613276E-03 x( 1) = 0.3225476896193923 x( 2) = 1.745761101158346 x( 3) = 4.536620296921128 x( 4) = 9.395070912301133 r( 1) = 0.000000000000000 r( 2) = 0.1797693134862000+309 GEN_LAGUERRE_RULE: Normal end of execution. August 30 2007 3:12:27.118 PM