Re: Parametric Numerical Integral
- To: mathgroup at smc.vnet.net
- Subject: [mg63120] Re: Parametric Numerical Integral
- From: "antononcube" <antononcube at gmail.com>
- Date: Thu, 15 Dec 2005 03:06:30 -0500 (EST)
- References: <dnopgm$2hb$1@smc.vnet.net>
- Sender: owner-wri-mathgroup at wolfram.com
If you use a definition that specifies the argument as a number
everything works fine:
In[1]:=
Clear[int];
int[(p_)?NumericQ] := NIntegrate[If[t == 0, 0, 1/(E^((1/2)*(-2 +
1/t)^2)*(Sqrt[2*Pi]*t))],
{t, -p, p}, MaxRecursion -> 50, WorkingPrecision -> 50]
In[3]:= << "NumericalMath`NLimit`"
NLimit[int[p], p -> Infinity]
In[4]:=
NIntegrate::slwcon:
Numerical integration converging too slowly; suspect one of the following:
singularity, value of the integration being 0, oscillatory integrand, or
insufficient WorkingPrecision. If your integrand is oscillatory try using
the option Method->Oscillatory in NIntegrate.
NIntegrate::slwcon:
Numerical integration converging too slowly; suspect one of the following:
singularity, value of the integration being 0, oscillatory integrand, or
insufficient WorkingPrecision. If your integrand is oscillatory try using
the option Method->Oscillatory in NIntegrate.
NIntegrate::slwcon:
Numerical integration converging too slowly; suspect one of the following:
singularity, value of the integration being 0, oscillatory integrand, or
insufficient WorkingPrecision. If your integrand is oscillatory try using
the option Method->Oscillatory in NIntegrate.
General::stop: Further output of NIntegrate::slwcon
will be suppressed during this calculation.
Out[4]= 0.639988
Anton Antonov,
Wolfram Research, Inc.