Sphere on a freely spinning turntable (again)..

After my initial post on this topic, I still felt that there is some more to the problem than what we've explored. I kept returning to this problem hoping to find something more and am glad I did.

We start with a following 'force' equation that we saw earlier.

$\mathbf{p}''=\alpha p_0 \mathbf{k}\times\widehat{\mathbf{p}}'$

Taking a dot product on both sides with $\mathbf{p}'$, we have

$\mathbf{p}'\cdot\mathbf{p}''=\alpha p_0 \mathbf{p}'\cdot(\mathbf{k}\times\widehat{\mathbf{p}}')$

Permuting the RHS and integrating the above, we see that

$\displaystyle \frac{1}{2}\left\vert{\mathbf{p}'}\right\vert^2+\alpha p_0 \int\mathbf{k}\cdot(\mathbf{p}'\times\widehat{\mathbf{p}}')\,dt=\text{const.}$

For any vector $\mathbf{q}$, with straightforward differentiation, we can show that

$\displaystyle \mathbf{q}'\times\widehat{\mathbf{q}}'=\frac{\mathbf{q}\cdot\mathbf{q}'}{\left\vert\mathbf{q}\right\vert^3}(\mathbf{q}\times\mathbf{q}')$

Using this result, the constant vector $\mathbf{m}$ from our earlier post and the fact that $\mathbf{p}\cdot\mathbf{k}=\sqrt{\delta/\alpha}$, we see that

$\displaystyle \mathbf{k}\cdot(\mathbf{p}'\times\widehat{\mathbf{p}}')=\frac{\mathbf{p}\cdot\mathbf{p}'}{p^3}\left(\mathbf{k}\cdot\mathbf{m}-\frac{\delta p_0}{p}\right)$

where $p=\left\vert\mathbf{p}\right\vert$.

Let $u=1/p$. It is easy to see that $u'=-(\mathbf{p}\cdot\mathbf{p}')/p^3$. This makes it almost trivial to integrate the above expression.

With the above, our 'energy' equation becomes

$\displaystyle E=\frac{1}{2}\left\vert{\mathbf{p}'}\right\vert^2-\alpha (\mathbf{k}\cdot\mathbf{m}) \frac{p_0}{p}+\frac{\delta\alpha}{2}\frac{p_0^2}{p^2}$

Now comes the interesting part. If we interpret the first term of the RHS as a measure of kinetic energy, the remaining terms gives us the 'effective potential' energy of the system. That is,

$\displaystyle V_{\text{eff}}(p)=-\alpha (\mathbf{k}\cdot\mathbf{m}) \frac{p_0}{p}+\frac{\delta\alpha}{2}\frac{p_0^2}{p^2}$

But that is exactly in the same form as that of the effective potential of the radial equation of an inverse square central force (Central forces). Therefore, the problem of a sphere on a freely spinning turntable is a Kepler problem in disguise which helps us in utilizing many known results.

For example, it known that the orbital period of a kepler problem can be written in terms of the 'energy' and the coefficient of the '$1/r$' term in the effective potential. Using the same, we can see that the orbital period $T$ for our problem is given by

$$\displaystyle T=\frac{2\pi\text{ }\alpha\text{ }p_0\text{ }(\mathbf{k}\cdot\mathbf{m})}{\left\vert2E\right\vert^{3/2}}$$

UPDATE (4-Aug-2026):

From (13) of Wessecker's paper, we can easily infer that

$L=\alpha p_0 p + \mathbf{k}\cdot(\mathbf{p}\times\mathbf{c})$

where $L$ is the reduced angular momentum about the $z$-axis which is shown to be a constant and $\mathbf{c}$ is the constant of integration in (10).

Now we define a vectorized form of $L$ such that,

$\mathbf{L}=\alpha p_0 p \mathbf{k}+\mathbf{p}\times\mathbf{c}$

We now show that $\mathbf{L}$ is a constant of motion by showing its components perpendicular and parallel to $\mathbf{k}$ are constants.

It's easy to see that $\mathbf{L}\cdot\mathbf{k}=L$ is a constant. Now,

$\mathbf{L}_{\perp}=\mathbf{k}\times(\mathbf{L}\times\mathbf{k})$

Using the fact that $\mathbf{c}\cdot\mathbf{k}=0$, we have

$\mathbf{k}\times(\mathbf{p}\times\mathbf{c})=(\mathbf{p}\cdot\mathbf{k})\mathbf{c}$

which is clearly a constant. Thus we see that $\mathbf{L}_{\perp}$, and hence, $\mathbf{L}$ is a constant.

We now start with the velocity equation which is (10) of Wessecker's paper expressed in terms of $\mathbf{p}$.

$\displaystyle\mathbf{p}'=\alpha p_0 \text{ }\mathbf{k}\times\hat{\mathbf{p}}+\mathbf{c}=\frac{\alpha p_0}{p} \text{ }\mathbf{k}\times\mathbf{p}+\mathbf{c}$.

Crossing the velocity equation with $\mathbf{p}$, with simple manipulations its easy to see in fact that $\mathbf{L}=\mathbf{m}$.

Dotting the velocity equation with itself and simplifying, we have,

$\displaystyle\vert \mathbf{p}'\vert^2=\vert\mathbf{c}\vert^2-\alpha^2 p_0^2+\frac{2\alpha p_0 L}{p}-\alpha\delta\frac{p_0^2}{p^2}$

Comparing with our energy equation, we see that

$\displaystyle E=\frac{1}{2}(\vert\mathbf{c}\vert^2-\alpha^2 p_0^2)=\frac{1}{2}\vert \mathbf{p}'\vert^2-\frac{\alpha p_0 L}{p}+\frac{\alpha\delta p_0^2}{2p^2}$

Hope you enjoyed the discussion. See ya later.

Until then

Yours Aye

Me

A Trigonometric Interpretation of an Infinite Product for the Lemniscate Constant

It is well known that Viete's infinite product formula for $\pi$ has a lemniscatic analog as showed by Levin in A Geometric Interpretation of an Infinite Product for the Lemniscate Constant.

$\displaystyle \frac{2}{\pi} = \sqrt{\frac{1}{2}} \cdot \sqrt{\frac{1}{2} + \frac{1}{2}\sqrt{\frac{1}{2}}} \cdot \sqrt{\frac{1}{2} + \frac{1}{2}\sqrt{\frac{1}{2} + \frac{1}{2}\sqrt{\frac{1}{2}}}} \cdots$

$\displaystyle \frac{2}{\varpi} = \sqrt{\frac{1}{2}} \cdot \sqrt{\frac{1}{2} + \frac{\frac{1}{2}}{\sqrt{\frac{1}{2}}}} \cdot \sqrt{\frac{1}{2} + \frac{\frac{1}{2}}{\sqrt{\frac{1}{2} + \frac{\frac{1}{2}}{\sqrt{\frac{1}{2}}}}}} \cdots ,$

Unfortunately, Levin's proof of this formula gets into Weierstrass with complex arguments to arrive at the result. This post is attempt to give a simpler proof with just trigonometric identities. We'll use the lemniscate sine and cosine functions but will reparametrize them in such a way that the proof remains simple.

The lemniscate elliptic functions, $\text{sl}(s)$ and $\text{cl}(s)$ are close analogues of the classical sine and cosine functions as can be seen in the figures in Relation to geometric shapes. The lemniscate curve was made immortal when Fagnano first found that they obey a doubling identity. That is,

$\displaystyle \text{sl}(2s)^2=\frac{4\text{sl}(s)^2(1-\text{sl}(s)^4)}{(1+\text{sl}(s)^4)^2}$

It can also easily seen that $\text{sl}'(s)^2=1-\text{sl}(s)^4$ where prime denotes differential w.r.t. $s$.

Surprisingly, this is all we need for our proof from the lemniscate. To begin, let's parametrise the lemniscate sine function with a parameter $t$ such that

$\text{sl}(s/2)=\sqrt{\tan(t)}$

Then using Fagnano's formula and double angle formulas of sine and cosine using tangents, we can quickly see that

$\text{sl}(s)=\sqrt{\sin(4t)}$

From here, it is straightforward that $\text{sl}'(s)=\cos(4t)$ and $\text{sl}'(s/2)=\sqrt{1-\tan^2(t)}$

Now with simple trigonometric identities, we have

$\displaystyle \frac{1}{\text{sl}'(s/2)}=\frac{\cos(t)}{\sqrt{\cos(2t)}}=\sqrt{\frac{1+\cos(2t)}{2\cos(2t)}}=\sqrt{\frac{1+\sec{2t}}{2}}$

But we know that $\text{sl}'(s)=\cos(4t)$ and therefore,

$\displaystyle \cos(2t)=\sqrt{\frac{1+\cos(4t)}{2}}=\sqrt{\frac{1+\text{sl}'(s)}{2}}$

Using this expression, we have

$\displaystyle\frac{1}{\text{sl}'(s/2)}=\sqrt{\frac{1+\left(\sqrt{\frac{1+\text{sl}'(s)}{2}}\right)^{-1}}{2}}$

To avoid this nested square roots, we introduce a new function such that $\hat{\text{sl}}(s)=\cos(2t)$. We then have,

$\displaystyle\text{sl}'(s/2)^{-1}=\sqrt{\frac{1+\hat{\text{sl}}(s)^{-1}}{2}}$ and $\displaystyle\hat{\text{sl}}(s)=\sqrt{\frac{1+\text{sl}'(s)}{2}}$

