MathGroup Archive 2006

[Date Index] [Thread Index] [Author Index]

Search the Archive

Re: singular point list (NIntegrate)

  • To: mathgroup at smc.vnet.net
  • Subject: [mg72214] Re: singular point list (NIntegrate)
  • From: "dimitris" <dimmechan at yahoo.com>
  • Date: Thu, 14 Dec 2006 05:49:45 -0500 (EST)
  • References: <elot8j$p9a$1@smc.vnet.net>

Something like?

nintSing1[f_, {x_, a_, b_, c_}, opts___] := Plus @@ (NIntegrate[f, {x,
#1[[1]], #1[[2]]}, opts] & ) /@  Partition[Range[a, b, c], 2, 1]

nintSing1[x^2*Sin[1/x], {x, -1, 3, 1/5}, MinRecursion -> 6,
MaxRecursion -> 12]
3.820553120346208

Or?

nintSing2[f_, {x_, a_, b_, c___}, opts___] := Plus @@ (NIntegrate[f,
{x, #1[[1]], #1[[2]]}, opts] & ) /@ Partition[{a, c, b}, 2, 1]

nintSing2[x^2*Sin[1/x], {x, -1, 3, 0, 1, 2}]
3.820553120346208

Dimitris

Ï/Ç wtplasar at ehu.es Ýãñáøå:
> Hi,
>
> I want to tell Mathematica to deal with the singular points in a
> numerical integration. In particular, I would like to give easily a
> large number of intermediate points. I would do it like that in a
> trivial example
>
>
> mybounds[xmin_,xmax_,divs_]:=
>     Prepend[Table[xmin+i*(xmax-xmin)/divs,{i,0,divs}],x]
>
>
> NIntegrate[x^2 Sin[1/x],Evaluate[mybounds[-1,3,20]]]
>
> but I think my way of doing it is not very time efficient because of
> the Evaluate command, so I wonder if there is a way round it.
> 
> Thanks,
> 
> Ruth


  • Prev by Date: RE: Mathematica video output without printing frames
  • Next by Date: RE: Arrange cells side by side
  • Previous by thread: Re: circular infinity
  • Next by thread: Re: Complex Numbers to list of points