The full problem

In a previous post I've solved issues of a farmer trying to maximise the area of their field using a given length of fence. Given that was solvable I thought, well what if rather than a wall the farmer was to use a circular lake?

This is also modelable but the answer isn’t as intuitive. I would argue the reason the answer isn’t as intuitive is that the more of the fence circle you use, the more the lake invades the pasture. There must be some benefit in using the edge of the lake though, at least my intuition leads me to believe.

Here I try to go about calculating the systems in the same way as in the previous post.

Let the reflex angle BAD represents the large simple area of the pasture.

The segment ABE represents half the remaining segment, bearing in mind that not all of this is usable.

The area BEFB is subtracted from this and represents the area covered by water that would otherwise make the remainder of the arc.

This can be calculated in a similar way to the circle against a wall example leading us to a formula for the total area of the pasture.

$\\ ABFGA = ABEGA - BEGDB - BFGDB\\\\ ABEGA = A_{acute}\overline{AB}^{2}\\\\ BEGDB = ABEGA - \frac{\overline{AD}\ \overline{BG}}{2}\\\\ BFGDB = CBFGC - \frac{\overline{CD}\ \overline{BG}}{2}\\\\ CBFGC = C_{acute}\overline{BC}^{2}$

We are assuming that the length BC is fixed and given to us. We need to know the radius of the pond/lake to work out the answer. The angle C is, however, dependent on how we arrange the fence.

The angle of C can be calculated using the fact that we can work out from the angle A and the length AB (which can be calculated given we know the length of fence) the length BD.

If we have BD and BC we can use trigonometry to find the angle C.

$\\ \sin = \frac{opposite}{hypotenuse }\\\\ \sin C = \frac{\overline{BD}}{\overline{BC}}$

$ABFGA = \frac{\overline{AD}\ \overline{BG}}{2} - \left ( C_{acute}\overline{BC}^{2} - \frac{\overline{CD}\ \overline{BG}}{2} \right )\\ For\ simplicity\ let:\\ \overline{AB} = R\\ \overline{BC} = r\\ A_{acute} = \theta\\ C_{acute} = \phi\\ \therefore \\\\ \overline{AD} = R\cos\frac{\theta}{2}\\\\ \overline{CD} = r\cos\frac{\phi}{2}\\\\ \overline{BD} = R\sin\frac{\theta}{2} = r\sin\frac{\phi}{2}\\\\ \overline{BG} = 2\overline{BD}\\\\ ABFGA = \frac{R\cos\frac{\theta}{2} \left ( 2 R\sin\frac{\theta}{2} \right ) }{2} - \left ( r^2\phi - \frac{r\cos\frac{\phi}{2} \left ( 2 R\sin\frac{\theta}{2} \right )}{2}\right )\\\\ Expanding\ out\\\\ ABFGA = R^2\sin\frac{\theta}{2}\cos\frac{\theta}{2} - \left ( r^2\phi - r\cos\frac{\phi}{2} R\sin\frac{\theta}{2}\right )\\\\$

One of the problems here is that we have trig functions in 2 angles, however we can convert them to one.

$ABFGA = R^2\sin\frac{\theta}{2}\cos\frac{\theta}{2} - \left ( r^2\phi - r\cos\frac{\phi}{2} R\sin\frac{\theta}{2}\right )\\\\ \overline{BD} = R\sin\frac{\theta}{2} = r\sin\frac{\phi}{2}\\\\ \frac{R}{r}\sin\frac{\theta}{2} = \sin\frac{\phi}{2}\\\\ \sin^2x+\cos^2x = 1\\\\ \cos^2x = 1-sin^2x\\\\ \cos^2\frac{\phi}{2} = 1- sin^2\frac{\phi}{2}\\\\ \cos^2\frac{\phi}{2} = 1- \left ( \frac{R}{r}\sin\frac{\theta}{2} \right )^2\\\\ ABFGA = R^2\sin\frac{\theta}{2}\cos\frac{\theta}{2} - \ r^2\sin^{-1}\left (\sqrt{1- \left ( \frac{R}{r}\sin\frac{\theta}{2} \right )^2 }\right ) \ rR\sin\frac{\theta}{2}\sqrt{1 \left ( \frac{R}{r}\sin\frac{\theta}{2} \right )^2} \\\\$

Along the way in this calculation we have a lot of choice in what we express ourself in terms of. Here we have choses to try and express everything as closely as possible in terms of theta. For the next step we can chose to represent the calculation either in terms of theta, or in terms of R. Both are linked by the equation:

