Re: Integrating over area of intersection
- To: mathgroup at smc.vnet.net
- Subject: [mg44191] Re: Integrating over area of intersection
- From: koopman at sfu.ca (Ray Koopman)
- Date: Sat, 25 Oct 2003 06:27:02 -0400 (EDT)
- References: <bn8esp$nle$1@smc.vnet.net> <bnap4l$4jv$1@smc.vnet.net>
- Sender: owner-wri-mathgroup at wolfram.com
"Toni Danza" <nospam at yoohoo.com> wrote in message
news:<bnap4l$4jv$1 at smc.vnet.net>...
> Adding:
> It would be really useful NOT to use the
> Nintegrate[If[x^2 + y^2 < fJ^2 && ...,1,0],{x,-inf,inf},{y,-inf,inf}]
> form, because it takes just waaaayyyy too long
>
> Ideally, I would like to extract the y boundary for every x to use as
> integral limits , that should speed it up...
If you're willing to do it numerically then this should get it:
With[{xmin = Max[-fJ,f1-fH,f2-fH], xmax = Min[fJ,f1+fH,f2+fH]},
If[xmin >= xmax, 0, NIntegrate[<whatever>, {x,xmin,xmax},
{y,-Sqrt[Min[fJ^2-x^2, fH^2-(x-f1)^2, fH^2-(x-f2)^2]],
Sqrt[Min[fJ^2-x^2, fH^2-(x-f1)^2, fH^2-(x-f2)^2]]}]]]