Goat problem

Short description: Recreational mathematics planar boundary and area problem

The goat problem is either of two related problems in recreational mathematics involving at least figuratively a tethered goat (horse, bull) grazing a circular area: the interior grazing problem and the exterior grazing problem. The former involves grazing the interior of a circular area, and the latter, grazing the exterior of a circular area.

The original problem was the exterior grazing problem and appeared in the 1748 edition of the English annual journal The Ladies' Diary: or, the Woman’s Almanack, designated as Question CCCIII attributed to Upnorensis (an unknown historical figure), stated thus:

Observing a horse tied to feed in a gentlemen’s park, with one end of a rope to his fore foot, and the other end to one of the circular iron rails, inclosing a pond, the circumference of which rails being 160 yards, equal to the length of the rope, what quantity of ground at most, could the horse feed?

The related problem involving area in the interior of a circle without reference to barnyard animals first appeared in 1894 in the first edition of the renown journal American Mathematical Monthly. Attributed to Charles E. Myers, it was stated as:

A circle containing one acre is cut by another whose center is on the circumference of the given circle, and the area common to both is one-half acre. Find the radius of the cutting circle.

The solutions in both cases are non-trivial but yield to straightforward application of trigonometry, analytical geometry or integral calculus. Both problems are intrinsically transcendental – they do not have closed-form analytical solutions in the Euclidean plane. The numerical answers must be obtained by an iterative approximation procedure. The goat problems do not yield any new mathematical insights; rather they are primarily exercises in how to artfully deconstruct problems in order to facilitate solution.

Three-dimensional analogues and planar boundary/area problems on other shapes, including the obvious rectangular barn and/or field, have been proposed and solved. A generalized solution for any smooth convex curve like an ellipse, and even unclosed curves, has been formulated.

Exterior grazing problem

thumb|400px|Goat tethered to silo at v, grazing an area under an involute

The question about the grazable area outside a circle is considered. This concerns a situation where the animal is tethered to a silo. The complication here is that the grazing area overlaps around the silo (i.e., in general, the tether is longer than one half the circumference of the silo): the goat can only eat the grass once, he can't eat it twice. The answer to the problem as proposed was given in the 1749 issue of the magazine by a Mr. Heath, and stated as 76,257.86 sq.yds. which was arrived at partly by "trial and a table of logarithms". The answer is not so accurate as the number of digits of precision would suggest. No analytical solution was provided.

A useful approximation

Let tether length R = 160 yds. and silo radius r = R/(2π) yds. The involute in the fourth quadrant is a nearly circular arc. One can imagine a circular segment with the same perimeter (arc length) would enclose nearly the same area; the radius and therefore the area of that segment could be readily computed. The arc length of an involute is given by [math]\displaystyle{ \tfrac{1}{2}r\theta^2 }[/math] so the arc length |FG| of the involute in the fourth quadrant is [math]\displaystyle{ \tfrac{1}{2} r\Big[\theta^2\Big]^{2\pi}_{\tfrac{3\pi}{2}-\varphi} }[/math]. Let c be the length of an arc segment of the involute between the y-axis and a vertical line tangent to the silo at θ = 3π/2; it is the arc subtended by Φ. [math]\displaystyle{ c\approx r }[/math] (while the arc is minutely longer than r, the difference is negligible). So [math]\displaystyle{ |FG|=\tfrac{1}{2}\cdot 25.46\cdot(6.28^2-4.71^2)+25.46=245.38 }[/math]. The arc length of a circular arc is [math]\displaystyle{ r\theta }[/math] and θ here is π/2 radians of the fourth quadrant, so [math]\displaystyle{ r\tfrac{\pi}{2}=245.38 }[/math], r the radius of the circular arc is [math]\displaystyle{ 156.21 }[/math] and the area of the circular segment bounded by it is [math]\displaystyle{ \tfrac{\pi r^2}{4}=19165 }[/math]. The area of the involute excludes half the area of the silo (1018.61) in the fourth quadrant, so its approx area is 18146, and the grazable area including the half circle of radius R, ([math]\displaystyle{ \tfrac{1}{2}\pi R^2 }[/math]) totals [math]\displaystyle{ 2\cdot 18146+40212=76505 }[/math]. That is 249 sq.yds. greater than the correct area of 76256, an error of just 0.33%. This method of approximating may not be quite so good for angles < 3π/2 of the involute.

