Solving Problems Submitted to MAA Journals (Part 6b)

The following problem appeared in Volume 97, Issue 3 (2024) of Mathematics Magazine.

Two points $P$ and $Q$ are chosen at random (uniformly) from the interior of a unit circle. What is the probability that the circle whose diameter is segment $\overline{PQ}$ lies entirely in the interior of the unit circle?

As discussed in the previous post, I guessed from simulation that the answer is $2/3$ . Naturally, simulation is not a proof, and so I started thinking about how to prove this.

My first thought was to make the problem simpler by letting only one point be chosen at random instead of two. Suppose that the point $P$ is fixed at a distance $t$ from the origin. What is the probability that the point $Q$ , chosen at random, uniformly, from the interior of the unit circle, has the desired property?

My second thought is that, by radial symmetry, I could rotate the figure so that the point $P$ is located at $(t,0)$ . In this way, the probability in question is ultimately going to be a function of $t$ .

There is a very nice way to compute such probabilities since $Q$ is chosen at uniformly from the unit circle. Let $A_t$ be the probability that the point $Q$ has the desired property. Since the area of the unit circle is $\pi(1)^2 = \pi$ , the probability of desired property happening is

$\displaystyle \frac{\hbox{area}(A_t)}{\pi}$ .

So, if I could figure out the shape of $A_t$ , I could compute this conditional probability given the location of the point $P$ .

But, once again, I initially had no idea of what this shape would look like. So, once again, I turned to simulation with Mathematica.

First, a technical detail that I ignored in the previous post. To generate points $(x,y)$ at random inside the unit circle, one might think to let $x = r \cos \theta$ and $y = r \sin \theta$ , where the distance from the origin $r$ is chosen at random between 0 and 1 and the angle $\theta$ is chosen at random from $0 to 2\pi$ . Unfortunately, this simple simulation generates too many points that are close to the origin and not enough that are close to the circle:

To see why this happened, let $R$ denote the distance of a randomly chosen point from the origin. Then the event $R < r$ is the same as saying that the point lies inside the circle centered at the origin with radius $r$ , so that the probability of this event should be

$F(r) = P(R < r) = \displaystyle \frac{\pi r^2}{\pi (1)^2} = r^2$ .

However, in the above simulation, $R$ was chosen uniformly from $[0,1]$ , so that $P(R < r) = r$ . All this to say, the above simulation did not produce points uniformly chosen from the unit circle.

To remedy this, we employ the standard technique of using the inverse of the above function $F(r)$ , which is clearly $F^{-1}(r) = \sqrt{r}$ . In other words, we will chose randomly chosen radius to have the form $R= \sqrt{U}$ , where $U$ is chosen uniformly on $[0,1]$ . In this way,

$P(R < r) = P( \sqrt{U} < r) = P(U < r^2) = r^2$ ,

as required. Making this modification (highlighted in yellow) produces points that are more evenly distributed in the unit circle; any bunching of points or empty spaces are simply due to the luck of the draw.

In the next post, I’ll turn to the simulation of $A_t$ .

Solving Problems Submitted to MAA Journals (Part 6a)

The following problem appeared in Volume 97, Issue 3 (2024) of Mathematics Magazine.

Two points $P$ and $Q$ are chosen at random (uniformly) from the interior of a unit circle. What is the probability that the circle whose diameter is segment $\overline{PQ}$ lies entirely in the interior of the unit circle?

It took me a while to wrap my head around the statement of the problem. In the figure, the points $P$ and $Q$ are chosen from inside the unit circle (blue). Then the circle (pink) with diameter $\overline{PQ}$ has center $M$ , the midpoint of $\overline{PQ}$ . Also, the radius of the pink circle is $MP=MQ$ .

The pink circle will lie entirely the blue circle exactly when the green line containing the origin $O$ , the point $M$ , and a radius of the pink circle lies within the blue circle. Said another way, the condition is that the distance $MO$ plus the radius of the pink circle is less than 1, or

$MO + MP < 1$ .

As a first step toward wrapping my head around this problem, I programmed a simple simulation in Mathematica to count the number of times that $MO + MP < 1$ when points $P$ and $Q$ were chosen at random from the unit circle.

In the above simulation, out of about 61,000,000 attempts, 66.6644% of the attempts were successful. This leads to the natural guess that the true probability is $2/3$ . Indeed, the 95% confidence confidence interval $(0.666524, 0.666764)$ contains $2/3$ , so that the difference of $0.666644$ from $2/3$ can be plausibly attributed to chance.

