MathGroup Archive 2000

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

Search the Archive

Re: Problem with the GaussianQuadratureWeights[n, a, b]


Hi,

setting the precision will help and 

GaussianQuadratureWeights[42, 0, 1, 32]

work fine.

Cheng Liu wrote:
> 
> Hi,
> 
>         I accidently run into this strange behavior.  I am using the
> package NumericalMath`GaussianQuadrature` and the function
> GaussianQuadratureWeights[n, a, b] to obtain the abscissas and weights
> for the n-point Gaussian quadrature formula on the interval (a, b).
> The specific interval in my problem is (0, 1),  I noticed that
> when 40 <= n <= 50, if n is odd, the function gives correct answer,
> but when n is even, the function gives the following error message:
> 
> In[43]:=
> GaussianQuadratureWeights[42, 0, 1]
> \!\(Power::"infy" \(\(:\)\(\ \)\)
>      "Infinite expression \!\(1\/0.`\) encountered."\)
> \!\(Power::"infy" \(\(:\)\(\ \)\)
>      "Infinite expression \!\(1\/0.`\) encountered."\)
> \!\(Power::"infy" \(\(:\)\(\ \)\)
>      "Infinite expression \!\(1\/0.`\) encountered."\)
> General::stop :
>      Further output of Power::infy will be suppressed during this calculation.
> f::indet :
>      Indeterminate expression - 1672.96 + ComplexInfinity +
>        ComplexInfinity + \[LeftSkeleton]6\[RightSkeleton] + ComplexInfinity +
>        ComplexInfinity encountered.
SNIPP SNAPP
> Thanks.
> 
> Cheng
> --


  • Prev by Date: Re: Integral involving erfc and j1
  • Next by Date: Re: Question: nonlinear differential equation with boundary conditions
  • Previous by thread: Problem with the GaussianQuadratureWeights[n, a, b]
  • Next by thread: Re: Re: Problem with the GaussianQuadratureWeights[n,a, b]