Re: Singularity-handling transformation employed by NIntegrate
- To: mathgroup at smc.vnet.net
- Subject: [mg71214] Re: Singularity-handling transformation employed by NIntegrate
- From: "antononcube" <antononcube at gmail.com>
- Date: Fri, 10 Nov 2006 06:38:29 -0500 (EST)
- References: <eiut46$h63$1@smc.vnet.net>
NIntegrate uses the so called IMT Rule (the Iri, Moriguti, Takasawa
rule) that changes the integration variable with this transformation
over the interval [xmin, xmax]:
Clear[\[Phi]]
\[Phi][{xmin_, xmax_}, t_] := xmin + (xmax - xmin)*E^(1 - 1/t);
\[Phi][t_] := \[Phi][{0, 1}, t]
It is instructive to look at the plot of the transformation:
Plot[\[Phi][t], {t, 0, 1}, PlotRange -> All, AxesOrigin -> {-0.05,
-0.05}]
and at the plot of its derivative (Jacobian):
Plot[Evaluate[D[\[Phi][t], t]], {t, 0, 1}, PlotRange -> All]
The IMT rule "flattens" the singularity, but the points are clustered
closer to the singular point. That is why in many(most) cases extra
precision is required for the evaluation of the integrand.
Anton Antonov,
Wolfram Research, Inc.
Andrew Moylan wrote:
> Hi all,
>
> Further to a previous question: Can anyone tell me what transformation
> NIntegrate employs when a singularity is detected (i.e., when the
> number of recursive subdivisions specified by SingularityDepth is
> reached)?
>
> Cheers,
> Andrew