I end with a quick programming note. This certainly isn’t the ideal way to perform the simulation. First, for a fast simulation, I should have programmed in C++ or Python instead of Mathematica. Second, the coordinates of $P$ and $Q$ are chosen from the unit square, so it’s quite possible for $P$ or $Q$ or both to lie outside the unit circle. Indeed, the chance that both $P$ and $Q$ lie in the unit disk in this simulation is $(\pi/4)^2 \approx 0.617$ , meaning that about $38.3\%$ of the simulations were simply wasted. So the only sense that this was a quick simulation was that I could type it quickly in Mathematica and then let the computer churn out a result. (I’ll talk about a better way to perform the simulation in the next post.)

Solving Problems Submitted to MAA Journals (Part 5e)

The following problem appeared in Volume 96, Issue 3 (2023) of Mathematics Magazine.

Evaluate the following sums in closed form:

$f(x) = \displaystyle \sum_{n=0}^\infty \left( \cos x - 1 + \frac{x^2}{2!} - \frac{x^4}{4!} \dots + (-1)^{n-1} \frac{x^{2n}}{(2n)!} \right)$

and

$g(x) = \displaystyle \sum_{n=0}^\infty \left( \sin x - x + \frac{x^3}{3!} - \frac{x^5}{5!} \dots + (-1)^{n-1} \frac{x^{2n+1}}{(2n+1)!} \right)$ .

By using the Taylor series expansions of $\sin x$ and $\cos x$ and flipping the order of a double sum, I was able to show that

$f(x) = -\displaystyle \frac{x \sin x}{2} \qquad \hbox{and} \qquad g(x) = \frac{x\cos x - \sin x}{2}$ .

I immediately got to thinking: there’s nothing particularly special about $\sin x$ and $\cos x$ for this analysis. Is there a way of generalizing this result to all functions with a Taylor series expansion?

Suppose

$h(x) = \displaystyle \sum_{k=0}^\infty a_k x^k$ ,

and let’s use the same technique to evaluate

$\displaystyle \sum_{n=0}^\infty \left( h(x) - \sum_{k=0}^n a_k x^k \right) = \sum_{n=0}^\infty \sum_{k=n+1}^\infty a_k x^k$

$= \displaystyle \sum_{k=1}^\infty \sum_{n=0}^{k-1} a_k x^k$

$= \displaystyle \sum_{k=1}^\infty k a_k x^k$

$= x \displaystyle \sum_{k=1}^\infty k a_k x^{k-1}$

$= x \displaystyle \sum_{k=1}^\infty \left(a_k x^k \right)'$

$= x \displaystyle \left[ (a_0)' + \sum_{k=1}^\infty \left(a_k x^k \right)' \right]$

$= x \displaystyle \sum_{k=0}^\infty \left(a_k x^k \right)'$

$= x \displaystyle \left( \sum_{k=0}^\infty a_k x^k \right)'$

$= x h'(x)$ .

To see why this matches our above results, let’s start with $h(x) = \cos x$ and write out the full Taylor series expansion, including zero coefficients:

$\cos x = 1 + 0x - \displaystyle \frac{x^2}{2!} + 0x^3 + \frac{x^4}{4!} + 0x^5 - \frac{x^6}{6!} \dots$ ,

so that

$x (\cos x)' = \displaystyle \sum_{n=0}^\infty \left( \cos x - \sum_{k=0}^n a_k x^k \right)$

$-x \sin x= \displaystyle \left(\cos x - 1 \right) + \left(\cos x - 1 + 0x \right) + \left( \cos x -1 + 0x + \frac{x^2}{2!} \right) + \left( \cos x -1 + 0x + \frac{x^2}{2!} + 0x^3 \right)$

$\displaystyle + \left( \cos x -1 + 0x + \frac{x^2}{2!} + 0x^3 - \frac{x^4}{4!} \right) + \left( \cos x -1 + 0x + \frac{x^2}{2!} + 0x^3 - \frac{x^4}{4!} + 0x^5 \right) \dots$

After dropping the zero terms and collecting, we obtain

$-x \sin x= \displaystyle 2 \left(\cos x - 1 \right) + 2 \left( \cos x -1 + \frac{x^2}{2!} \right) + 2 \left( \cos x -1 + \frac{x^2}{2!} - \frac{x^4}{4!} \right) \dots$

$-x \sin x = 2 f(x)$

$\displaystyle -\frac{x \sin x}{2} = f(x)$ .

A similar calculation would apply to any even function $h(x)$ .

We repeat for

