Re: Parametric Numerical Integral
- To: mathgroup at smc.vnet.net
- Subject: [mg63131] Re: [mg63095] Parametric Numerical Integral
- From: Bob Hanlon <hanlonr at cox.net>
- Date: Thu, 15 Dec 2005 03:06:39 -0500 (EST)
- Reply-to: hanlonr at cox.net
- Sender: owner-wri-mathgroup at wolfram.com
Off[NIntegrate::slwcon];
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];
int[10]
0.6183797223358392902748797062033600519932
Needs["NumericalMath`NLimit`"];
NLimit[int[p], p->Infinity, WorkingPrecision->50]
0.63998805892006372723475974276826659
Bob Hanlon
>
> From: "Nikola Letic" <nletic at yahoo.com>
To: mathgroup at smc.vnet.net
> Date: 2005/12/14 Wed AM 04:35:57 EST
> Subject: [mg63131] [mg63095] Parametric Numerical Integral
>
> Hello,
> I have the following problem:
> I define a parametric integral as follows:
> int[p_]:=NIntegrate[If[t == 0, 0,1/(E^((1/2)*(-2 +
> 1/t)^2)*(Sqrt[2*Pi]*t))],{t, -p, p}, MaxRecursion -> 50,WorkingPrecision ->
> 50]
>
> it works fine for numerical values of p such as int[10], but when I try
> something like this:
>
> << NumericalMath`NLimit`
>
> NLimit[int[p], p -> Infinity]
>
> it ends up with error messages:
>
> NIntegrate::nlim: t = 1.00000 p is not a valid limit of integration.
>
> Now I tought that NLimit should also deal only with numerical values, but I
> was wrong. Seems that the way I defined this function above has serious
> limitations in its further usage.
>
> Does any of you Mathematica Gurus know an easy workaround for this?
>
> Thanks in advance
>
> Nikola Letic
>
>
>