20 January 2020 01:01:47 PM
LAGUERRE_RULE
C version
Compute a Gauss-Laguerre quadrature rule for
Integral ( A <= x < +oo ) exp(-b*(x-a)) f(x) dx
of order N.
The user specifies N, A, B, and FILENAME.
N is the number of points;
A is the left endpoint, usually 0.0:
B is a scale factor, usually 1.0;
FILENAME is used to generate 3 files:
* filename_w.txt - the weight file
* filename_x.txt - the abscissa file.
* filename_r.txt - the region file.
N = 4
A = 0
B = 1
FILENAME = "lag_o4"
Creating quadrature files.
Root file name is "lag_o4"
Weight file will be "lag_o4_w.txt"
Abscissa file will be "lag_o4_x.txt"
Region file will be "lag_o4_r.txt"
LAGUERRE_RULE:
Normal end of execution.
20 January 2020 01:01:47 PM