If it matters, there is a constructive way to obtain a quick and very accurate estimate of [math]\displaystyle{ \varphi }[/math]: draw a diagonal from point [math]\displaystyle{ \tfrac{3}{2} \pi }[/math] on the circumference of the pond to its intersection on the y-axis. The length of the diagonal is 120yds. because it is [math]\displaystyle{ \tfrac{3}{4} }[/math] of the tether. So the other leg of the triangle, the hypotenuse as drawn, is [math]\displaystyle{ \sqrt{120^2-r^2}=117.27 }[/math] yds. So [math]\displaystyle{ \varphi\approx \sin^{-1}(\tfrac{r}{117.27})=.219 }[/math] radians, rounded to three places.

Solution by integrating with polar coordinates

Find the area between a circle and its involute over an angle of 2π to −2π excluding any overlap. In Cartesian coordinates, the equation of the involute is transcendental; doing a line integral there is hardly feasible. A more felicitous approach is to use polar coordinates (z,θ). Because the "sweep" of the area under the involute is bounded by a tangent line (see diagram and derivation below) which is not the boundary ([math]\displaystyle{ \overline {qF} }[/math]) between overlapping areas, the decomposition of the problem results in four computable areas: a half circle whose radius is the tether length (A₁); the area "swept" by the tether over an angle of 2π (A₂); the portion of area A₂ from θ = 0 to the tangent line segment [math]\displaystyle{ \overline{tF} }[/math] (A₃); and the wedge area qFtq (A₄). So, the desired area A is A₁ + (A₂ − A₃ + A₄) · 2. The area(s) required to be computed are between two quadratic curves, and will necessarily be an integral or difference of integrals.

The primary parameters of the problem are [math]\displaystyle{ R }[/math], the tether length defined to be 160yds, and [math]\displaystyle{ r }[/math], the radius of the silo. There is no necessary relationship between [math]\displaystyle{ R }[/math] and [math]\displaystyle{ r }[/math], but here [math]\displaystyle{ r=\tfrac{160}{2\pi} }[/math] is the radius of the circle whose circumference is [math]\displaystyle{ R }[/math]. If one defines the point of tethering [math]\displaystyle{ v }[/math] (see diagram, above) as the origin with the circle representing the circumference of the pond below the x-axis, and [math]\displaystyle{ F }[/math] on the y-axis below the circle representing the point of intersection of the tether when wound clockwise and counterclockwise, let [math]\displaystyle{ t }[/math] be a point on the circle such that the tangent at [math]\displaystyle{ t }[/math] intersects [math]\displaystyle{ F }[/math], and [math]\displaystyle{ |\overset{\frown}{vt}| }[/math] + [math]\displaystyle{ |\overline {tF}| }[/math] Is the length of the tether. Let [math]\displaystyle{ q }[/math] be the point of intersection of the circumference of the pond on the y-axis (opposite to [math]\displaystyle{ v }[/math]) below the origin. Then let acute [math]\displaystyle{ \angle tFq }[/math] be [math]\displaystyle{ \varphi }[/math].

The area under the involute is a function of [math]\displaystyle{ R^3 }[/math] because it is an integral over a quadratic curve. The area has a fixed boundary defined by the parameter [math]\displaystyle{ r }[/math] (i.e. the circumference of the silo). In this case the area is inversely proportional to [math]\displaystyle{ r }[/math], i.e. the larger [math]\displaystyle{ r }[/math], the smaller the area of the integral, and the circumference is a linear function of [math]\displaystyle{ r }[/math] ([math]\displaystyle{ 2\pi r }[/math]). So we seek an expression for the area under the involute [math]\displaystyle{ E = f(R^3/r) }[/math].