$h(x) = \sin x = 0 + x + 0x^2 - \displaystyle \frac{x^3}{3!} + 0x^4 + \frac{x^5}{5!} + 0x^6 - \frac{x^7}{7!} \dots$ ,

so that

$x (\sin x)' = (\sin x - 0) + (\sin x - 0 - x) + (\sin x - 0 - x + 0x^2)$

$+ \displaystyle \left( \sin x - 0 - x + 0x^2 + \frac{x^3}{3!} \right) + \left( \sin x - 0 - x + 0x^2 + \frac{x^3}{3!} + 0x^4 \right)$

$+ \displaystyle \left( \sin x - 0 - x + 0x^2 + \frac{x^3}{3!} + 0x^4 - \frac{x^5}{5!} \right) + \left( \sin x - 0 - x + 0x^2 + \frac{x^3}{3!} + 0x^4 - \frac{x^5}{5!} + 0 x^6 \right) \dots$ ,

$x\cos x - \sin x = 2(\sin x - x) + \displaystyle 2\left(\sin x - x + \frac{x^3}{3!} \right) + 2 \left( \sin x - x + \frac{x^3}{3!} - \frac{x^5}{5!} \right) \dots$

$x \cos x - \sin x = 2 g(x)$

$\displaystyle \frac{x \cos x - \sin x}{2} = g(x)$ .

A similar argument applies for any odd function $h(x)$ .

Solving Problems Submitted to MAA Journals (Part 5d)

The following problem appeared in Volume 96, Issue 3 (2023) of Mathematics Magazine.

Evaluate the following sums in closed form:

$f(x) = \displaystyle \sum_{n=0}^\infty \left( \cos x - 1 + \frac{x^2}{2!} - \frac{x^4}{4!} \dots + (-1)^{n-1} \frac{x^{2n}}{(2n)!} \right)$

and

$g(x) = \displaystyle \sum_{n=0}^\infty \left( \sin x - x + \frac{x^3}{3!} - \frac{x^5}{5!} \dots + (-1)^{n-1} \frac{x^{2n+1}}{(2n+1)!} \right)$ .

In the previous two posts, I showed that

$f(x) = - \displaystyle \frac{x \sin x}{2} \qquad \hbox{and} \qquad g(x) = \displaystyle \frac{x \cos x - \sin x}{2}$ ;

the technique that I used was using the Taylor series expansions of $\sin x$ and $\cos x$ to write $f(x)$ and $g(x)$ as double sums and then interchanging the order of summation.

In the post, I share an alternate way of solving for $f(x)$ and $g(x)$ . I wish I could take credit for this, but I first learned the idea from my daughter. If we differentiate $g(x)$ , we obtain

$g'(x) = \displaystyle \sum_{n=0}^\infty \left( [\sin x]' - [x]' + \left[\frac{x^3}{3!}\right]' - \left[\frac{x^5}{5!}\right]' \dots + \left[(-1)^{n-1} \frac{x^{2n+1}}{(2n+1)!}\right]' \right)$

$= \displaystyle \sum_{n=0}^\infty \left( \cos x - 1 + \frac{3x^2}{3!} - \frac{5x^4}{5!} \dots + (-1)^{n-1} \frac{(2n+1)x^{2n}}{(2n+1)!} \right)$

$= \displaystyle \sum_{n=0}^\infty \left( \cos x - 1 + \frac{3x^2}{3 \cdot 2!} - \frac{5x^4}{5 \cdot 4!} \dots + (-1)^{n-1} \frac{(2n+1)x^{2n}}{(2n+1)(2n)!} \right)$

$= \displaystyle \sum_{n=0}^\infty \left( \cos x - 1 + \frac{x^2}{2!} - \frac{x^4}{4!} \dots + (-1)^{n-1} \frac{x^{2n}}{(2n)!} \right)$

$= f(x)$ .

Something similar happens when differentiating the series for $f(x)$ ; however, it’s not quite so simple because of the $-1$ term. I begin by separating the $n=0$ term from the sum, so that a sum from $n =1$ to $\infty$ remains:

$f(x) = \displaystyle \sum_{n=0}^\infty \left( \cos x - 1 + \frac{x^2}{2!} - \frac{x^4}{4!} \dots + (-1)^{n-1} \frac{x^{2n}}{(2n)!} \right)$

$= (\cos x - 1) + \displaystyle \sum_{n=1}^\infty \left( \cos x - 1 + \frac{x^2}{2!} - \frac{x^4}{4!} \dots + (-1)^{n-1} \frac{x^{2n}}{(2n)!} \right)$ .

I then differentiate as before:

