Re: Integrating over area of intersection
- To: mathgroup at smc.vnet.net
- Subject: [mg44187] Re: Integrating over area of intersection
- From: "Toni Danza" <nospam at yoohoo.com>
- Date: Sat, 25 Oct 2003 06:26:47 -0400 (EDT)
- References: <bn8esp$nle$1@smc.vnet.net> <bnaodq$4g8$1@smc.vnet.net>
- Sender: owner-wri-mathgroup at wolfram.com
Yes, I have tried that, but it takes literally over an hour to finish... In the process it keeps complaining about lack of convergence. NIntegrate[If[...],...] is attempting to integrate a discontinuous function, causing mayhem in convergence algorithms.... OTOH, the integration limits are readily available, all I need to do is figure out how to plug the result of Reduce[region,{x,y}] into the limits of NIntegrate. The function will be nice and continuous here, should be no problem Any ideas? "Steve Luttrell" <luttrell at _removemefirst_westmal.demon.co.uk> wrote in message news:bnaodq$4g8$1 at smc.vnet.net... > Here is an example of the sort of thing you can do. > > First of all read in the Calculus`Integration` package which gives you the > Boole function for doing integrals over regions defined by inequalities: > > << Calculus`Integration` > > Now integrate the function x^2*y^2 (as an example) over the region of > interest (as an example fJ = 3, fH = 2, f1 = 1/2, f2 = -4^(-1)}): > > With[{fJ = 3, fH = 2, f1 = 1/2, f2 = -4^(-1)}, > Integrate[Boole[x^2 + y^2 < fJ^2 && > (x - f1)^2 + y^2 < fH^2 && (x - f2)^2 + y^2 < > fH^2]*x^2*y^2, {x, -Infinity, Infinity}, > {y, -Infinity, Infinity}]] > > which gives the result: > > -((56359*Sqrt[247])/163840) + (11/3)*ArcCos[3/16] + > (35/6)*ArcSin[Sqrt[13/2]/4] > > -- > Steve Luttrell > West Malvern, UK > > "Toni Danza" <nospam at yoohoo.com> wrote in message > news:bn8esp$nle$1 at smc.vnet.net... > > OK, I have three functions that are defined within their respective > circles. > > I would like to integrate over the intersection of the three circles. > > > > Here's what I have done: > > define region of integration: > > region = x^2 + y^2 < fJ^2 && (x - f1)^2 + y^2 < fH^2 && (x - f2)^2 + > y^2 > > < fH^2 > > > > Then I try to solve for the intersection using > > > > Reduce[region,{x,y}] > > > > and the result is something like (only works with numerical parameters...) > > > > -0.4<x<0.3 && sqrt(....)< y <sqrt(...) || -0.3<x<-0.2 && sqrt(....)< y > > <sqrt(...) > > > > How do I use this result to do integration over the region? > > > >