First, the area A₁ is a half circle of radius [math]\displaystyle{ R }[/math] so [math]\displaystyle{ A_1=\tfrac{1}{2}\pi \cdot R^2. }[/math]

Next, find the angle [math]\displaystyle{ \varphi }[/math] which will be used in the limits of the integrals below. Let [math]\displaystyle{ x = |\overline{tF}| }[/math]. [math]\displaystyle{ \varphi }[/math] is complementary to the opposite angle of the triangle whose right angle is at point t; and also complementary to that angle in the third quadrant of the circle. [math]\displaystyle{ |\overline{tF}| }[/math] is the unrolled arc [math]\displaystyle{ \overset{\frown}{vqt} }[/math], so its arclength is [math]\displaystyle{ r\theta_x }[/math]. So [math]\displaystyle{ \theta_x=\tfrac{x}{r} }[/math] and [math]\displaystyle{ \theta_x=\tfrac{3\pi}{2}-\varphi }[/math], so [math]\displaystyle{ \varphi=\tfrac{3\pi}{2} - \tfrac{x}{r} }[/math]. Finally,[math]\displaystyle{ \tan(\varphi)=\tfrac{r}{x} }[/math] and the following equation is obtained: [math]\displaystyle{ r = x \cdot \tan {(\tfrac{3}{2}\pi - \tfrac{x}{r})} }[/math]. That is a transcendental equation that can only be solved by trial-and-error, polynomial expansion, or an iterative procedure like Newton–Raphson. [math]\displaystyle{ \varphi\approx 0.21900[+0,-0.00003] }[/math].

Next compute the area between the circumference of the pond and involute. Compute the area in the tapering "tail" of the involute, i.e. the overlapped area (note, on account of the tangent tF, that this area includes the wedge section, area A₄, which will have to be added back in during the final summation). Recall that the area of a circular sector is [math]\displaystyle{ \tfrac {1}{2}r^2 \theta }[/math] if the angle is in radians. Imagine an infinitely thin circular sector from [math]\displaystyle{ o }[/math] to [math]\displaystyle{ p }[/math] subtended by an infinitely small angle [math]\displaystyle{ \Delta\theta }[/math]. Tangent to [math]\displaystyle{ o }[/math], there is a corresponding infinitely thin sector of the involute from [math]\displaystyle{ m }[/math] to [math]\displaystyle{ n }[/math] subtending the same infinitely small angle [math]\displaystyle{ \Delta\theta }[/math]. The area of this sector is [math]\displaystyle{ \tfrac{1}{2}z^2\Delta\theta }[/math] where [math]\displaystyle{ z }[/math] is the radius at some angle [math]\displaystyle{ \theta_i }[/math], which is [math]\displaystyle{ r\theta_i }[/math], the arc length of the circle so far "unwrapped" at angle [math]\displaystyle{ \theta_i }[/math]. The area under the involute is the sum of all the infinitely many infinitely thin sectors [math]\displaystyle{ i }[/math] through some angle [math]\displaystyle{ \theta_n }[/math]. This sum is

[math]\displaystyle{ \lim_{n\to\infty} \sum_{i=1}^n \tfrac{1}{2}(r\theta_i)^2 \Delta\theta = \int_{\frac {3\pi}{2}-\varphi}^{2\pi} \frac{1}{2}(r\theta)^2 \ d\theta = \frac{1}{2}r^2 \left[\frac {\theta^3}{3}\right]_{\frac {3\pi}{2}-\varphi}^{2\pi} = \frac {4}{3}r^2\pi^3 - \frac{1}{2}r^2\frac{(\frac {3\pi}{2}-\varphi)^3}{3}. }[/math]

