Re: Nested numerical integration
- To: mathgroup at smc.vnet.net
- Subject: [mg99205] Re: Nested numerical integration
- From: "Sjoerd C. de Vries" <sjoerd.c.devries at gmail.com>
- Date: Wed, 29 Apr 2009 06:42:01 -0400 (EDT)
- References: <gt6flm$rtj$1@smc.vnet.net>
The answer seems to be correct, at least, it is the same as when
calculated with the synbolic Integrate instead of Nintegrate. You can
remove the error in the following way:
In[266]:= ni[w_?NumberQ] := NIntegrate[1/(s - w), {s, 1, 5}]
In[268]:= NIntegrate[w^2*ni[w], {w, -5, -1}]
Out[268]= 25.8364
In[270]:= Integrate[w^2*Integrate[1/(s - w), {s, 1, 5}], {w, -5, -1}]
Out[270]= -16 - 84 Log[3] + (250 Log[5])/3
In[250]:= -16 - 84 Log[3] + (250 Log[5])/3 // N
Out[250]= 25.8364
Cheers -- Sjoerd
On Apr 28, 10:44 am, tsg.mo... at googlemail.com wrote:
> Hi, I'd like to compute a nested integral like,
>
> NIntegrate[ w^2 * NIntegrate[1/(s-w), {s, 1, 5}], {w, -5, -1}]
>
> Unfortunately, mathematica gives me an error (NIntegrate::inumr) and
> it also outputs the wrong value. Is there a way to do these types of
> integrals?