Re: Piecewise Integration Woes
- To: mathgroup at smc.vnet.net
- Subject: [mg63399] Re: [mg63388] Piecewise Integration Woes
- From: pdesai1 at umbc.edu
- Date: Thu, 29 Dec 2005 02:57:30 -0500 (EST)
- References: <200512280855.DAA28668@smc.vnet.net>
- Sender: owner-wri-mathgroup at wolfram.com
> mathgroup at smc.vnet.net; Tue, 27 Dec 2005 17:33:45 -0600 (CST)
> From: Bruce Colletti <vze269bv at verizon.net>
To: mathgroup at smc.vnet.net
> Subject: [mg63399] [mg63388] Piecewise Integration Woes
>
> Re Mathematica 5.2 under WinXP.
>
> h, g, and (presumably f) are continuous functions.
>
> So, what causes the NIntegrate::slwcon warning? Although continuous,
> g(h(x)) isn't always differentiable. Could this be the cause?
>
> Thankx.
>
> Bruce
>
> ------------
>
>
> T = {{16, 113.32}, {20, 68.11}, {900, 0.15}, {1000, 0.13}};
>
> f = Interpolation@T;
>
> h[r_] := Piecewise[{{T[[1, 2]], r < T[[
> 1, 1]]}, {f@r, T[[1, 1]] <= r < T[[-1, 1]]}}, T[[-1, 2]]];
>
> g[x_] := Piecewise at {{0.99, 29 <= x}, {0.99(x - 14.5)/(29 -
> 14.5), 14.5 <= x < 29}}
>
> NIntegrate[g[h@r], {r, 14, 40}, MaxRecursion -> 30]
>
>
Have you tried PiecewiseExpand
NIntegrate[g[h@r]//PiecewiseExpand, {r, 14, 40}, MaxRecursion -> 30]
it seems to work without spitting any error
Version 5.1.1 for windows
Hope this helps
Pratik