[Date Index]
[Thread Index]
[Author Index]
Re: (revision) NIntegrate that upper limit is infinite
*To*: mathgroup at smc.vnet.net
*Subject*: [mg71963] Re: (revision) NIntegrate that upper limit is infinite
*From*: "dimitris" <dimmechan at yahoo.com>
*Date*: Wed, 6 Dec 2006 06:03:43 -0500 (EST)
*References*: <el12c7$5j5$1@smc.vnet.net>
I think Maxim and Daniel cover completely the subject, so
you are sure that your integral diverges.
Here I provide another short of evidence.
G[a_] := (-(a - Cos[0.5])^(-3/2)/(3*10^2))*Sin[0.5]^2 + (1/(1 -
0.1*0.1))*
((0.005*(1 + Cos[0.5]) + 0.001/1000)*(1 +
Cos[0.5])^(1/2)*(ArcTan[(a - Cos[0.5])^(1/2)/(1 + Cos[0.5])^(1/2)] -
Pi/2) +
(0.005*(1 - Cos[0.5]) - 0.001/1000)*(1 -
Cos[0.5])^(1/2)*ArcTanh[(1 - Cos[0.5])^(1/2)/(a - Cos[0.5])^(1/2)])
P[a_] := N[LegendreP[-2^(-1) + I, a]]
First see how the setting of MaxRecursion affects the estimated value
for the integral.
It seems that there is not a value that the integral converges (even
with setting so small the required PrecisionGoal)
Block[{Message}, ({#,Re[NIntegrate[G[a]*P[a], {a, 1, Infinity},
SingularityDepth -> 1000, MaxRecursion -> #1,
PrecisionGoal -> 4]]} & ) /@ Range[6, 34,
2]]//TableForm[#,TableHeadings->{None,{"maxrecursion","estimated
value"}},TableAlignments->Center]&
Also take a look of the following table:
(TableForm[#1, TableHeadings -> {None, {"upper limit", "estimated
value"}}, TableAlignments -> Center] & )[
({#1, Chop[NIntegrate[G[a]*P[a], {a, 1, #1}, SingularityDepth ->
1000, PrecisionGoal -> 4, MaxRecursion -> 24]]} & ) /@
{1000, 10000, 20000, 50000, 100000, 200000, 300000, 500000, 700000,
1000000, 2000000, 2500000, 3000000}]
Best Regards
Dimitris
Evanescence wrote:
> Dear all:
> I am very sorry for my article contain a lot of mistake before.
> So I check it and state my questions again.
> I am really sorry that cause a lot of confuse for everybody.
> My questions is as follows:
> First I definite a function that is:
> G[a_] := ((-(a -Cos[0.5])^(-3/2))/(3*(10)^2))*(Sin[0.5])^(2) +
> (1/(1-0.1*0.1))*((0.005*(1 + Cos[0.5]) + 0.001/1000)*(1 +
> Cos[0.5])^(1/2)*(ArcTan[(a - Cos[0.5])^(1/2)/(1 + Cos[0.5])^(1/2)] -
> Pi/2) + (0.005*(1 - Cos[0.5]) - 0.001/1000)*(1 -
> Cos[0.5])^(1/2)*(ArcTanh[(1 - Cos[0.5])^(1/2)/(a - Cos[0.5])^(1/2)]))
>
> another function is:
> P[a_]:=Re[N[LegendreP[(-1/2)+I,a]]] where I=(-1)^(1/2)
>
> then
> NIntegrate[G[a]*P[a], {a, 1, Infinity}]
>
> but get error message is:
> NIntegrate's singularity handling has failed at point
> {a}={2.7013362243395366*(10^150)}for the specified
> precision goal.
> Try using larger values for any of $MaxExtraPrecision
> or the options WorkingPrecision, or SingularityDepth and MaxRecursion.
>
> NIntegrate::"inum":
> "Integrand G[a]*P[a] is not numerical at {a}=
> {2.7013362243395366`*(10^150)}.
>
> Please solve my confuse ,thank you very much!
>
> By the way
> G[Infinity]*P[Infinity]=0 ¡A G[1] is singular point¡AP[1]=1
> and I try Plot the graph for G[a]¡AP[a]¡Aand G[a]*P[a] ¡Afrom the
> graph of G[a]*P[a]
> it should convergence for a-->Infinity.
Prev by Date:
**RE: Re: Run option**
Next by Date:
**Re: Points sampled by N[Derivative[]]**
Previous by thread:
**Re: (revision) NIntegrate that upper limit is infinite**
Next by thread:
**Re: NETLink:**
| |