$(2 \pi - \theta)R=S$

Which brings us back to the original question. The length of the fence, S, is fixed. That means that changing the angle and the radius of the circle we make with the fence are interlinked.

Substituting:

We get:

$ABFGA = \frac{S}{2\pi-\theta}^2\sin\frac{\theta}{2}\cos\frac{\theta}{2} - \ r^2\sin^{-1}\left (\sqrt{1- \left ( \frac{S}{r\left (2\pi-\theta \right )}\sin\frac{\theta}{2} \right )^2 }\right ) + r\frac{S}{2\pi-\theta}\sin\frac{\theta}{2}\sqrt{1-\left ( \frac{\frac{S}{2\pi-\theta}}{r}\sin\frac{\theta}{2} \right )^2} \\\\$ $ABFGA = \frac{S}{2\pi-\theta}^2\sin\frac{\theta}{2}\cos\frac{\theta}{2} - \ r^2\sin^{-1}\left (\sqrt{1- \left ( \frac{S}{r\left (2\pi-\theta \right )}\sin\frac{\theta}{2} \right )^2 }\right ) + r\frac{S}{2\pi-\theta}\sin\frac{\theta}{2}\sqrt{1-\left ( \frac{S}{r\left ( 2\pi-\theta \right )}\sin\frac{\theta}{2} \right )^2} \\\\$

Which is starting to get a bit ridiculous. However, every variable other than theta is now constant. Which makes the differentiation process long but procedural.

The following identity makes the task a little simpler by breaking it into 3 separate differentiations.

$\frac{d}{du}(u+v)=\frac{du}{dx}+\frac{dv}{dx}$

$ABFGA = \frac{S}{2\pi-\theta}^2\sin\frac{\theta}{2}\cos\frac{\theta}{2} - \ r^2\sin^{-1}\left (\sqrt{1- \left ( \frac{S}{r\left (2\pi-\theta \right )}\sin\frac{\theta}{2} \right )^2 }\right ) + r\frac{S}{2\pi-\theta}\sin\frac{\theta}{2}\sqrt{1-\left ( \frac{S}{r\left ( 2\pi-\theta \right )}\sin\frac{\theta}{2} \right )^2} \\\\ let:\\ ABFGA = \alpha - \beta + \gamma\\\\ \alpha = \frac{S}{2\pi-\theta}^2\sin\frac{\theta}{2}\cos\frac{\theta}{2}\\\\ \beta = r^2\sin^{-1}\left (\sqrt{1- \left ( \frac{S}{r\left (2\pi-\theta \right )}\sin\frac{\theta}{2} \right )^2 }\right )\\\\ \gamma = r\frac{S}{2\pi-\theta}\sin\frac{\theta}{2}\sqrt{1-\left ( \frac{S}{r\left ( 2\pi-\theta \right )}\sin\frac{\theta}{2} \right )^2}\\\\ \frac{d}{d\theta}\left ( ABFGA\right )= \frac{d\alpha}{d\theta} - \frac{d\beta}{\theta} + \frac{d\gamma}{d\theta}$

Redefining theta

For the sake of my sanity I am redefining theta here to two theta to make life easier:

This need to be back propagated through the equations:

$\alpha = \frac{S}{2\pi-2\theta}^2\sin\theta\cos\theta\\\\ Using\ the\ product\ rule:\\ F(x) = f(x)\cdot g(x)\\ F'(x) = f'(x)g(x) + f(x)g'(x)\\ \alpha = uvw\\\\ u = \frac{S}{2\pi-2\theta}^2\\ v = \sin\theta\\ w =\cos\theta\\\\ \frac{d\alpha}{d\theta} = \frac{du}{d\theta}vw+u\frac{dv}{d\theta}w+uv\frac{dw}{d\theta}\\\\ Using\ the\ reciplical\ rule:\\ g(x)= 1/f(x)\\ g'(x) = - \frac{f'(x)}{f(x)^2}\\ u' = \frac{S^2}{2\left (\pi -\theta \right )^3}\\ v' = \cos\theta\\ w' = -\sin\theta\\$

$\\\frac{d\alpha}{d\theta} = \frac{S^2}{2\left (\pi -\theta \right )^3}\sin\theta\cos\theta+\frac{S}{2\pi-2\theta}^2\cos^2\theta - \frac{S}{2\pi-2\theta}^2sin^2\theta\\ \frac{d\alpha}{d\theta} = S^2\left( \frac{\sin\theta\cos\theta}{2\left (\pi -\theta \right )^3}+ \frac{\cos\theta}{2\pi-2\theta}^2 - \frac{sin\theta}{2\pi-2\theta}^2\right )\\ \frac{d\alpha}{d\theta} = S^2\left( \frac{\sin\theta\cos\theta}{2\left (\pi -\theta \right )^3}+ \frac{\cos^2\theta}{4 \left( \pi-\theta\right )^2} - \frac{sin^2\theta}{4 \left( \pi- \theta \right )^2}\right )\\ \frac{d\alpha}{d\theta} = \frac{S^2}{4\left( \pi-\theta\right )^2}\left( \frac{2\sin\theta\cos\theta}{\pi -\theta }+ \cos^2\theta - sin^2\theta\right )\\$

You can see this is blowing up in it’s complexity.

Let us give beta a go, this one is made more complicated by the inverse trig function.

We have used it before but I think it helps here to make it more explicit. Here is the chain rule:

$\beta = r^2\sin^{-1}\left (\sqrt{1- \left ( \frac{S}{r\left (2\pi-2\theta \right )}\sin\theta \right )^2 }\right )\\\\ \frac{d}{dx}\sin^{-1}x = \frac{1}{\sqrt{1-x^2}}\\ let:\\ \delta = \sqrt{1- \left ( \frac{S}{r\left (2\pi-2\theta \right )}\sin\theta \right )^2 }\\ \therefore \beta = r^2\sin^-1\left(\delta \right )\\\\ \frac{d\beta}{d\theta} = r^2 \frac{1}{\sqrt{1-\delta^2}}\frac{d\delta}{d\theta}\\ \frac{d\delta}{d\theta} = (\frac{1}{2}) \left(1- \left ( \frac{S}{r\left (2\pi-2\theta \right )}\sin\theta \right )^2 \right )^{-\frac{1}{2}} \left (-2 \right ) \left (1- \frac{S\sin\theta}{2r\left (\pi-\theta \right )} \right ) \left (\cos \theta \frac{S}{2r\left (\pi-\theta \right )} + \sin\theta\frac{S}{2r}\frac{-1}{\left (\pi - \theta\right)^2}\right )$

There is a huge amount going on here, but one point of note is that things are a little simplified by the use of delta inside the square root is squared which does make things a little simpler.

First we need to simplify our equation for the differential of delta:

$\\\frac{d\delta}{d\theta} = \frac{1}{\sqrt{1- \left ( \frac{S}{2r\left (\pi-\theta \right )}\sin\theta \right )^2 }} \left (\frac{S\sin\theta}{2r\left (\pi-\theta \right )}-1 \right ) \left(\frac{S}{2r\left (\pi-\theta \right )} \right ) \left (\cos \theta - \sin\theta\frac{1}{\pi - \theta}\right )\\\\ \frac{d\delta}{d\theta} = \left(\sqrt{\frac{S^2}{4r^2\left (\pi-\theta \right )^2}} \right ) \frac{1}{\sqrt{1- \frac{S^2}{4r^2\left (\pi-\theta \right )^2}\sin^2\theta }} \left (\frac{S\sin\theta}{2r\left (\pi-\theta \right )}-1 \right ) \left (\cos \theta - \sin\theta\frac{1}{\pi - \theta}\right )\\\\ \frac{d\delta}{d\theta} = \frac{1}{\sqrt{\frac{S^2}{4r^2\left (\pi-\theta \right )^2}- \sin^2\theta }} \left (\frac{S\sin\theta}{2r\left (\pi-\theta \right )}-1 \right ) \left (\cos \theta - \sin\theta\frac{1}{\pi - \theta}\right )\\\\$

$\\ \frac{d\delta}{d\theta} = \frac{1}{\sqrt{\frac{S^2-2r^2\sin^2\theta \left (\pi-\theta \right )^2}{4r^2\left (\pi-\theta \right )^2} }} \left (\frac{S\sin\theta}{2r\left (\pi-\theta \right )}-1 \right ) \left (\cos \theta - \sin\theta\frac{1}{\pi - \theta}\right )\\\\ \frac{d\delta}{d\theta} = \sqrt{\frac{4r^2\left (\pi-\theta \right )^2}{S^2-2r^2\sin^2\theta \left (\pi-\theta \right )^2} } \left (\frac{S\sin\theta}{2r\left (\pi-\theta \right )}-1 \right ) \left (\cos \theta - \sin\theta\frac{1}{\pi - \theta}\right )\\\\ \frac{d\delta}{d\theta} = \frac{2r\left (\pi-\theta \right )}{\sqrt{{S^2-2r^2\sin^2\theta \left (\pi-\theta \right )^2} }} \left (\frac{S\sin\theta}{2r\left (\pi-\theta \right )}-1 \right ) \left (\cos \theta - \sin\theta\frac{1}{\pi - \theta}\right )\\\\$

$\\ \frac{d\delta}{d\theta} = \frac{2r\left (\pi-\theta \right )}{\sqrt{{S^2-2r^2\sin^2\theta \left (\pi-\theta \right )^2} }} \left (\frac{S\sin\theta}{2r\left (\pi-\theta \right )}-1 \right ) \left (\cos \theta - \frac{\sin\theta}{\pi - \theta}\right )\\\\ \frac{d\delta}{d\theta} = \frac{\left ( S\sin\theta - 2r\left (\pi-\theta \right )\right )\left (2r\left(\pi-\theta \right )\cos \theta 2r\sin\theta \right )}{\sqrt{{S^2-2r^2\sin^2\theta \left (\pi\theta \right )^2} }}$

$\frac{d\beta}{d\theta} = r^2 \frac{1}{\sqrt{1-\delta^2}}\frac{d\delta}{d\theta}\\\\ \frac{d\beta}{d\theta} = r^2 \frac{1}{\frac{S}{2r\left (\pi-\theta \right )}\sin\theta } \frac{2r\left (\pi-\theta \right )}{\sqrt{{S^2-2r^2\sin^2\theta \left (\pi-\theta \right )^2} }} \left (\frac{S\sin\theta}{2r\left (\pi-\theta \right )}-1 \right ) \left (\cos \theta - \frac{\sin\theta}{\pi - \theta}\right )\\\\\\ \frac{d\beta}{d\theta} = r^2 \frac{2r\left (\pi-\theta \right )}{S\sin\theta } \frac{2r\left (\pi-\theta \right )}{\sqrt{{S^2-2r^2\sin^2\theta \left (\pi-\theta \right )^2} }} \left (\frac{S\sin\theta}{2r\left (\pi-\theta \right )}-1 \right ) \left (\cos \theta - \frac{\sin\theta}{\pi - \theta}\right )\\\\$

Onto differentiating the final term of our equation Gamma.

$\gamma = \frac{rS}{2\left ( \pi-\theta\right )}\sin\theta\sqrt{1-\left ( \frac{S}{2r\left ( \pi-\theta\right )}\sin\theta\right )^2}\\\\ \gamma = \frac{S}{4\left ( \pi-\theta\right )^2}\sin\theta \sqrt{4r^2\left ( \pi-\theta\right )^2-\left (S\sin\theta\right )^2}\\\\ a = \frac{S}{4\left ( \pi-\theta\right )^2}\sin\theta\\\\ b = \sqrt{4r^2\left ( \pi-\theta\right )^2-\left (S\sin\theta\right )^2}\\\\ \gamma = a \cdot b\\\\ \frac{d \gamma}{d \theta}=\frac{d(a \cdot b)}{d \theta} = \frac{da}{d \theta} \cdot b + \frac{db}{d \theta} \cdot a \\\\ a_1 = (\pi - \theta)^{-2} \\\\ a_2 = sin \theta \\\\ a = \frac{S}{4} \cdot a_1 \cdot a_2\\\\ \frac{da_1}{d \theta} = 2(\pi - \theta)^{-3}\\\\ \frac{da_2}{d \theta} = cost \theta \\\\ \frac{da}{d \theta} = \frac{S}{4} \cdot \left ( \frac{da_1}{d \theta} \cdot a_2 + a_1 \cdot \frac{da_2}{d \theta}\right )$

$a_1 = (\pi - \theta)^{-2} \\\\ a_2 = sin \theta \\\\ a = \frac{S}{4} \cdot a_1 \cdot a_2\\\\ \frac{da_1}{d \theta} = 2(\pi - \theta)^{-3}\\\\ \frac{da_2}{d \theta} = cos \theta \\\\ \frac{da}{d \theta} = \frac{S}{4} \cdot \left ( \frac{da_1}{d \theta} \cdot a_2 + a_1 \cdot \frac{da_2}{d \theta}\right )\\\\ \frac{da}{d \theta} = \frac{S}{4} \cdot \left( 2(\pi - \theta)^{-3} \cdot sin \theta + cos \theta (\pi - \theta)^{-2} \right ) \\\\ \frac{da}{d \theta} = \frac{S}{4(\pi - \theta)^3} \cdot \left( 2 sin \theta + cos \theta (\pi - \theta) \right )$

$b = \sqrt{4r^2\left ( \pi-\theta\right )^2-\left (S\sin\theta\right )^2}\\\\ b^2 = 4r^2\left ( \pi-\theta\right )^2-\left (S\sin\theta\right )^2 \\ \frac{d b^2}{d \theta} = -(2)4r^2(\pi - \theta) - 2S \cos \theta (S \sin \theta) \\\\ \frac{d b^2}{d \theta} = -8r^2(\pi - \theta) - 2S^2 \sin \theta \cos \theta \\\\ b = f(b^2) \implies \frac{d b}{d \theta} = \frac{d b^2}{d \theta} \cdot f'(b^2) \\\\ \frac{db}{d \theta} = \left ( 8r^2(\pi - \theta) - 2S^2 \sin \theta \cos \theta \right ) \cdot \left ( 4r^2\left ( \pi\theta\right )^2-\left (S\sin\theta\right )^2 \right )^{-\frac{1}{2}} \\\\ \frac{db}{d \theta} = \frac{-8r^2(\pi - \theta) - 2S^2 \sin \theta \cos \theta }{\sqrt{4r^2\left ( \pi-\theta\right )^2-\left (S\sin\theta\right )^2} } \\\\$

$\frac{d \gamma}{d \theta}=\frac{d(a \cdot b)}{d \theta} = \frac{da}{d \theta} \cdot b + \frac{db}{d \theta} \cdot a \\\\ a = \frac{S}{4\left ( \pi-\theta\right )^2}\sin\theta\\\\ b = \sqrt{4r^2\left ( \pi-\theta\right )^2-\left (S\sin\theta\right )^2}\\\\ \frac{da}{d \theta} = \frac{S}{4(\pi - \theta)^3} \cdot \left( 2 sin \theta + cos \theta (\pi - \theta) \right )\\\\ \frac{db}{d \theta} = \frac{-8r^2(\pi - \theta) - 2S^2 \sin \theta \cos \theta }{\sqrt{4r^2\left ( \pi-\theta\right )^2-\left (S\sin\theta\right )^2} } \\\\ \frac{d \gamma}{d \theta}= \frac{S}{4(\pi - \theta)^3} \cdot \left( 2 sin \theta + cos \theta (\pi - \theta) \right ) \cdot \sqrt{4r^2\left ( \pi-\theta\right )^2-\left (S\sin\theta\right )^2} + \frac{-8r^2(\pi - \theta) - 2S^2 \sin \theta \cos \theta }{\sqrt{4r^2\left ( \pi-\theta\right )^2-\left (S\sin\theta\right )^2} } \cdot \frac{S}{4\left ( \pi-\theta\right )^2}\sin\theta \\\\$

$\frac{d \gamma}{d \theta}= \frac{S}{4(\pi - \theta)^3} \cdot \left( 2 sin \theta + cos \theta (\pi - \theta) \right ) \cdot \sqrt{4r^2\left ( \pi-\theta\right )^2-\left (S\sin\theta\right )^2} + \frac{-8r^2(\pi - \theta) - 2S^2 \sin \theta \cos \theta }{\sqrt{4r^2\left ( \pi-\theta\right )^2-\left (S\sin\theta\right )^2} } \cdot \frac{S}{4\left ( \pi-\theta\right )^2}\sin\theta \\\\ \frac{d \gamma}{d \theta}= \frac{S\left( 2 sin \theta + cos \theta (\pi - \theta) \right ) \cdot \sqrt{4r^2\left ( \pi-\theta\right )^2-\left (S\sin\theta\right )^2}}{4(\pi - \theta)^3} \cdot + \frac{-2Sr^2(\pi - \theta) - 2S^3 \sin^2 \theta \cos \theta }{\left ( \pi-\theta\right )^2\sqrt{4r^2\left ( \pi-\theta\right )^2-\left (S\sin\theta\right )^2} } \\\\$

That is where the story ends for now. The sum of the three calculated differentials should be set to zero and the equation solved for theta. Given the nature of the calculation it seems that the solution is transcendental and so there is little hope of getting an analytical solution.

Perhaps there are some tools that are not yet in my mathematical toolbox which will yield a more fruitful search for answers.

Search This Blog

The full problem

Comments

Post a Comment

Popular posts from this blog

Structural engineering with cardboard

The twelve fold way

Removing links from blogger comment automatically