The bounds of the integral represent the area under the involute in the fourth quadrant between [math]\displaystyle{ \overline{tF} }[/math] and [math]\displaystyle{ \overline{vG} }[/math]. The angle is measured on the circle, not on the involute, so it is less than [math]\displaystyle{ \tfrac{3\pi}{2} }[/math] by some angle designated [math]\displaystyle{ \varphi }[/math]. [math]\displaystyle{ \varphi }[/math] Is not given, and must be determined indirectly. Unfortunately, there is no way to simplify the latter term representing the lower bound of the eval expression because [math]\displaystyle{ \varphi }[/math] is not a rational fraction of [math]\displaystyle{ \pi }[/math], so it may as well be substituted and evaluated at once (factoring out [math]\displaystyle{ \pi }[/math] preemptively): [math]\displaystyle{ \tfrac{1}{6}r^2\pi^3\cdot (\tfrac{27}{8}-K(\tfrac{\varphi}{\pi})) }[/math] which for expository reasons can be rewritten [math]\displaystyle{ A_2-A_3=\tfrac{R^3}{6r}-\tfrac{R^3}{48r}\cdot (\tfrac{27}{8}-K(\tfrac{\varphi}{\pi})) }[/math]. It seems apropo to merge a factor of [math]\displaystyle{ \tfrac{1}{8} }[/math] into the constant term to get a common denominartor for the terms, so [math]\displaystyle{ A_2-A_3=\tfrac{R^3}{6r}-\tfrac{R^3}{6r}\cdot (\tfrac{27}{64}-K(\tfrac{\varphi}{8\pi})) }[/math]. [math]\displaystyle{ K }[/math] is dominated by a linear term from the integration, so may be written, [math]\displaystyle{ K=\tfrac{27}{32\pi}\varphi-Z }[/math] where [math]\displaystyle{ Z }[/math] is a non-zero positive but negligible quantity.

A₄ is the area of the peculiar wedge [math]\displaystyle{ \lambda^{tFqt} }[/math]. That area is the area of a right triangle with vertex t, minus the area of a sector bounded by [math]\displaystyle{ \overset{\frown}{tq} }[/math]. [math]\displaystyle{ \lambda=\tfrac{rx}{2}-\tfrac{1}{2}r^2\theta }[/math] where x is |tF| and θ is the angle opposite to Φ in the right angle triangle. So, [math]\displaystyle{ \lambda=\tfrac{1}{2}r(R-r(\tfrac{\pi}{2}+\varphi))-\tfrac{1}{2}r^2(\tfrac{\pi}{2}-\varphi)=\tfrac{1}{2}(rR-r^2\pi) }[/math]. If [math]\displaystyle{ R=2\pi r }[/math], then the area [math]\displaystyle{ \lambda }[/math] of the wedge is [math]\displaystyle{ \tfrac{1}{2}\pi r^2 }[/math] by reduction.

The final summation A₁ + (A₂ − A₃ + A₄) · 2 is [math]\displaystyle{ \tfrac{1}{2}\pi\cdot R^2 + (\tfrac{R^3}{6r} - \tfrac{R^3}{6r}\cdot (\tfrac{27}{64}-(\tfrac{27}{32\pi}\varphi-Z)) + \tfrac{1}{2}(rR-r^2\pi))\cdot 2 }[/math]. All imprecision in the calculation is now uncertainty in [math]\displaystyle{ \varphi }[/math] and the residual [math]\displaystyle{ Z }[/math]. [math]\displaystyle{ \varphi\doteq (\tfrac{R}{r}-\tfrac{\pi}{2}-\varphi)^{-1} }[/math]. That's useful for elucidating the relationships between the parameters. [math]\displaystyle{ \varphi }[/math] is transcendental, so the definition is a recurrence relation. The initial guess [math]\displaystyle{ \varphi }[/math] is a small fraction of [math]\displaystyle{ \pi }[/math]. The numerical answer is [math]\displaystyle{ A=76256 }[/math] rounded up to the nearest square yard. It is worth noting that [math]\displaystyle{ \lim_{\varphi\rightarrow \frac{\pi}{2}} A=\tfrac{1}{2}\pi\cdot R^2 + \tfrac{R^3}{6r} }[/math], which is the answer given for the case where the tether length is half the circumference (or any length such that [math]\displaystyle{ \tfrac{R}{r}\leq\pi }[/math]) of the silo, or no overlap to account for. The goat can eat all but 5% of the area of the great circle defined by its tether length, and half the area it cannot eat is within the perimeter of the pond/silo. The only imprecision in the calculation is that no closed-form representation for [math]\displaystyle{ \varphi }[/math] can be derived from the geometry presented. But small inaccuracies in [math]\displaystyle{ \varphi }[/math] when [math]\displaystyle{ \varphi\ll1 }[/math] don't significantly affect the final result.

