MathGroup Archive 2001

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

Search the Archive

Re: A tough Integral

bobbym1953 wrote:
> Does anyone know how I can get the following Integral to at least 60 places,
> using Mathematica?
> Int(1/(cos(x)+x^2)) between x=0 and x=infinity. Both Integrate and NIntegrate
> seemed helpless.
> Thanks,
> Angela
> bobbym1953 at

60 places?! I wish you luck.

In[16]:= InputForm[nii = NIntegrate[1/(Cos[x]+x^2),  {x,0,Infinity},
WorkingPrecision->200, PrecisionGoal->70, MaxRecursion->12]]

   Numerical integration converging too slowly; suspect one of the
    singularity, value of the integration being 0, oscillatory
integrand, or
    insufficient WorkingPrecision. If your integrand is oscillatory try
    the option Method->Oscillatory in NIntegrate.

   NIntegrate failed to converge to prescribed accuracy after 13
     recursive bisections in x near x = 18538.4.

Out[16]//InputForm= 1.838017695010507757407163`14.843992684648171

There might be ways to improve on this, but I doubt you'll come close to
the desired precision (the effort tends to increase exponentially in the
number of digits). One possibility might be to show that the integral
from some reasonably small t to infinity differs by no more than
10^-(60) from the integral of 1/x^2 from t to infinity. Then put all
your resources into numerically integrating from 0 to t. But offhand I
do not know if a suitable value for t exists, or, if so, how to find it
(best I can readily get is 10^20, which is too big to be terribly

Daniel Lichtblau
Wolfram Research

  • Prev by Date: Re: complexity of Mathematica NullSpace[]
  • Next by Date: RE: Regressions and the Mathematica buttons
  • Previous by thread: Re: A tough Integral
  • Next by thread: Re: A tough Integral