MathGroup Archive 2001

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

Search the Archive

Re: Re: Strange error message from NIntegrate


Use NumberForm. Example:

In[1]:=
Table[NumberForm[Exp[-3.], j], {j, 2,
 15}]
Out[1]=
{"0.05", "0.0498", "0.04979", "0.049787", "0.0497871", "0.04978707", \
"0.049787068", "0.0497870684", "0.04978706837", "0.049787068368", \
"0.0497870683679", "0.04978706836786", "0.049787068367864", \
"0.0497870683678639"}

Tomas Garza
Mexico City
----- Original Message -----
From: "Henrik Holm" <h.holm at spray.no>
To: mathgroup at smc.vnet.net
Subject: [mg28204] [mg28178] Re: Strange error message from NIntegrate


> [BobHanlon at aol.com]
>
> >  You made the power of x a list (brackets vice parenthreses).
>
> Thank you, that seems to be the problem.
>
> >  There are still problems with the numerical integration.  However, you
can
> >
> >      - approximate the answer with a finite upper limit of integration
> >      - change the number of intervals (GaussPoints)
> >      - change the method of numerical integration
> >
> >  NIntegrate[x^(3/2)*Exp[-x]*BesselI[3/2,x/2], {x,0, 1000}]
> >
> >  1.0030037040845279
> >
> >  NIntegrate[x^(3/2)*Exp[-x]*BesselI[3/2,x/2], {x,0, Infinity},
> >    GaussPoints -> 10]
> >
> >  1.0030037040849122
> >
> >  NIntegrate[x^(3/2)*Exp[-x]*BesselI[3/2,x/2], {x,0, Infinity},
> >    Method -> Trapezoidal]
> >
> >  1.0030037040847963
>
> Also. thanks for these tips.  How can I ask Mathematica to
> show the answer with more than a few digits?  (It seems that 6
> digits is standard.)
>
>
> Henrik Holm
>



  • Prev by Date: Serious bug in MatchQ causes crashes and incorrect results
  • Next by Date: Re: Text and 3D graphics overlay
  • Previous by thread: Re: Strange error message from NIntegrate
  • Next by thread: Re: Re: Strange error message from NIntegrate