Solution by ratio of arc length

Just as the area below a line is proportional to the length of the line between boundaries, and the area of a circular sector is a ratio of the arc length ([math]\displaystyle{ L=r\theta }[/math]) of the sector ([math]\displaystyle{ A=\tfrac{r}{2}\cdot L }[/math]), the area between an involute and its bounding circle is also proportional to the involute's arc length[math]\displaystyle{ L=\tfrac{1}{2}r\theta^2 }[/math]: [math]\displaystyle{ A=\tfrac{r\theta}{3}\cdot L =\tfrac{r^2\theta^3}{6} }[/math] for [math]\displaystyle{ 0\lt \theta\lt \theta_i }[/math]. So the total grazing area is [math]\displaystyle{ A=A_1+(A_2-A_3+A_4)\cdot 2 }[/math]. [math]\displaystyle{ A_1=\tfrac{1}{2}\pi R^2 }[/math]. [math]\displaystyle{ A_2-A_3=\left[\tfrac{r^2\theta^3}{6}\right]_{\tfrac{3\pi}{2}-\varphi}^{2\pi} }[/math].[math]\displaystyle{ A_4=\tfrac{1}{2}\pi r^2 }[/math].

Interior grazing problem

Interior grazing problem with goat tethered at Q

Let [math]\displaystyle{ P }[/math] be the center of a unit circle. A goat/bull/horse is tethered at point [math]\displaystyle{ Q }[/math] on the circumference. How long does the rope [math]\displaystyle{ r }[/math] need to be to allow the animal to graze on exactly one half of the circle's area (white area in diagram, in plane geometry, called a lens)?

A useful approximation

The analytical solutions are arduous and yield a transcendental formula, which must be iterated to produce a numerical approximation. A rough and quick approximation may be suitable for most practical applications and is obtainable by a few simple observations.

First, a chord drawn through the center [math]\displaystyle{ P }[/math] of the enclosing circle perpendicular to the radius at [math]\displaystyle{ Q }[/math] divides the area in half. A radius of the cutting circle equal to the radius of the enclosing circle (i.e. unit length) will sweep much but not all of the area below the chord. A radius of the cutting circle that reaches either end of the chord (i.e. length √2 ≈ 1.414) will sweep out the entire area below the chord plus a thin circular segment above the chord. So a plausible radius length is the mean of the two. But the swept area is expanding quadratically, not linearly, with the radius, so the harmonic mean rather than the arithmetic mean is preferable, i.e. 1.189... This is accurate to within ~2.5% and quite usable in a practical situation, and the procedure is easy to replicate.

A much better approximation may be had by the observation that the relationship between tether length and grazable area is substantially linear in the relevant range [1, √2]:

[math]\displaystyle{ A = 2.2 r - 0.98. }[/math]

At [math]\displaystyle{ A = \pi / 2 }[/math], the approximated tether length [math]\displaystyle{ r }[/math] is 1.15945, an error of just 0.06%.

Solution by calculating the lens area

The area reachable by the animal is in the form of an asymmetric lens, delimited by the two circular arcs.

The area [math]\displaystyle{ A }[/math] of a lens with two circles of radii [math]\displaystyle{ R, r }[/math] and distance between centers [math]\displaystyle{ d }[/math] is

[math]\displaystyle{ \begin{align} A = {} & r^2 \cos^{-1}\left(\frac{d^2 + r^2 - R^2}{2 d r}\right) + R^2 \cos^{-1}\left(\frac{d^2 - r^2 + R^2}{2 d R}\right) \\ & {}- \frac{1}{2} \sqrt{(d + r - R) (d - r + R) (-d + r + R) (d + r + R)}, \end{align} }[/math]

