19 March 2020 02:29:13 PM

gen_laguere_exactness
  C++ version

  Compiled on Mar 19 2020 at 14:21:43.

  Investigate the polynomial exactness of a generalized Gauss-Laguerre
  quadrature rule by integrating exponentially weighted
  monomials up to a given degree over the [0,+oo) interval.

  The rule may be defined on another interval, [A,+oo)
  in which case it is adjusted to the [0,+oo) interval.

  The quadrature file rootname is "gen_lag_o4_a0.5".

  The requested maximum monomial degree is = 10

  The requested value of ALPHA = 0.5

gen_laguere_exactness: User input:
  Quadrature rule X file = "gen_lag_o4_a0.5_x.txt".
  Quadrature rule W file = "gen_lag_o4_a0.5_w.txt".
  Quadrature rule R file = "gen_lag_o4_a0.5_r.txt".
  Maximum degree to check = 10
  Weighting exponent ALPHA = 0.5
  OPTION = 0, integrate x^alpha*exp(-x)*f(x)

  Spatial dimension = 1
  Number of points  = 4

  The quadrature rule to be tested is
  a generalized Gauss-Laguerre rule
  ORDER = 4
  with A =     0
  and  ALPHA = 0.5

  Standard rule:
    Integral ( A <= x < +oo ) x^alpha exp(-x) f(x) dx
    is to be approximated by
    sum ( 1 <= I <= ORDER ) w(i) * f(x(i)).

  Weights W:

  w[ 0] =       0.4530087465586076
  w[ 1] =       0.3816169601717996
  w[ 2] =      0.05079462757224078
  w[ 3] =    0.0008065911501100311

  Abscissas X:

  x[ 0] =       0.5235260767382691
  x[ 1] =        2.156648763269094
  x[ 2] =        5.137387546176711
  x[ 3] =        10.18243761381592

  Region R:

  r[ 0] =                        0
  r[ 1] =                    1e+30

  A generalized Gauss-Laguerre rule would be able to exactly
  integrate monomials up to and including degree = 7

          Error          Degree

                         0   0
      3.34067341817812e-16   1
     4.008808101813744e-16   2
     4.581494973501421e-16   3
     5.429919968594276e-16   4
     9.872581761080503e-16   5
     1.579613081772881e-15   6
     2.203357734575505e-15   7
       0.01053064582476666   8
       0.05043625105545014   9
        0.1330978618904358  10

gen_laguere_exactness:
  Normal end of execution.

19 March 2020 02:29:13 PM