$f'(x) = (\cos x - 1)' + \displaystyle \sum_{n=1}^\infty \left( [\cos x - 1]' + \left[ \frac{x^2}{2!} \right]' - \left[ \frac{x^4}{4!} \right]' \dots + \left[ (-1)^{n-1} \frac{x^{2n}}{(2n)!} \right]' \right)$

$= -\sin x + \displaystyle \sum_{n=1}^\infty \left( -\sin x + \frac{2x}{2!} - \frac{4x^3}{4!} \dots + (-1)^{n-1} \frac{(2n) x^{2n-1}}{(2n)!} \right)$

$= -\sin x + \displaystyle \sum_{n=1}^\infty \left( -\sin x + \frac{2x}{2 \cdot 1!} - \frac{4x^3}{4 \cdot 3!} \dots + (-1)^{n-1} \frac{(2n) x^{2n-1}}{(2n)(2n-1)!} \right)$

$= -\sin x + \displaystyle \sum_{n=1}^\infty \left( -\sin x + x - \frac{x^3}{3!} + \dots + (-1)^{n-1} \frac{x^{2n-1}}{(2n-1)!} \right)$

$= -\sin x - \displaystyle \sum_{n=1}^\infty \left( \sin x - x + \frac{x^3}{3!} + \dots - (-1)^{n-1} \frac{x^{2n-1}}{(2n-1)!} \right)$ .

At this point, we reindex the sum. We make the replacement $k = n - 1$ , so that $n = k+1$ and $k$ varies from $k=0$ to $\infty$ . After the replacement, we then change the dummy index from $k$ back to $n$ .

$f'(x) = -\sin x - \displaystyle \sum_{k=0}^\infty \left( \sin x - x + \frac{x^3}{3!} + \dots - (-1)^{(k+1)-1} \frac{x^{2(k+1)-1}}{(2(k+1)-1)!} \right)$

$= -\sin x - \displaystyle \sum_{k=0}^\infty \left( \sin x - x + \frac{x^3}{3!} + \dots - (-1)^{k} \frac{x^{2k+1}}{(2k+1)!} \right)$

$= -\sin x - \displaystyle \sum_{n=0}^\infty \left( \sin x - x + \frac{x^3}{3!} + \dots - (-1)^{n} \frac{x^{2n+1}}{(2n+1)!} \right)$

With a slight alteration to the $(-1)^n$ term, this sum is exactly the definition of $g(x)$ :

$f'(x)= -\sin x - \displaystyle \sum_{n=0}^\infty \left( \sin x - x + \frac{x^3}{3!} + \dots - (-1)^1 (-1)^{n-1} \frac{x^{2n+1}}{(2n+1)!} \right)$

$= -\sin x - \displaystyle \sum_{n=0}^\infty \left( \sin x - x + \frac{x^3}{3!} + \dots + (-1)^{n-1} \frac{x^{2n+1}}{(2n+1)!} \right)$

$= -\sin x - g(x)$ .

Summarizing, we have shown that $g'(x) = f(x)$ and $f'(x) = -\sin x - g(x)$ . Differentiating $f'(x)$ a second time, we obtain

$f''(x) = -\cos x - g'(x) = -\cos x - f(x)$

$f''(x) + f(x) = -\cos x$ .

This last equation is a second-order nonhomogeneous linear differential equation with constant coefficients. A particular solution, using the method of undetermined coefficients, must have the form $F(x) = Ax\cos x + Bx \sin x$ . Substituting, we see that

$[Ax \cos x + B x \sin x]'' + A x \cos x + Bx \sin x = -\cos x$

$-2A \sin x - Ax \cos x + 2B \cos x - B x \sin x + Ax \cos x + B x \sin x = -\cos x$

$-2A \sin x + 2B \cos x = -\cos x$

We see that $A = 0$ and $B = -1/2,$ which then lead to the particular solution

$F(x) = -\displaystyle \frac{1}{2} x \sin x$

Since $\cos x$ and $\sin x$ are solutions of the associated homogeneous equation $f''(x) + f(x) = 0$ , we conclude that

$f(x) = c_1 \cos x + c_2 \sin x - \displaystyle \frac{1}{2} x \sin x$ ,

where the values of $c_1$ and $c_2$ depend on the initial conditions on $f$ . As it turns out, it is straightforward to compute $f(0)$ and $f'(0)$ , so we will choose $x=0$ for the initial conditions. We observe that $f(0)$ and $g(0)$ are both clearly equal to 0, so that $f'(0) = -\sin 0 - g(0) = 0$ as well.