       Re: Area of inside contour of continuous function?

• To: mathgroup at smc.vnet.net
• Subject: [mg84238] Re: Area of inside contour of continuous function?
• From: "Steve Luttrell" <steve at _removemefirst_luttrell.org.uk>
• Date: Fri, 14 Dec 2007 23:17:24 -0500 (EST)
• References: <fjun8o\$24v\$1@smc.vnet.net>

```This does the sort of thing that you want:

In:= f[x_,y_,a_,b_]:=1-a x^2-b y^2;

In:= Integrate[f[x,y,a,b]Boole[f[x,y,a,b]>k
f[0,0,a,b]],{x,-\[Infinity],\[Infinity]},{y,-\[Infinity],\[Infinity]},Assumptions->a>0&&b>0&&0<k<1]

Out= -(((-1+k^2) \[Pi])/(2 Sqrt[a b]))

Steve Luttrell
West Malvern, UK

"Gareth Russell" <russell at njit.edu> wrote in message
news:fjun8o\$24v\$1 at smc.vnet.net...
> Hi Group,
>
> I have a two-parameter smooth continuous function f[x_,y_] that is
> unimodal (it's a log-likeihood function), and I know maxf, the height
> of the peak. I would like to calculate the area inside a contour that
> is a given number of units u below the peak. What is the easiest way to
> do this?
>
> My own idea is to construct a separate function with discontinuities
> along the lines of
>
> g[x_,y_]:=If[f[x,y]>(maxf-u),1,0]
>
> And do numerical integration over a region big enough to contain the
> contour (which of course, I can see from a ContourPlot).
>
> But am I missing something much easier? Searching the mathgroup
> archives I came across how to find the area of a discrete polygon, so I
> realize that another method would be to extract the contour data from
> the ContourPlot object and apply that, but it seems like a bit of hack!
>
> Thanks,
>
> Gareth
>
> --
> Gareth Russell
> NJIT
>
>

```

• Prev by Date: ParametricPlot3d not showing up when I use the manipulate command
• Next by Date: Re: Re: How to express this function.
• Previous by thread: Area of inside contour of continuous function?
• Next by thread: alignment in TabView (v6) labels