Now comes the interesting part. From our definitions of the lemniscate functions in terms of $t$,

$\displaystyle \frac{\text{sl}(s)}{2\text{sl}(s/2)}=\frac{\sqrt{\sin(4t)}}{2\sqrt{\tan(t)}}=\cos(t)\sqrt{\cos(2t)}$

Also,

$\displaystyle\frac{\hat{\text{sl}}(s)}{\text{sl}'(s/2)}=\frac{\cos(2t)}{\sqrt{1-\tan^2(t)}}=\cos(t)\sqrt{\cos(2t)}$

Therefore, we see that

$\displaystyle \frac{\text{sl}(s)}{2\text{sl}(s/2)}=\frac{\hat{\text{sl}}(s)}{\text{sl}'(s/2)}$

Repeated application of this relation $n$ times,

$\displaystyle \frac{\text{sl}(s)}{2^n\text{sl}(s/2^n)}=\frac{\hat{\text{sl}}(s)}{\text{sl}'(s/2)}\frac{\hat{\text{sl}}(s/2)}{\text{sl}'(s/4)}\cdots\frac{\hat{\text{sl}}(s/2^{n-1})}{\text{sl}'(s/2^n)}$

As $n \to \infty$ and noting that $\lim_{n \to \infty}2^n\text{sl}(s/2^n)=s$, we have

$\displaystyle \frac{\text{sl}(s)}{s}=\hat{\text{sl}}(s)\cdot\text{sl}'(s/2)^{-1}\cdot\hat{\text{sl}}(s/2)\cdot\text{sl}'(s/4)^{-1}\cdot\hat{\text{sl}}(s/4)\cdot\text{sl}'(s/8)^{-1}\cdots$

Based on the relations we saw before, if $a$ and $b$ are two consecutive terms on the RHS above, we see that

$\displaystyle b=\sqrt{\frac{1+a^{-1}}{2}}$

We have $\text{sl}(\varpi/2)=1$, $\text{sl}'(\varpi/2)=0$ and $\hat{\text{sl}}(\varpi/2)=1/\sqrt{2}$. Using these values in the infinite product above, we get the lemniscatic analogue of Viete's formula.

Hope you enjoyed this. See ya later.

Until then

Yours Aye

Me

Blending the Birthday Problem and the Coupon Collector Problem

The Birthday Problem and Coupon Collector problem are well known problems in probability that are like duals of each other. While the former deals with 'first collision', the latter deals with the 'complete collection'. I believe they are too famous even warrant further introduction here.

In this post, we attempt to solve a blend of these two problems. We start with the following question: If we are collecting coupons from among $n$ different types of coupons, what is the expected numbers of coupons that are needed to complete the collection given that the first collision occurred at the $c$'th coupon.

Let $K$ be the random variable indicating the number of coupons required to complete the collection and $C$ be the random variable indicated the number of coupons for the first collision. Then, we are looking for $\mathbb{E}(K|C=c)$.

This is not very hard. First, we already have collected $c$ coupons and because the first collision occurred at the $c$'th coupon, we have '$c-1$' distinct coupons. From here, we could proceed just like in the case of the original coupon collector problem. Therefore, we have,

$\begin{align}\mathbb{E}(K|C=c)&=c+\frac{n}{n-(c-1)}+\frac{n}{n-(c-2)}+\cdots+\frac{n}{1}\\&=c+n\text{ }H_{n-c+1} \\\end{align}$

where $H_n$ is the $n$'th Harmonic number. Note that the formula is valid for $1<c\leq n$.

We now deal with the (more challenging) 'reverse' question. That is, given we have completed the collection in $K=k$ coupons, how many cards were average would the first collision have occurred. That is, we look for $\mathbb{E}(C|K=k)$.

After trying out different approaches, I settled on the following approach. Let's first count the number of valid sequences. Let $N(n,j)$ be the number of sequences which collects $n$ different coupons exactly on the $j$'th draw. Obviously, $N(n,j)=0$ for $j<n$ and $N(n,n)=n!$.

With some thought, we could see that

$N(n,j)=n\text{ }N(n-1,j-1)+(n-1)\text{ }N(n,j-1)$

This works because, to collect $n$ exactly at the $j$'th draw, we could've (1) collected $n-1$ exactly at the $j-1$'th draw and the final coupon (for which we have $n$ choices) becomes the new coupon to complete the set (or) taken a sequence which completed the set of $n$ coupons in $j-1$ draws, choose any of the $n-1$ coupons except the last and insert it the penultimate position.

The $N$ function also has 'direct' formula via Stirling numbers of the second kind but we prefer this recurrence. We'll see why in a moment. Now, the expected value we are after can be written as,

$\displaystyle\mathbb{E}(C|K=k)=\frac{T(n,k)}{N(n,k)}$

where $T(n,j)$ is the sum of the first-collision-index of all the valid $N(n,j)$ sequences.

We don't have to look further to solve $T(n,j)$. Note that the neither of the two operations we described above for deriving the recurrence of $N(n,j)$ affects the first-collision-index of that sequence because the last coupon cannot be part of first collision (by definition, it has to be a new coupon to complete the set). Therefore, the same recurrence works for $T(n,j)$ as well!

The only tricky part is the boundary. Consider $T(n,n)$ which is the sum of all first collision indices where the set of $n$ coupons is collected in exactly $n$ draws. But those sequences must contain all distinct elements and therefore cannot have first collision. However, staying true to the process, we attribute a first collision index of $n$ to these sequences as we have collected $n$ coupons and the process have concluded. Therefore, we define, $T(n,n)=n \cdot n!$.

For example, using these recurrences, I get $\displaystyle\mathbb{E}(C|K=8,n=5)=612/175\approx 3.497$.

Now that we are on this subject, we can also check for $N(n,c,k)$, the number of sequences of coupons (from $n$ different coupons) that is 'complete' in the $k$'th draw and has its first collision on the $c$'th draw.

To compute $N(n,c,k)$, we must (1) choose '$c-1$' numbers in order (2) choose one among the '$c-1$' numbers for the $c$'th draw (3) choose one of the remaining '$n-c+1$' numbers to the last that completes the collection which means by now $c$ are already present in the sequence (4) choose $i$ of the remaining '$k-c-1$' positions for the '$n-c$' unseen numbers (5) partition the '$n-c$' numbers in $i$ non-empty subsets (6) use the '$c-1$' numbers at the start in the remaining positions without any constraint. Then,

$\displaystyle N(n,c,k)=\frac{n!}{(n-c+1)!}\cdot (c-1)\cdot(n-c+1)\cdot\sum_{i=0}^{k-c-1}\binom{k-c-1}{i}\cdot (n-c)!\left\{ \begin{matrix} i \\ n-c \end{matrix} \right\} \cdot (c-1)^{k-c-1-i}$

This can be simplified to finally give,

$\displaystyle N(n,c,k)=n!\cdot (c-1)\cdot\sum_{i=0}^{k-c-1}\binom{k-c-1}{i}\cdot \left\{ \begin{matrix} i \\ n-c \end{matrix} \right\} \cdot (c-1)^{k-c-1-i}$

I tried hard to find a recurrence or a generating function but failed miserably. I then delegated to CoPilot, Gemini and ChatGPT, and had a great time pointing one's error to the other. Finally, CoPilot (though I used it more than the others) figured out

$\displaystyle \sum_{n \geq 0,c > 0,k > c}N(n,c,k)\frac{x^n}{n!}\frac{y^{k-c-1}}{(k-c-1)!}z^{n-c}=\frac{x^2}{(1-xe^y)^2}\exp(y+xz(e^y-1))$

Gemini gave me the following python code that can be used to check the validity of the formula.

(code starts here..)

import itertools

def count_sequences(n, c, k):

    """

    n: total number of distinct digits available

    c: position of the first repetition (1-indexed)

    k: total length of the sequence

    """

    digits = list(range(n))

    valid_count = 0

    

    # Brute force through all n^k possibilities

    for seq in itertools.product(digits, repeat=k):

        # 1. Find the first repetition

        first_repeat_pos = 0

        seen = set()

        for i, val in enumerate(seq):

            if val in seen:

                first_repeat_pos = i + 1

                break

            seen.add(val)

        

        if first_repeat_pos != c:

            continue

            

        # 2. Check if the k-th digit is the one that completes the n-set

        if len(set(seq)) == n and len(set(seq[:-1])) == n - 1:

            valid_count += 1

            

    return valid_count

# Specific Example

n, c, k = 6, 4, 9

result = count_sequences(n, c, k)

print(f"count_sequences(n={n}, c={c}, k={k}) = {result}")

(... and ends here.)

Hope you enjoyed the discussion.

Until then

Yours Aye

Me

A classic Pursuit problem

One the most famous and classic pursuit problem goes like this: A pursuer, initially at point $(0, b)$, is trying to catch a pursuee, initially at point $(a,0)$ and running along the $x$-axis in the positive direction, by always running directly towards him. If they both run at constant speed and the ratio of pursuer's speed to that of the pursuee is $1:k$ (with $k<1$), how far do the pursuer runs in the pursuit?

This problem is discussed in The American Mathematical Monthly (P3573, Vol. 40, No. 7, (Aug. - Sep., 1933)) and Pi Mu Epsilon Journal (Vol. 6, No. 9, Fall 1978, Pg. 545) among others (for example Pursuit Curve and Pursuit).

In both AMA and PME, the solution is obtained (albeit for a special case where the pursuee starts at the origin) by solving a 'challenging' differential equation which gives the exact equation of the path traversed by the pursuer.

However, if the distance traversed by the pursuer is all we want, do we still have to go through the trouble of solving such a complicated differential equation?

A pursuit curve with $k=4/5$ (Note that Geogebra is unable to render the full curve).

One of the key points to note in this problem is that the respective distance travelled by the participants in given interval of time is directly proportional to the ratio of the speeds.

We make an assumption that the pursuer ends up catching the pursuee (say at $C$ as shown above). While this may be obvious from the full solution, for our purposes, this will remain an assumption.

Consider the infinitesimal distance $\,ds$, say between $F$ and $H$ as shown below, covered by the pursuer in time $\,dt$. Drawing tangents at these points, we know that the pursuee was at position $J$ and $K$ respectively.

As the distance covered by them are in proportion to their speeds, we know that $|JK|=k\cdot|FH|=k\,ds$.

Let $l$ denote the relative distance between the pursuer and pursuee. Then, $\,dl=|FJ|-|HK|$. To calculate this, we complete the right $\triangle FMH$ and draw a line passing through $J$ and perpendicular to $FJ$ such that it intersects line $HK$ at $L$.

We now rotate line $FJ$ about $H$ such that $F$ moves to $F'$ (not shown in image) in line $HK$. Also, because the rotation is infinitesimally small, $J$ moves to $L$ in line $HK$.

Then, $|FJ|\approx |F'L|=|FH|+|HK|-|LK|$. That is, $|FH|=|FJ|-|HK|+|LK| \implies \,ds=\,dl+|LK|$

Because both $\angle FJL$ and $\angle HLJ$ are approximately right angled, $\angle HFM = \angle LJK$. Therefore, by the similarity of $\triangle FMH$ and $\triangle JLK$ and the fact that $|JK|=k\,ds$, we have,

$|LK|=k|HM|=k\,dx$

where $\,dx$ is the infinitesimal movement of pursuer along the $x$-direction.

Combining everything we have so far, we see that $\,ds=\,dl+k\,dx$.

We have now shown that $\,d(s-l-kx)=0$ (or) $s - l - kx = \text{const.}$

We can choose $s$ to be the distance to be travelled by the pursuer to reach $C$ and $x$ to be the distance along the $x$-axis between the pursuer and $C$.

We know that at point $C$, $s=l=x=0$ which means that the constant in the equation above is $0$ giving us the relation that $s=l+kx$.

Before we celebrate our success, we notice a glaring inconvenience in the above relation. While we will know that $l=l_0$ is the relative distance between the two parties at the start of the pursuit and $s=s_0$ is the distance we are trying to find, we do not know $x_0$. However, this is simple to resolve.

Let $x_0=h_0+p_0$ where $h_0$ is the horizontal distance between the pursuer and the pursuee at the start (which should be known for the problem to be defined) and $p_0$ is the distance between the pursuee and $C$.

By our definition of point $C$, we know that the pursuer travels distance $s_0$ in the time the pursuee travels distance $p_0$. Therefore, $p_0=k\cdot s_0$. Using this and solving for $s_0$, we get (after dropping subscripts),

$\displaystyle s=\frac{l+k\cdot h}{1-k^2}$

which solves the problem of finding the distance traversed by the pursuer in a more general setting where the they are not necessarily in the same vertical.

Hope you enjoyed this post.

Until then

Yours Aye

Me

Arclength of the Cycloid(s) from velocity vectors

Cycloid, the curve traced by a point on the rim of a circle rolling on a straight line, probably needs no introduction. Amongst its other distinctions, it remains one of the first few curves to be rectified. In this post we attempt to rectify the cycloid(s) with kinematics.

The cycloid's trajectory is a combination of simple motions - rolling and translation. Because the circle rolls without slipping on the straight line, its translational velocity is equal to its rotational velocity.

In fact, at the point at which the circle touches the base, the velocities are exactly equal in magnitude and opposite in direction making this point the instantaneous center of rotation.

The green vectors in the image below shows the translational and rotational velocities, and the velocity of the tracing point is represented by the red vector. Further, the rotational velocity vector is further resolved along the perpendicular components $PG$ and $PO'$.

From the geometry of the configuration, we can easily see that $|PF|=2\cdot |HE|$.

Because $PE$ represents the circle's infinitesimal arclength (in time $dt$) and $HE$ is the component of that arclength along $PG$, we also know that $|HE|=d(|PG|)$. Therefore, the infinitesimal arclength $ds$ of the cycloid traversed by the tracing point $P$ in time $dt$ is,

$ds=|PF|=2\cdot |HE|=2\cdot d(|PG|)$

As the change in arclength is equal to twice the change in $PG$ in any given time interval, the total arclength of the semiarch of the cycloid is equal to twice the diameter of the circle.

We can use the almost same reasoning for hypocycloids which are curves traced by a point on a rolling  circle of radius (say) $b$ inside within a larger circle of radius (say) $a$. For integer $q=a/b$, we get nice closed hypocycloids. The case for $q=5$ is shown below.

A pentoid ($q=5$)

Note that because the circle rolls without slipping, the angular velocity of rotation of the rolling circle is $q$ times its angular velocity of its revolution. Consequently, the movement of the tracing point is a composition of (i) rotation of vector of length $a-b$ by angle $t$ in the counter-clockwise direction about the origin and (ii) rotation of a vector of length $b$ by an angle $(q-1)t$ in the clockwise direction about the origin.

The above in turn means that the velocity vectors of both the rotations noted above are the same and are proportional to $(a-b)$ times $\,dt$ (or equivalently $b$ times $(q-1)\,dt$) denoted by green vectors below.

Here again, from the geometry of the configuration, it is immediately clear that the infinitesimal arclength traversed by the tracing point in time $dt$ is $ds=|PF|=2|HE|$.

But, because the smaller circle rotates $q$ times faster, the magnitude of the velocity vector of any point on its circumference will be proportional to $b\cdot q\,dt=a\,dt$ which is represented by the cyan vector in the image above.

Therefore, $\displaystyle \frac{|HE|}{|H'E'|}=\frac{a-b}{a} \implies |HE|=\frac{a-b}{a}|H'E'|=\frac{a-b}{a}\,d(|PG|)$

We then finally have $\displaystyle ds=2\frac{a-b}{a}d(|PG|)$

Then, the arclength to be traversed by the tracing point to reach the top of the arch is $2(a-b)/a$ times its distance to the topmost point of the rolling circle.

While the demonstration given for the cycloid in Wikipedia is even nicer method to determine the arclength, it is always fun to arrive at the same results with different methods. As shown below, its is not very difficult to extend the demonstration quoted above to hypocycloids as well.

Hope you enjoyed this post.

Until then

Yours Aye

Me

A short note on the Median ratio property in Spherical Geometry

We know that in a plane triangle, the median is divided by the centroid in the ratio $2:1$. In this post, we address the question of how this relation generalizes in Spherical geometry.

Consider a spherical triangle $ABC$ with $D$, $E$ and $F$ being the midpoints of $BC$, $CA$ and $AB$. We know from the spherical analogy of Ceva's theorem, that the three medians $AD$, $BE$ and $CF$ meet at a common point, say $G$.

Using the volume and sine-of-side relation, we know that $|ABG|=|BCG|=|CAG|$.

Then, using our results from this post, we have,

$\displaystyle \frac{\sin AG}{\sin GD}=\frac{\sin BC}{\sin BD}$

Using the fact that $D$ is the midpoint of the arc $BC$,

$\displaystyle \frac{\sin AG}{\sin GD}=\frac{\sin BC}{\sin BD}=\frac{2\sin BD \cos BD}{\sin BD}=2\cos(BC/2)$

This shows that, unlike the planar case, the median ratio is dependent on the side on which the median rests ($BC$ in this case). However, by the same relation, we see that it is independent of the vertex of the median. That is, no matter the position of $A$ and as long as the $BC$ is fixed, the ratio remains invariant.

Hope you enjoyed the post.

Yours Aye

Me

A nice analogy in Spherical Geometry with a nicer application to Ceva

It is well known that the Law of sines has a spherical counterpart. But from the Wiki article, we can infer something even better. That is,

$\displaystyle V=\frac{1}{6}\sin{a}\sin{b}\sin{C}=\frac{1}{6}\sin{b}\sin{c}\sin{A}=\frac{1}{6}\sin{c}\sin{a}\sin{B}$

where $V$ is the volume of tetrahedron formed by the vertices of the spherical triangle along with the origin. The above can be re-written as

$\displaystyle V=\frac{1}{6}\sin{a}\sin{h_a}=\frac{1}{6}\sin{b}\sin{h_b}=\frac{1}{6}\sin{c}\sin{h_c}$

where $h_a$ is the altitude passing through $A$ and so on. This way, we can see the volume described above plays a role very similar to that of area in plane geometry. From this, we note that for two spherical triangles with the same height,

$V_1/V_2=\sin{b_1}/\sin{b_2}$

This fact can be used to prove Ceva's theorem, a powerful theorem in plane geometry. Though it is well known that Ceva's theorem also works for spherical geometry, I was not able to find a clear proof which we address in post.

Note that even though the volume and sine-of-side gives a nice analogy with the planar case, we lose the ability to add or subtract the sides and area which is very useful in the planar case. However, this does not seem to be a serious limitation.

While the proof given in Wiki article uses 'subtracting two areas', we can achieve the result even without that. Using the notations used in the Wiki article (but assuming a spherical triangle) and denoting the volume of a tetrahedron formed by points $X$, $Y$, $Z$ and the origin by $|XYZ|$,

$\displaystyle \frac{\sin OC}{\sin OF}=\frac{|OCB|}{|OBF|}$ and $\displaystyle \frac{\sin FB}{\sin AB}=\frac{|OFB|}{|OAB|}$

Combining the two, we have $\displaystyle \frac{\sin OC}{\sin OF}=\frac{|OBC|}{|OAB|}\frac{\sin AB}{\sin FB}$

We can use a similar argument to show $\displaystyle \frac{\sin OC}{\sin OF}=\frac{|OCA|}{|OAB|}\frac{\sin AB}{\sin AF}$

Using the above two, we have $\displaystyle \frac{|OBC|}{|OCA|}=\frac{\sin FB}{\sin AF}$

We can replicate this to the other sides as well. Then,

$\displaystyle \frac{\sin AF}{\sin FB} \cdot \frac{\sin BD}{\sin DC} \cdot \frac{\sin CE}{\sin EA}=\frac{|OCA|}{|OBC|} \cdot \frac{|OAB|}{|OCA|}\cdot \frac{|OBC|}{|OAB|}=1$

The existence of Ceva's theorem can then be used to show the concurrency of medians, altitudes and angle bisectors of spherical triangles (strictly speaking, we need the converse of Ceva's, but at this point I'm taking it for granted).

Hope you enjoyed this post.

Yours Aye

Me