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

• To: mathgroup at smc.vnet.net
• Subject: [mg26444] Re: Problem with the GaussianQuadratureWeights[n, a, b]
• From: Jens-Peer Kuska <kuska at informatik.uni-leipzig.de>
• Date: Thu, 21 Dec 2000 01:51:45 -0500 (EST)
• Organization: Universitaet Leipzig
• References: <91pjn3\$5vp@smc.vnet.net>
• Sender: owner-wri-mathgroup at wolfram.com

```Hi,

setting the precision will help and

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]:=
> \!\(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