which simplifies in case of [math]\displaystyle{ R = d = 1 }[/math] and one half of the circle area to

[math]\displaystyle{ \frac{1}{2} \pi = r^2 \cos^{-1}\left(\frac{1}{2} r\right) + \cos^{-1}\left(1 - \frac{1}{2} r^2\right) - \frac{1}{2} r \sqrt{4 - r^2}. }[/math]

The equation can only be solved iteratively and results in [math]\displaystyle{ r = 1.1587\ldots }[/math] (sequence A133731 in the OEIS).

Solution using integration

By integrating over the right half of the lens area with

[math]\displaystyle{ \frac{1}{4} \pi = \int_0^{\sqrt{r^2 - \frac{r^4}{4}}} \left(\sqrt{r^2 - t^2} + \sqrt{1 - t^2} - 1\right)\, dt }[/math]

the transcendental equation

[math]\displaystyle{ r = \frac{\pi}{\pi r - \sqrt{4 - r^2} + \left(\frac{4}{r} - 2 r\right) \sin^{-1}\left(\frac{1}{2} r\right)} }[/math]

follows, with the same solution.

Closed-form solution

By using complex analysis methods, Ingo Ullisch obtained a closed-form solution as the cosine of a ratio of two contour integrals:

[math]\displaystyle{ r = 2 \cos\left(\frac{1}{2} \frac{\oint_{|z - 3 \pi / 8| = \pi / 4} z / (\sin z - z \cos z - \pi / 2) \,dz}{\oint_{|z - 3 \pi / 8| = \pi / 4} 1 / (\sin z - z \cos z - \pi / 2)\, dz}\right). }[/math]

3-dimensional extension

3-dimensional case with unit sphere on top and goat sphere below

In the three-dimensional case, point [math]\displaystyle{ Q }[/math] lies on the surface of a unit sphere, and the problem is to find radius [math]\displaystyle{ r }[/math] of the second sphere so that the volume of the intersection body equals exactly half the volume of the unit sphere.

The volume of the unit sphere reachable by the animal has the form of a three-dimensional lens with differently shaped sides and defined by the two spherical caps.

The volume [math]\displaystyle{ V }[/math] of a lens with two spheres of radii [math]\displaystyle{ R, r }[/math] and distance between the centers [math]\displaystyle{ d }[/math] is

[math]\displaystyle{ V = \frac{\pi (R + r - d)^2 \left(d^2 + 2dr - 3r^2 + 2dR + 6rR - 3R^2\right)}{12d}, }[/math]

which simplifies in case of [math]\displaystyle{ R = d = 1 }[/math] and one half of the sphere volume to

[math]\displaystyle{ \frac{1}{2} \cdot \frac{4}{3}\pi = -\frac{1}{4} \pi r^4 + \frac{2}{3} \pi r^3, }[/math]

leading to a solution of [math]\displaystyle{ r = 1.2285\ldots }[/math]

It can be demonstrated that, with increasing dimensionality, the reachable area approaches one half the sphere at the critical length [math]\displaystyle{ r = \sqrt{2} }[/math]. If [math]\displaystyle{ r \lt \sqrt{2} }[/math], the area covered approaches almost none of the sphere; if [math]\displaystyle{ r \gt \sqrt{2} }[/math], the area covered approaches the sphere's entire area.

See also

• Mrs. Miniver's problem, another problem of equalizing areas of circular lunes and lenses

References

• Raymond Clare Archibald: "Involutes of a circle and a pasturage problem." American Mathematical Monthly, 1921 (28), pp. 328–329.
• Marshall Fraser: "A tale of two goats." Mathematics Magazine, 1982 (55), pp. 221–227.
• Jean Jacquelin: "Le problème de l’hyperchèvre." Quadrature, 2003, 49, pp. 6–12. ISSN 1142-2785

External links

• Weisstein, Eric W.. "Goat Problem". http://mathworld.wolfram.com/GoatProblem.html. 
• Mathematician Solves Centuries-Old Grazing Goat Problem Exactly in Quanta Magazine