MathGroup Archive 2007

[Date Index] [Thread Index] [Author Index]

Search the Archive

NIntegrate - Gaussian quadrature more exact than thought


NIntegrate[x^4, {x, 0, 1}, Method->GaussKronrod, GaussPoints->2]

gives 0.2, the exact result, although the Gauss quadrature should be
inexact on a polynomial of degree 2n+2 where n is the number of
GaussPoints. More exactly, we expected the same results as here:

<< NumericalMath`GaussianQuadrature`
gw = GaussianQuadratureWeights[2, 0, 1]
f[{x_, y_}] := x^4 y
Total[f /@ gw]
0.194444

This is inexact, OK.

Why is NIntegrate so good?
Something I may have missed.

Thank you for your help.

Janos



  • Prev by Date: change filename automatically in Export Command
  • Next by Date: Problems with simplify and Sqrt[a^2]
  • Previous by thread: NIntegrate - Gaussian quadrature more exact than thought
  • Next by thread: Re: NIntegrate - Gaussian quadrature more exact than thought