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: Interrupt Service Routines in Mathematica v.5.2?
  • Next by Date: Re: Problems...
  • Previous by thread: Re: Interrupt Service Routines in Mathematica v.5.2?
  • Next by thread: NIntegrate - Gaussian quadrature more exact than thought