Groups | Search | Server Info | Keyboard shortcuts | Login | Register [http] [https] [nntp] [nntps]

Groups > comp.soft-sys.math.mathematica > #2280

Re: Numerical/analytical integration around a pole

Show all headers | View raw

On May 10, 7:31 am, "A. Lapraitis" <ffcita...@gmail.com> wrote:
> Hello,
>
> I am trying to integrate a 1/(polynomial) function around its pole.
> Could someone explain, why Integrate and NIntegrate give different
> results?
>
> Thanks a lot in advance!
> A. L.
>
> In[83]:= y[x_] := .24583 + 0.0175121 x - 0.793943 x^2 + 0.534908 x^3 -
> 0.650198 x^4;
>
> (* Locating the pole *)
> In[86]:= FindRoot[y[x] == 0, {x, .6}]
> Out[86]= {x -> 0.601945}
>
> In[84]:= Integrate[1/(y[x] + I/1000), {x, .5, .7}]
> Out[84]= 5.25406 - 15.0137 I
>
> In[85]:= NIntegrate[1/(y[x] + I/1000), {x, .5, .7}]
> Out[85]= 0.323456 - 3.37701 I



I get the same result for:

Integrate[1/(y[x] + I/1000), {x, .5, .7}]

and for

NIntegrate[1/(y[x] + I/1000), {x, .5, .7}]

In both cases I get:
0.323456 - 3.37701*I


Gary Wardall

Back to comp.soft-sys.math.mathematica | Previous | Next — Previous in thread | Find similar | Unroll thread

Thread

Numerical/analytical integration around a pole "A. Lapraitis" <ffcitatos@gmail.com> - 2011-05-10 12:31 +0000
  Re: Numerical/analytical integration around a pole Gary Wardall <gwardall@gmail.com> - 2011-05-11 08:26 +0000

csiph-web