Mathematics · Book 4 · Bachelor Year 2

University Mathematics — Year 2

University Mathematics — Year 2 · Bachelor Year 2

18Curves

Geometry now turns differential. A curve is a point moving through space; calculus gives us its velocity and acceleration, and geometry asks what is intrinsic — independent of how fast we traverse the trajectory. The answers are arc length, which measures the trajectory itself, and curvature, which measures how it bends. In dimension 33 a second invariant, torsion, measures how the curve twists out of its plane. The bookkeeping device for all of this is the moving Frenet frame.

18.1 Parametrized arcs

Definition 18.1 (Parametrized arc)

A parametrized arc of class Ck\mathcal{C}^k (k1k \geq 1) is a map γ ⁣:IRn\gamma \colon I \to \R^n of class Ck\mathcal{C}^k on an interval II. A point γ(t)\gamma(t) is regular if γ(t)0\gamma'(t) \neq 0, and the arc is regular if all its points are. The line through γ(t)\gamma(t) directed by γ(t)\gamma'(t) is the tangent line at a regular point.

Definition 18.2 (Change of parameter)

A change of parameter of class Ck\mathcal{C}^k is a Ck\mathcal{C}^k diffeomorphism θ ⁣:JI\theta \colon J \to I between intervals (θ0\theta' \neq 0 everywhere). The arcs γ\gamma and γθ\gamma \circ \theta are called equivalent; a geometric arc (or curve) is an equivalence class. Notions invariant under change of parameter — the trajectory, the tangent line, arc length, curvature — are called geometric.

Remark 18.3

The trajectory alone does not determine the geometric arc: the parametrizations t(cost,sint)t \mapsto (\cos t, \sin t) on [0,2π][0, 2\pi] and on [0,4π][0, 4\pi] have the same image but travel the circle once and twice. A geometric arc remembers the multiplicity and orientation of the traversal, not its speed.

Example 18.4 (Speed changes nothing geometric)

Parametrize the unit circle by

γ(t)=(cost2, sint2),t[0,2π].\gamma(t) = (\cos t^2,\ \sin t^2), \qquad t \in \intcc0{\sqrt{2\pi}} .

The speed γ(t)=2t\norm{\gamma'(t)} = 2t grows linearly, yet

L=02π2t ⁣dt=2π,L = \int_0^{\sqrt{2\pi}}2t\,\dd t = 2\pi ,

the same length as at constant speed — as Theorem 18.7 promises, via the change of parameter tt2t \mapsto t^2. The tangent line, the curvature computed from Proposition 18.17, and every other geometric quantity agree as well; only t=0t = 0 deserves a glance, where γ(0)=0\gamma'(0) = 0 makes this parametrization irregular although the trajectory is a perfect circle. Geometric statements tolerate bad parametrizations badly: reparametrize first, conclude second.

Example 18.5

The arc γ(t)=(t2,t3)\gamma(t) = (t^2, t^3) is C\mathcal{C}^\infty but not regular: γ(0)=(0,0)\gamma'(0) = (0, 0). Its trajectory, the semicubical parabola y2=x3y^2 = x^3, has a cusp at the origin: smoothness of the parametrization does not prevent a geometric singularity where the velocity vanishes. This is why the regularity hypothesis γ0\gamma' \neq 0 is not cosmetic.

18.2 Arc length

Definition 18.6 (Arc length)

Let γ ⁣:[a,b]Rn\gamma \colon [a, b] \to \R^n be a C1\mathcal{C}^1 arc. Its length is

L(γ)=abγ(t) ⁣dt,L(\gamma) = \int_a^b \norm{\gamma'(t)}\, \dd t ,

where \norm{\cdot} is the Euclidean norm. The arc-length function based at t0t_0 is s(t)=t0tγ(u) ⁣dus(t) = \int_{t_0}^t \norm{\gamma'(u)}\,\dd u.

Theorem 18.7 (Length is geometric; polygonal characterization)

  1. If θ ⁣:[c,d][a,b]\theta \colon [c, d] \to [a, b] is a C1\mathcal{C}^1 change of parameter, then L(γθ)=L(γ)L(\gamma \circ \theta) = L(\gamma).
  2. L(γ)L(\gamma) is the supremum of the lengths of inscribed polygons:

    L(γ)=sup{i=1mγ(ti)γ(ti1)  :  a=t0<t1<<tm=b}.L(\gamma) = \sup\Bigl\{\, \sum_{i=1}^{m} \norm{\gamma(t_i) - \gamma(t_{i-1})} \;:\; a = t_0 < t_1 < \dots < t_m = b \,\Bigr\}.

Proof. 1. By the change of variables t=θ(u)t = \theta(u) (Year 1 volume, valid since θ\theta is C1\mathcal{C}^1 monotone),

cd(γθ)(u) ⁣du=cdγ(θ(u))θ(u) ⁣du=abγ(t) ⁣dt,\int_c^d \norm{(\gamma\circ\theta)'(u)}\,\dd u = \int_c^d \norm{\gamma'(\theta(u))}\,\abs{\theta'(u)}\,\dd u = \int_a^b \norm{\gamma'(t)}\,\dd t ,

where we used (γθ)=θ(γθ)(\gamma\circ\theta)' = \theta'\cdot (\gamma'\circ\theta) and, if θ\theta is decreasing, the sign of θ\theta' is absorbed by the reversal of the bounds.

2. For any subdivision, γ(ti)γ(ti1)=ti1tiγ(t) ⁣dt\gamma(t_i) - \gamma(t_{i-1}) = \int_{t_{i-1}}^{t_i}\gamma'(t)\,\dd t, so by the triangle inequality for integrals γ(ti)γ(ti1)ti1tiγ\norm{\gamma(t_i) - \gamma(t_{i-1})} \leq \int_{t_{i-1}}^{t_i}\norm{\gamma'}: every polygon is shorter than L(γ)L(\gamma), hence supL(γ)\sup \leq L(\gamma).

For the reverse inequality, let ε>0\varepsilon > 0. Since γ\gamma' is continuous on the compact [a,b][a,b], it is uniformly continuous: there is δ>0\delta > 0 with γ(t)γ(u)ε\norm{\gamma'(t) - \gamma'(u)} \leq \varepsilon whenever tuδ\abs{t - u} \leq \delta. Take a subdivision of step δ\leq \delta. On each piece, for t[ti1,ti]t \in [t_{i-1}, t_i],

γ(ti)γ(ti1)=ti1tiγ(t) ⁣dt=(titi1)γ(ti1)+Ri,Riε(titi1),\gamma(t_i) - \gamma(t_{i-1}) = \int_{t_{i-1}}^{t_i}\gamma'(t)\,\dd t = (t_i - t_{i-1})\,\gamma'(t_{i-1}) + R_i, \qquad \norm{R_i} \leq \varepsilon\,(t_i - t_{i-1}),

since Ri=ti1ti(γ(t)γ(ti1)) ⁣dtR_i = \int_{t_{i-1}}^{t_i}(\gamma'(t) - \gamma'(t_{i-1}))\,\dd t. Hence

γ(ti)γ(ti1)(titi1)γ(ti1)ε(titi1).\norm{\gamma(t_i) - \gamma(t_{i-1})} \geq (t_i - t_{i-1})\norm{\gamma'(t_{i-1})} - \varepsilon (t_i - t_{i-1}) .

Summing, and comparing (titi1)γ(ti1)\sum (t_i - t_{i-1})\norm{\gamma'(t_{i-1})} with abγ\int_a^b\norm{\gamma'} (a Riemann sum of the continuous function γ\norm{\gamma'}, within ε(ba)\varepsilon(b - a) of the integral for δ\delta small enough by uniform continuity again), we get a polygon of length L(γ)2ε(ba)\geq L(\gamma) - 2\varepsilon(b - a). Letting ε0\varepsilon \to 0 proves the claim.

Example 18.8 (Archimedes and the inscribed polygons)

For the unit circle, the inscribed regular nn-gon has length Ln=2nsinπnL_n = 2n\sin\frac\pi n, and the expansion sinx=xx36+O(x5)\sin x = x - \frac{x^3}6 + O(x^5) gives

Ln=2ππ33n2+O(1n4):L_n = 2\pi - \frac{\pi^3}{3n^2} + O\Bigl(\frac1{n^4}\Bigr) :

the polygonal lengths of Theorem 18.7 converge quadratically. Numerically: L6=6L_6 = 6 (the hexagon, giving the crude π>3\pi > 3), while L96=192sinπ966.28206L_{96} = 192\sin\frac{\pi}{96} \approx 6.28206 against 2π6.283192\pi \approx 6.28319 — the error 0.001130.00113 agrees with the predicted π3/(3962)0.00112\pi^3/(3\cdot96^2) \approx 0.00112. This is why Archimedes, doubling the hexagon five times to 9696 sides, could bracket π\pi to three digits by hand: each doubling divides the error by four. The supremum in the polygonal characterization is not just attained in the limit; it is attained fast, because a smooth curve separates from its chords only at second order.

Theorem 18.9 (Arc-length parametrization)

Let γ ⁣:IRn\gamma \colon I \to \R^n be a regular Ck\mathcal{C}^k arc (k1k \geq 1). The arc-length function ss is a Ck\mathcal{C}^k diffeomorphism from II onto an interval JJ, and γ~=γs1\tilde\gamma = \gamma \circ s^{-1} satisfies γ~=1\norm{\tilde\gamma'} = 1 everywhere. Up to translation of the parameter and orientation, this arc-length (or unit-speed) parametrization is unique.

Proof. s(t)=γ(t)>0s'(t) = \norm{\gamma'(t)} > 0 and ss' is Ck1\mathcal{C}^{k-1} (composition of the Ck1\mathcal{C}^{k-1} map γ\gamma' with the norm, smooth away from 00), so ss is Ck\mathcal{C}^k, strictly increasing, a bijection onto J=s(I)J = s(I), and its inverse is Ck\mathcal{C}^k by the inverse function theorem in dimension 11 (Year 1 volume). Then

γ~(σ)=γ(t)s(t)=γ(t)γ(t),t=s1(σ),\tilde\gamma'(\sigma) = \frac{\gamma'(t)}{s'(t)} = \frac{\gamma'(t)}{\norm{\gamma'(t)}}, \qquad t = s^{-1}(\sigma),

a unit vector. If γ^=γθ\hat\gamma = \gamma\circ\theta is another unit-speed parametrization, then θ=1\abs{\theta'} = 1, so θ=±1\theta' = \pm 1 constant (continuity), i.e. θ(u)=±u+c\theta(u) = \pm u + c.

Remark 18.10

Arc length is the parameter that separates geometry from dynamics. A trajectory can be traversed with any speed profile — the physics of the motion — but every parametrization-invariant question (shape, bending, osculation) has a canonical clock, the distance travelled. This is why all curvature formulas below are defined at unit speed and then translated to arbitrary parametrizations by the chain rule: the translation factors are powers of v=sv = s', and tracking them correctly is the entire content of Proposition 18.17.

Example 18.11 (Circle and helix)

For the circle γ(t)=(Rcost,Rsint)\gamma(t) = (R\cos t, R\sin t), γ=R\norm{\gamma'} = R, so s=Rts = Rt and the length of a full turn is 2πR2\pi R. For the helix γ(t)=(acost, asint, bt)\gamma(t) = (a\cos t,\ a\sin t,\ bt) with a>0a > 0, γ(t)=a2+b2\norm{\gamma'(t)} = \sqrt{a^2 + b^2} is constant: the helix is traversed at constant speed, and s=ta2+b2s = t\sqrt{a^2 + b^2}.

Example 18.12 (Arc length in polar coordinates)

A polar curve r=r(θ)r = r(\theta) is the arc γ(θ)=(rcosθ, rsinθ)\gamma(\theta) = (r\cos\theta,\ r\sin\theta), with

γ(θ)=(rcosθrsinθ, rsinθ+rcosθ),γ(θ)2=r2+r2\gamma'(\theta) = (r'\cos\theta - r\sin\theta,\ r'\sin\theta + r\cos\theta), \qquad \norm{\gamma'(\theta)}^2 = r'^2 + r^2

(the cross terms cancel): the polar length element is r2+r2 ⁣dθ\sqrt{r^2 + r'^2}\,\dd\theta. For the cardioid r=1+cosθr = 1 + \cos\theta:

r2+r2=(1+cosθ)2+sin2θ=2+2cosθ=4cos2θ2,r^2 + r'^2 = (1 + \cos\theta)^2 + \sin^2\theta = 2 + 2\cos\theta = 4\cos^2\tfrac\theta2,

and on [π,π]\intcc{-\pi}{\pi} the half-angle cosine is nonnegative, so

L=ππ2cosθ2 ⁣dθ=[4sinθ2]ππ=4(4)=8:L = \int_{-\pi}^{\pi}2\cos\tfrac\theta2\,\dd\theta = \Bigl[4\sin\tfrac\theta2\Bigr]_{-\pi}^{\pi} = 4 - (-4) = 8 :

like the cycloid arch of Exercise 18.1, a curve built from circles has a rational length, with no π\pi anywhere. The half-angle factorization is the standard trick for lengths of circle-generated curves; when it fails (the ellipse), the length is a genuinely new function — an elliptic integral, beyond elementary closed forms.

18.3 Curvature in the plane

Throughout this section, arcs are C2\mathcal{C}^2 and regular in the oriented Euclidean plane. We parametrize by arc length and write T(s)=γ~(s)T(s) = \tilde\gamma'(s) for the unit tangent, and N(s)N(s) for the unit vector directly orthogonal to T(s)T(s) (rotation of TT by +π/2+\pi/2).

Theorem 18.13 (Plane Frenet formulas)

Let γ~\tilde\gamma be a unit-speed C2\mathcal{C}^2 arc in the oriented plane. There is a continuous function κ\kappa, the (algebraic) curvature, such that

T(s)=κ(s)N(s),N(s)=κ(s)T(s).T'(s) = \kappa(s)\, N(s), \qquad N'(s) = -\kappa(s)\, T(s).

Proof. Since T(s)2=1\norm{T(s)}^2 = 1 for all ss, differentiating the scalar product gives 2T(s),T(s)=02\langle T'(s), T(s)\rangle = 0: T(s)T'(s) is orthogonal to T(s)T(s), hence collinear with N(s)N(s) (dimension 22); write T(s)=κ(s)N(s)T'(s) = \kappa(s) N(s) with κ(s)=T(s),N(s)\kappa(s) = \langle T'(s), N(s)\rangle, continuous. Likewise NNN' \perp N, so N=λTN' = \lambda T; and differentiating T,N=0\langle T, N\rangle = 0 gives T,N+T,N=κ+λ=0\langle T', N\rangle + \langle T, N'\rangle = \kappa + \lambda = 0.

Definition 18.14

When κ(s)0\kappa(s) \neq 0, the radius of curvature is R(s)=1/κ(s)R(s) = 1/\abs{\kappa(s)} and the center of curvature is γ~(s)+1κ(s)N(s)\tilde\gamma(s) + \frac{1}{\kappa(s)} N(s); the circle with that center and radius R(s)R(s) is the osculating circle, the best circular approximation of the curve at γ~(s)\tilde\gamma(s).

Example 18.15 (The osculating circle of the exponential)

For y=exy = \eu^x at the point (0,1)(0, 1): f(0)=f(0)=1f'(0) = f''(0) = 1, so by the graph formula below,

κ(0)=1(1+1)3/2=122,R=22.\kappa(0) = \frac{1}{(1 + 1)^{3/2}} = \frac1{2\sqrt2}, \qquad R = 2\sqrt2 .

The unit tangent is T=(1,1)2T = \frac{(1, 1)}{\sqrt2}, the direct normal N=(1,1)2N = \frac{(-1, 1)}{\sqrt2}, and the center of curvature is

(0,1)+22(1,1)2=(2, 3):(0, 1) + 2\sqrt2\cdot\frac{(-1, 1)}{\sqrt2} = (-2,\ 3) :

the osculating circle has equation (x+2)2+(y3)2=8(x + 2)^2 + (y - 3)^2 = 8. As a check of the “best circular approximation” claim: solving the circle’s equation for yy near (0,1)(0,1) and expanding gives y=1+x+x22+O(x3)y = 1 + x + \frac{x^2}2 + O(x^3) — exactly the second-order Taylor expansion of ex\eu^x. The osculating circle matches value, slope and second derivative; an ordinary tangent circle would match only the first two.

Example 18.16 (The evolute of a circle is its center)

For the circle of radius RR traversed counterclockwise, κ=1/R\kappa = 1/R and NN points toward the center, so the center of curvature γ~+1κN\tilde\gamma + \frac1\kappa N is the center of the circle, for every ss: the osculating circle of a circle is the circle itself, and the locus of centers of curvature collapses to a point. This degenerate case calibrates Exercise 18.6: there the evolute’s velocity is κκ2N-\frac{\kappa'}{\kappa^2}N, which vanishes identically precisely when κ\kappa is constant.

Proposition 18.17 (Curvature in an arbitrary parametrization)

For a regular C2\mathcal{C}^2 plane arc γ(t)=(x(t),y(t))\gamma(t) = (x(t), y(t)),

κ(t)=x(t)y(t)y(t)x(t)(x(t)2+y(t)2)3/2,\kappa(t) = \frac{x'(t)\,y''(t) - y'(t)\,x''(t)} {\bigl(x'(t)^2 + y'(t)^2\bigr)^{3/2}} ,

in particular κ=y(1+y2)3/2\kappa = \dfrac{y''}{(1 + y'^2)^{3/2}} for a graph y=f(x)y = f(x).

Proof. Write v(t)=γ(t)=s(t)v(t) = \norm{\gamma'(t)} = s'(t), so γ=vT\gamma' = vT (composing the unit-speed data with ss). Differentiating,

γ=vT+vTs=vT+v2κN.\gamma'' = v'T + v\,T'\cdot s' = v'T + v^2\kappa N .

Now take the determinant (in the canonical oriented basis) of (γ,γ)(\gamma', \gamma''): since det(T,T)=0\det(T, T) = 0 and det(T,N)=1\det(T, N) = 1,

det(γ,γ)=det(vT, vT+v2κN)=v3κ.\det(\gamma', \gamma'') = \det(vT,\ v'T + v^2\kappa N) = v^3\kappa .

The left side is xyyxx'y'' - y'x'', and v3=(x2+y2)3/2v^3 = (x'^2 + y'^2)^{3/2}. The graph case is the parametrization t(t,f(t))t \mapsto (t, f(t)).

Example 18.18 (Circle, line, parabola)

A line has κ=0\kappa = 0 (and conversely: T=0T' = 0 means TT constant, so γ~(s)=γ~(0)+sT\tilde\gamma(s) = \tilde\gamma(0) + sT, a line). The circle of radius RR traversed counterclockwise has κ=1/R\kappa = 1/R: with γ(t)=(Rcost,Rsint)\gamma(t) = (R\cos t, R\sin t), the formula gives κ=R2/R3\kappa = R^2/R^3. For the parabola y=x2/2y = x^2/2: κ(x)=1/(1+x2)3/2\kappa(x) = 1/(1 + x^2)^{3/2}, maximal at the vertex — the parabola is most sharply bent where it turns around.

Remark 18.19 (Common pitfalls around curvature)

(i) The algebraic curvature of a plane arc changes sign when the orientation of the arc or of the plane is reversed: only κ\abs\kappa and R=1/κR = 1/\abs\kappa are purely geometric. A circle traversed clockwise has κ=1/R\kappa = -1/R. (ii) The graph formula κ=f/(1+f2)3/2\kappa = f''/(1 + f'^2)^{3/2} silently chooses the parametrization by xx; applying it to a curve that is not a graph near the point (vertical tangent) is the classical blunder. (iii) At a point where γ=0\gamma' = 0 nothing is defined — neither TT nor κ\kappa — and the trajectory may genuinely break (Example 18.5); always check regularity before differentiating the unit tangent. (iv) In space, κ=T0\kappa = \norm{T'} \geq 0 by convention: there is no sign to get wrong, but also no sign to exploit — inflection-type information moves into the torsion. (v) Finally, κ\kappa is a derivative with respect to arc length: for a non-unit-speed parametrization, forgetting the factor v3v^3 in Proposition 18.17 is the most frequent error in practice.

The parabola y = x2/2, its moving Frenet frame (T, N), and the osculating circle at the vertex (radius 1, since (0) = 1). The frame turns as the point moves; curvature is the rate of that turning per unit of arc length.
Figure 18.1. The parabola y=x2/2y = x^2/2, its moving Frenet frame (T,N)(T, N), and the osculating circle at the vertex (radius 11, since κ(0)=1\kappa(0) = 1). The frame turns as the point moves; curvature is the rate of that turning per unit of arc length.

Theorem 18.20 (Curvature determines the curve)

Let κ ⁣:JR\kappa \colon J \to \R be continuous. There exists a unit-speed C2\mathcal{C}^2 arc in the plane with curvature κ\kappa, and it is unique up to a direct isometry (rotation followed by translation).

Proof. Existence. Fix s0Js_0 \in J and set φ(s)=s0sκ(u) ⁣du\varphi(s) = \int_{s_0}^s \kappa(u)\,\dd u, then

γ~(s)=(s0scosφ(u) ⁣du, s0ssinφ(u) ⁣du).\tilde\gamma(s) = \Bigl(\int_{s_0}^s \cos\varphi(u)\,\dd u,\ \int_{s_0}^s \sin\varphi(u)\,\dd u\Bigr).

Then T(s)=(cosφ(s),sinφ(s))T(s) = (\cos\varphi(s), \sin\varphi(s)) is a unit vector, N(s)=(sinφ,cosφ)N(s) = (-\sin\varphi, \cos\varphi), and

T(s)=φ(s)(sinφ,cosφ)=κ(s)N(s):T'(s) = \varphi'(s)\,(-\sin\varphi, \cos\varphi) = \kappa(s)\,N(s) :

the arc is unit-speed with curvature κ\kappa.

Uniqueness. Let γ1,γ2\gamma_1, \gamma_2 be unit-speed arcs with the same curvature. Each unit tangent lifts to an angle function, by an explicit construction: view TjT_j as the complex number zj=aj+ibjz_j = a_j + \iu b_j of modulus 11, pick φj(0)\varphi_j(0) with zj(0)=eiφj(0)z_j(0) = \eu^{\iu\varphi_j(0)}, and set

φj(s)=φj(0)+0sdet(Tj,Tj)(u) ⁣du.\varphi_j(s) = \varphi_j(0) + \int_0^s\det\bigl(T_j, T_j'\bigr)(u)\,\dd u .

From zj=1\abs{z_j} = 1: Re(zjzj)=0\operatorname{Re}(\conj{z_j}\,z_j') = 0, so zjzj=idet(Tj,Tj)=iφj\conj{z_j}\,z_j' = \iu\det(T_j, T_j') = \iu\,\varphi_j', i.e. zj=iφjzjz_j' = \iu\varphi_j'z_j; then

(zjeiφj)=eiφj(zjiφjzj)=0,\bigl(z_j\,\eu^{-\iu\varphi_j}\bigr)' = \eu^{-\iu\varphi_j}\bigl(z_j' - \iu\varphi_j'z_j\bigr) = 0,

so zj=eiφjz_j = \eu^{\iu\varphi_j} throughout: Tj=(cosφj,sinφj)T_j = (\cos\varphi_j, \sin\varphi_j) with φj\varphi_j of class C1\mathcal C^1. Moreover det(Tj,Tj)=det(Tj,κNj)=κ\det(T_j, T_j') = \det(T_j, \kappa N_j) = \kappa, so φj=κ\varphi_j' = \kappa. Hence φ2=φ1+c\varphi_2 = \varphi_1 + c for a constant cc: T2T_2 is T1T_1 rotated by the fixed angle cc, so integrating, γ2=ρ(γ1)+w\gamma_2 = \rho(\gamma_1) + w where ρ\rho is the rotation of angle cc and ww a constant vector.

Remark 18.21

This is the one-dimensional prototype of a fundamental theorem of geometry: a complete set of local invariants (here, one function) classifies the object up to rigid motion. The three-dimensional version below needs two invariants.

Example 18.22 (Constant curvature means circle)

Take κκ0>0\kappa \equiv \kappa_0 > 0 in the existence formula: φ(s)=κ0s\varphi(s) = \kappa_0 s and

γ~(s)=(sinκ0sκ0, 1cosκ0sκ0):\tilde\gamma(s) = \Bigl(\frac{\sin\kappa_0s}{\kappa_0},\ \frac{1 - \cos\kappa_0s}{\kappa_0}\Bigr) :

the circle of radius 1/κ01/\kappa_0 centered at (0,1/κ0)(0, 1/\kappa_0), traversed at unit speed. By the uniqueness half of the theorem, every unit-speed arc of constant curvature κ0\kappa_0 is a piece of a circle of radius 1/κ01/\kappa_0 (or a line if κ0=0\kappa_0 = 0) — the converse of the computation in Example 18.18, and the plane case of Exercise 18.9.

Example 18.23 (Reconstructing a curve from its curvature)

Which unit-speed curve has radius of curvature R(s)=1+s2R(s) = 1 + s^2? Following the existence proof with κ(s)=11+s2\kappa(s) = \frac1{1+s^2} and s0=0s_0 = 0: φ(s)=arctans\varphi(s) = \arctan s, so

T(s)=(cosarctans, sinarctans)=(11+s2, s1+s2),T(s) = (\cos\arctan s,\ \sin\arctan s) = \Bigl(\frac{1}{\sqrt{1+s^2}},\ \frac{s}{\sqrt{1+s^2}}\Bigr),

and integrating,

γ~(s)=(ln(s+1+s2), 1+s21).\tilde\gamma(s) = \Bigl(\ln\bigl(s + \sqrt{1 + s^2}\bigr),\ \sqrt{1 + s^2} - 1\Bigr).

Setting x=ln(s+1+s2)x = \ln(s + \sqrt{1+s^2}), i.e. s=sinhxs = \sinh x, the second coordinate is coshx1\cosh x - 1: the curve is the catenary y=coshx1y = \cosh x - 1. This closes the loop with Exercise 18.3, where we computed R=cosh2x=1+sinh2x=1+s2R = \cosh^2 x = 1 + \sinh^2 x = 1 + s^2 directly: the fundamental theorem guarantees the catenary is the only curve with this curvature profile, up to a direct isometry.

Remark 18.24 (Where curvature is used next)

The decomposition γ=vT+v2κN\gamma'' = v'T + v^2\kappa N obtained in the proof of Proposition 18.17 is the kinematics of every curved motion: tangential versus centripetal acceleration. Curvature returns for surfaces (Chapter 19) through the curvature of curves drawn on them, and the envelope calculus of this chapter’s weekend problem — evolutes, caustics — is the geometric optics of wavefronts. The Year 3 volume takes the intrinsic point of view up again for submanifolds of Rn\R^n.

18.4 Frenet frame in space

Now let γ~ ⁣:JR3\tilde\gamma \colon J \to \R^3 be a unit-speed C3\mathcal{C}^3 arc that is biregular: T(s)0T'(s) \neq 0 for all ss. Then κ(s)=T(s)>0\kappa(s) = \norm{T'(s)} > 0 defines the curvature (no sign in space: there is no preferred normal orientation), and we set:

N(s)=T(s)κ(s)(principal normal),B(s)=T(s)N(s)(binormal),N(s) = \frac{T'(s)}{\kappa(s)} \quad\text{(principal normal)}, \qquad B(s) = T(s) \wedge N(s) \quad\text{(binormal)} ,

so that (T,N,B)(T, N, B) is a direct orthonormal frame, the Frenet frame. The plane through γ~(s)\tilde\gamma(s) spanned by T,NT, N is the osculating plane.

Theorem 18.25 (Frenet formulas in space)

There is a continuous function τ\tau, the torsion, with

T=κN,N=κT+τB,B=τN.T' = \kappa N, \qquad N' = -\kappa T + \tau B, \qquad B' = -\tau N .

Proof. The first formula is the definition of NN. Each of the vectors T,N,BT, N, B has constant norm 11 and they are pairwise orthogonal; differentiating the six relations X,Y=δXY\langle X, Y\rangle = \delta_{XY} shows the matrix of (T,N,B)(T', N', B') in the basis (T,N,B)(T, N, B) is antisymmetric: indeed X,Y+X,Y=0\langle X', Y\rangle + \langle X, Y'\rangle = 0 and X,X=0\langle X', X\rangle = 0. Its (N,T)(N, T) entry is N,T=N,T=κ\langle N', T\rangle = -\langle N, T'\rangle = -\kappa, and its (T,B)(T, B)-column entry T,B=κN,B=0\langle T', B\rangle = \kappa\langle N, B\rangle = 0. Naming the remaining free entry τ=N,B\tau = \langle N', B\rangle gives exactly the three displayed formulas: antisymmetry fills in B,N=τ\langle B', N\rangle = -\tau and B,T=0\langle B', T\rangle = 0. Continuity of τ=N,B\tau = \langle N', B\rangle is clear since NN' and BB are continuous (γ~\tilde\gamma is C3\mathcal{C}^3, so N=T/κN = T'/\kappa is C1\mathcal{C}^1).

Example 18.26 (The Darboux vector)

The three Frenet formulas compress into one. Set ω(s)=τT+κB\omega(s) = \tau\,T + \kappa\,B (the Darboux vector). Using BT=NB \wedge T = N, TN=BT \wedge N = B, NB=TN \wedge B = T:

ωT=κN=T,ωN=τBκT=N,ωB=τN=B:\omega \wedge T = \kappa\,N = T', \qquad \omega \wedge N = \tau\,B - \kappa\,T = N', \qquad \omega \wedge B = -\tau\,N = B' :

each frame vector evolves by X=ωXX' = \omega \wedge X, the kinematic signature of an instantaneous rotation with angular velocity vector ω\omega. The frame spins at rate ω=κ2+τ2\norm\omega = \sqrt{\kappa^2 + \tau^2} about the moving axis ω\omega; curvature is the component of the spin about the binormal, torsion the component about the tangent. For the helix, ω\omega is a constant vector along the cylinder’s axis — which is exactly why the helix’s frame precesses evenly. The antisymmetry of the Frenet matrix, exploited in Exercise 18.7, is the matrix form of this single geometric fact.

Proposition 18.27 (Torsion measures planarity)

A biregular arc is contained in a plane if and only if τ0\tau \equiv 0; in that case the plane is the (constant) osculating plane.

Proof. If τ0\tau \equiv 0, then B=0B' = 0, so BB is a constant unit vector B0B_0, and

 ⁣d ⁣dsγ~(s),B0=T(s),B0=0:\frac{\dd}{\dd s}\langle \tilde\gamma(s), B_0\rangle = \langle T(s), B_0\rangle = 0 :

γ~,B0\langle \tilde\gamma, B_0\rangle is constant, so the arc lies in a plane orthogonal to B0B_0. Conversely, if the arc lies in a plane PP, then TT and TT' (hence NN) are parallel to the direction of PP for all ss; so B=TNB = T \wedge N is one of the two unit normals of PP, and being continuous it is constant; then 0=B=τN0 = B' = -\tau N with N0N \neq 0 forces τ0\tau \equiv 0.

Example 18.28 (A tilted circle has zero torsion)

The arc γ(t)=(cost, sint2, sint2)\gamma(t) = \bigl(\cos t,\ \tfrac{\sin t}{\sqrt2},\ \tfrac{\sin t}{\sqrt2}\bigr) lies in the plane y=zy = z, and is the unit circle of that plane (check: γ(t)=1\norm{\gamma(t)} = 1 and the plane’s orthonormal basis (1,0,0)(1,0,0), (0,12,12)(0, \tfrac1{\sqrt2}, \tfrac1{\sqrt2}) exhibits the standard parametrization). Without any Frenet computation, Proposition 18.27 predicts τ0\tau \equiv 0, and the fixed binormal must be the plane’s unit normal ±(0,12,12)\pm(0, \tfrac1{\sqrt2}, -\tfrac1{\sqrt2}). Torsion does not measure being “tilted in space”; it measures leaving a plane. Only the helix’s nonzero bb below produces genuine torsion.

Example 18.29 (The helix)

For the helix γ(t)=(acost,asint,bt)\gamma(t) = (a\cos t, a\sin t, bt), a>0a > 0, we computed s=cts = ct with c=a2+b2c = \sqrt{a^2 + b^2}. Then

T=1c(asint, acost, b),T ⁣dt ⁣ds=1c2(acost,asint,0),T = \frac1c(-a\sin t,\ a\cos t,\ b), \qquad T' \cdot \frac{\dd t}{\dd s} = \frac{1}{c^2}(-a\cos t, -a\sin t, 0),

so κ=a/c2=a/(a2+b2)\kappa = a/c^2 = a/(a^2 + b^2) and N=(cost,sint,0)N = (-\cos t, -\sin t, 0): the principal normal points horizontally toward the axis. Next B=TN=1c(bsint,bcost,a)B = T \wedge N = \frac1c(b\sin t, -b\cos t, a), and B ⁣dt ⁣ds=bc2(cost,sint,0)=τNB' \frac{\dd t}{\dd s} = \frac{b}{c^2}(\cos t, \sin t, 0) = -\tau N gives

 κ=aa2+b2,τ=ba2+b2. \boxed{\ \kappa = \frac{a}{a^2 + b^2}, \qquad \tau = \frac{b}{a^2 + b^2}. \ }

Both invariants are constant — and one can show, conversely, that the only biregular curves with constant κ>0\kappa > 0 and constant τ\tau are helices (circles when τ=0\tau = 0). Note the signs: b>0b > 0 gives a right-handed helix with positive torsion.

Example 18.30 (Curvature and torsion without arc length; the twisted cubic)

Reparametrizing by arc length is usually impossible in closed form, so the invariants must be extracted from the raw derivatives. Write v=γ=sv = \norm{\gamma'} = s'; then γ=vT\gamma' = vT and, as in the proof of Proposition 18.17,

γ=vT+v2κN,γγ=v3κ(TN)=v3κB.\gamma'' = v'T + v^2\kappa N, \qquad \gamma' \wedge \gamma'' = v^3\kappa\,(T \wedge N) = v^3\kappa\,B .

Taking norms (κ0\kappa \geq 0 in space):

κ=γγv3.\kappa = \frac{\norm{\gamma' \wedge \gamma''}}{v^3} .

Differentiating γ\gamma'' once more and converting N=v(κT+τB)N' = v(-\kappa T + \tau B) (chain rule through ss), the only BB-component comes from the last term:

γ=(vv3κ2)T+(vvκ+(v2κ))N+v3κτB,\gamma''' = \bigl(v'' - v^3\kappa^2\bigr)T + \bigl(v'v\kappa + (v^2\kappa)'\bigr)N + v^3\kappa\tau\,B,

so that, pairing with γγ=v3κB\gamma' \wedge \gamma'' = v^3\kappa B,

det(γ,γ,γ)=γγ, γ=v6κ2τ,i.e.τ=det(γ,γ,γ)γγ2.\det(\gamma', \gamma'', \gamma''') = \langle\gamma' \wedge \gamma'',\ \gamma'''\rangle = v^6\kappa^2\tau, \qquad\text{i.e.}\qquad \tau = \frac{\det(\gamma', \gamma'', \gamma''')}{\norm{\gamma' \wedge \gamma''}^2} .

Application to the twisted cubic γ(t)=(t, t2, t3)\gamma(t) = (t,\ t^2,\ t^3) at t=0t = 0: γ=(1,0,0)\gamma' = (1, 0, 0), γ=(0,2,0)\gamma'' = (0, 2, 0), γ=(0,0,6)\gamma''' = (0, 0, 6), so v=1v = 1,

γγ=(0,0,2),κ(0)=2,det(γ,γ,γ)=12,τ(0)=124=3.\gamma' \wedge \gamma'' = (0, 0, 2), \qquad \kappa(0) = 2, \qquad \det(\gamma', \gamma'', \gamma''') = 12, \qquad \tau(0) = \frac{12}{4} = 3 .

Closing insight: both formulas are ratios in which the speed vv cancels to the exact degree needed — κ\kappa scales like a second derivative per unit length, τ\tau like the mixed volume of three derivatives per squared area — which is why they are geometric while γ\gamma'' itself is not.

Remark 18.31 (Fundamental theorem for space curves)

As in the plane, the pair (κ,τ)(\kappa, \tau) with κ>0\kappa > 0 determines a biregular arc up to direct isometry of R3\R^3: the Frenet formulas form a linear differential system for the frame (T,N,B)(T, N, B), to which the Cauchy–Lipschitz theory of Chapter 16 applies; orthonormality of the solution frame is preserved because the coefficient matrix is antisymmetric (same Gram-matrix argument as in Exercise 18.7), and the curve is recovered by integrating TT. We leave the details to the reader as a substantial but instructive exercise.

18.5 Local study: position with respect to the tangent

Proposition 18.32 (Local shape at a regular point)

Let γ\gamma be a plane arc of class Ck\mathcal{C}^k at t0t_0, with pp the smallest index with γ(p)(t0)0\gamma^{(p)}(t_0) \neq 0 and qq the smallest index >p> p with γ(q)(t0)\gamma^{(q)}(t_0) not collinear with γ(p)(t0)\gamma^{(p)}(t_0) (assuming both exist, qkq \leq k). In the basis (u,v)=(γ(p)(t0),γ(q)(t0))(u, v) = (\gamma^{(p)}(t_0), \gamma^{(q)}(t_0)) centered at γ(t0)\gamma(t_0), Taylor–Young gives coordinates

X(t)(tt0)pp!,Y(t)(tt0)qq!.X(t) \sim \frac{(t - t_0)^p}{p!}, \qquad Y(t) \sim \frac{(t - t_0)^q}{q!} .

The local picture depends only on the parities of pp and qq:

pp odd, qq evenordinary pointcurve crosses no line, stays on one side of tangent
pp odd, qq oddinflection pointcurve crosses its tangent
pp even, qq oddcusp of the first kindboth branches on opposite sides of tangent
pp even, qq evencusp of the second kindboth branches on the same side

Proof. Taylor–Young at order qq (the function γ\gamma is Cq\mathcal{C}^q near t0t_0):

γ(t)γ(t0)=j=pq(tt0)jj!γ(j)(t0)+o((tt0)q).\gamma(t) - \gamma(t_0) = \sum_{j=p}^{q} \frac{(t-t_0)^j}{j!}\,\gamma^{(j)}(t_0) + o\bigl((t-t_0)^q\bigr).

By the choice of pp and qq, each γ(j)(t0)\gamma^{(j)}(t_0) with pj<qp \leq j < q is collinear with uu; collecting components in the basis (u,v)(u, v): X(t)=(tt0)pp!(1+o(1))X(t) = \frac{(t-t_0)^p}{p!}(1 + o(1)) and Y(t)=(tt0)qq!(1+o(1))Y(t) = \frac{(t-t_0)^q}{q!}(1 + o(1)). The sign table of XX and YY for tt0t \gtrless t_0 — governed exactly by the parities — gives the four pictures: for instance if pp is even, X>0X > 0 on both sides (both branches leave in the direction +u+u: a cusp), and the side of the tangent line (signY\operatorname{sign} Y) flips with qq odd.

Example 18.33

For γ(t)=(t2,t3)\gamma(t) = (t^2, t^3) at t0=0t_0 = 0 (Example 18.5): γ(0)=(2,0)\gamma'' (0)= (2, 0), γ(0)=(0,6)\gamma'''(0) = (0, 6), so p=2p = 2, q=3q = 3: a cusp of the first kind, the familiar picture of the semicubical parabola. For γ(t)=(t,t3)\gamma(t) = (t, t^3) at 00: p=1p = 1, q=3q = 3: inflection — the cubic crosses its tangent.

Remark 18.34 (Perspectives within this volume)

Curves feed the next chapters in three ways. Drawn on a surface, they define its tangent planes and its first fundamental form (Chapter 19), and their lengths are computed by restricting the ambient metric — the chapter ahead is largely this chapter relativized. The envelope calculus of the weekend problem meets double integrals in Chapter 20, where the astroid’s area is recomputed by Green’s formula (Exercise 20.5) — one curve, two theories, matching answers. And the Frenet system already used the linear differential equations of Chapter 16 (existence, uniqueness, and the orthogonality-preservation argument of Exercise 18.7): the fundamental theorem of curves is a differential equations theorem wearing geometric clothes.

Remark 18.35 (Method: running the local study)

In practice the classification is a four-step routine. One, differentiate at t0t_0 until the first nonzero derivative appears: its index is pp, its value the vector uu. Two, keep differentiating until a derivative not collinear with uu appears: index qq, vector vv. Three, read the parities (p,q)(p, q) in the table. Four, draw: the curve leaves along +u+u if pp is odd (along uu then back along uu if pp is even), on the side of vv dictated by the sign of YY. Two cautions. The frame (u,v)(u, v) is generally not orthonormal — the table describes positions relative to the tangent line, not angles or distances, so do not read curvature off the picture. And intermediate derivatives collinear with uu are allowed between ranks pp and qq (they only shift the expansion of XX); what must not happen is stopping at the first nonzero derivative and guessing q=p+1q = p + 1: for γ(t)=(t2,t4+t5)\gamma(t) = (t^2, t^4 + t^5) the naive guess q=3q = 3 is wrong, as γ(3)(0)\gamma^{(3)}(0) is still collinear with γ(0)\gamma''(0) — this is precisely Exercise 18.5.

18.6 Exercises

Exercise 18.1

Compute the length of one arch of the cycloid γ(t)=(tsint, 1cost)\gamma(t) = (t - \sin t,\ 1 - \cos t), t[0,2π]t \in [0, 2\pi]. (Use 1cost=2sin2(t/2)1 - \cos t = 2\sin^2(t/2).)

Solution

Solution of Exercise 18.1.

γ(t)=(1cost, sint)\gamma'(t) = (1 - \cos t,\ \sin t), so

γ(t)2=(1cost)2+sin2t=22cost=4sin2t2,\norm{\gamma'(t)}^2 = (1 - \cos t)^2 + \sin^2 t = 2 - 2\cos t = 4\sin^2\tfrac t2 ,

and γ(t)=2sint2\norm{\gamma'(t)} = 2\sin\frac t2 (nonnegative on [0,2π][0, 2\pi]). Hence

L=02π2sint2 ⁣dt=[4cost2]02π=8:L = \int_0^{2\pi} 2\sin\tfrac t2\,\dd t = \Bigl[-4\cos\tfrac t2\Bigr]_0^{2\pi} = 8 :

one arch of the cycloid has length 88 (for a wheel of radius 11) — a famous result of Wren, with no π\pi in sight.

Exercise 18.2

Compute the curvature of the ellipse γ(t)=(acost, bsint)\gamma(t) = (a\cos t,\ b\sin t) (a>b>0a > b > 0) and locate the points of maximal and minimal curvature.

Solution

Solution of Exercise 18.2.

With x=acostx = a\cos t, y=bsinty = b\sin t: x=asintx' = -a\sin t, y=bcosty' = b\cos t, x=acostx'' = -a\cos t, y=bsinty'' = -b\sin t, so by Proposition 18.17

κ(t)=xyyx(x2+y2)3/2=absin2t+abcos2t(a2sin2t+b2cos2t)3/2=ab(a2sin2t+b2cos2t)3/2.\kappa(t) = \frac{x'y'' - y'x''}{(x'^2 + y'^2)^{3/2}} = \frac{ab\sin^2 t + ab\cos^2 t} {(a^2\sin^2 t + b^2\cos^2 t)^{3/2}} = \frac{ab}{(a^2\sin^2 t + b^2\cos^2 t)^{3/2}} .

The denominator is minimal when sint=0\sin t = 0 (value b3b^3, points (±a,0)(\pm a, 0)) and maximal when cost=0\cos t = 0 (value a3a^3, points (0,±b)(0, \pm b)), since a>ba > b. Hence κ\kappa is maximal at the ends of the major axis, κmax=a/b2\kappa_{\max} = a/b^2, and minimal at the ends of the minor axis, κmin=b/a2\kappa_{\min} = b/a^2: the ellipse bends most sharply at the tips of its long axis.

Exercise 18.3

Show that the arc length of the graph of f(x)=coshxf(x) = \cosh x over [0,x][0, x] equals sinhx\sinh x, and compute the curvature of this curve (the catenary). Verify that R(x)=1/κ(x)=cosh2xR(x) = 1/\kappa(x) = \cosh^2 x.

Solution

Solution of Exercise 18.3.

For the graph γ(x)=(x,coshx)\gamma(x) = (x, \cosh x): γ(x)=1+sinh2x=coshx\norm{\gamma'(x)} = \sqrt{1 + \sinh^2 x} = \cosh x, so the arc length from 00 to xx is 0xcoshu ⁣du=sinhx\int_0^x \cosh u\,\dd u = \sinh x. Curvature of a graph (Proposition 18.17):

κ(x)=f(x)(1+f(x)2)3/2=coshxcosh3x=1cosh2x,\kappa(x) = \frac{f''(x)}{(1 + f'(x)^2)^{3/2}} = \frac{\cosh x}{\cosh^3 x} = \frac{1}{\cosh^2 x} ,

so R(x)=cosh2xR(x) = \cosh^2 x, as announced. Note the neat coincidence R(x)=1+s(x)2R(x) = 1 + s(x)^2 with s=sinhxs = \sinh x the arc length: the catenary’s radius of curvature grows with the square of the arc length from the vertex.

Exercise 18.4 ★★

(Logarithmic spiral) Let γ(t)=et(cost, sint)\gamma(t) = e^{t}(\cos t,\ \sin t), tRt \in \R. Show that the angle between γ(t)\gamma(t) and γ(t)\gamma'(t) is constant, compute the arc length of γ\gamma on (,0](-\infty, 0] (finite!), and the curvature.

Solution

Solution of Exercise 18.4.

γ(t)=et(costsint, sint+cost)\gamma'(t) = e^t(\cos t - \sin t,\ \sin t + \cos t), so

γ(t),γ(t)=e2t(cost(costsint)+sint(sint+cost))=e2t,\langle \gamma(t), \gamma'(t)\rangle = e^{2t} \bigl(\cos t(\cos t - \sin t) + \sin t(\sin t + \cos t)\bigr) = e^{2t},

while γ(t)=et\norm{\gamma(t)} = e^t and γ(t)=et2\norm{\gamma'(t)} = e^t\sqrt 2. Hence

cos(γ,γ)=e2tetet2=12:\cos\angle\bigl(\gamma, \gamma'\bigr) = \frac{e^{2t}}{e^t \cdot e^t\sqrt2} = \frac{1}{\sqrt2} :

the tangent always makes the angle π/4\pi/4 with the radius — the equiangular property of the logarithmic spiral. Arc length on (,0](-\infty, 0]:

0γ(t) ⁣dt=20et ⁣dt=2,\int_{-\infty}^0 \norm{\gamma'(t)}\,\dd t = \sqrt2\int_{-\infty}^0 e^t\,\dd t = \sqrt 2 ,

finite although the spiral winds infinitely many times around the origin. Curvature: with xyyxx'y'' - y'x'' computed from γ=et(2sint, 2cost)\gamma'' = e^t(-2\sin t,\ 2\cos t),

xyyx=e2t(2cost(costsint)+2sint(sint+cost))=2e2t,x'y'' - y'x'' = e^{2t}\bigl(2\cos t(\cos t - \sin t) + 2\sin t(\sin t + \cos t)\bigr) = 2e^{2t},

so κ(t)=2e2t(et2)3=1et2\kappa(t) = \dfrac{2e^{2t}}{(e^t\sqrt2)^3} = \dfrac{1}{e^t\sqrt2}: the curvature is 1/(2γ)1/(\sqrt2\, \norm{\gamma}), decaying as the spiral grows.

Exercise 18.5 ★★

Determine pp, qq and the local shape (ordinary, inflection, cusp) of γ(t)=(t2, t4+t5)\gamma(t) = (t^2,\ t^4 + t^5) at t=0t = 0, and of γ(t)=(t3, t4)\gamma(t) = (t^3,\ t^4) at t=0t = 0.

Solution

Solution of Exercise 18.5.

First arc: γ(t)=(t2, t4+t5)\gamma(t) = (t^2,\ t^4 + t^5). Derivatives at 00: γ=(2,0)0\gamma'' = (2, 0) \neq 0, so p=2p = 2. Then γ(3)(0)=(0,0)\gamma^{(3)}(0) = (0, 0), γ(4)(0)=(0,24)\gamma^{(4)}(0) = (0, 24), not collinear with (2,0)(2, 0): q=4q = 4. Both even: cusp of the second kind — both branches leave in the direction +u=(1,0)+u = (1,0) and stay on the same side of the tangent. (Indeed y=x2±x5/2y = x^2 \pm x^{5/2} on the two branches: same sign for small xx.)

Second arc: γ(t)=(t3,t4)\gamma(t) = (t^3, t^4). γ(0)=γ(0)=0\gamma'(0) = \gamma''(0) = 0, γ(3)(0)=(6,0)\gamma^{(3)}(0) = (6, 0): p=3p = 3, odd. Next γ(4)(0)=(0,24)\gamma^{(4)}(0) = (0, 24): q=4q = 4, even. Odd–even: ordinary point — despite the vanishing velocity, the trajectory y=x4/3y = x^{4/3} crosses the origin smoothly, staying above its tangent y=0y = 0.

Exercise 18.6 ★★

Let γ\gamma be a unit-speed plane arc with κ(s)>0\kappa(s) > 0 for all ss, and let c(s)=γ(s)+1κ(s)N(s)c(s) = \gamma(s) + \frac{1}{\kappa(s)}N(s) be the center of curvature (the curve cc is the evolute). Assuming κ\kappa is C1\mathcal{C}^1, show that c(s)=κ(s)κ(s)2N(s)c'(s) = -\frac{\kappa'(s)}{\kappa(s)^2}N(s): the evolute is tangent to the normal lines of γ\gamma.

Solution

Solution of Exercise 18.6.

Differentiate c(s)=γ(s)+1κ(s)N(s)c(s) = \gamma(s) + \dfrac{1}{\kappa(s)}N(s) using the plane Frenet formulas (Theorem 18.13):

c(s)=T(s)κ(s)κ(s)2N(s)+1κ(s)(κ(s)T(s))=κ(s)κ(s)2N(s),c'(s) = T(s) - \frac{\kappa'(s)}{\kappa(s)^2}N(s) + \frac{1}{\kappa(s)}\,\bigl(-\kappa(s)T(s)\bigr) = -\frac{\kappa'(s)}{\kappa(s)^2}\,N(s) ,

the tangent terms cancelling exactly. So the velocity of the evolute is carried by N(s)N(s), which directs the normal line of γ\gamma at γ(s)\gamma(s) — and the point c(s)c(s) lies on that very normal line: the evolute is the envelope of the normals. (Where κ=0\kappa' = 0 the evolute has a singular point; this is what produces the cusps of the evolute of an ellipse.)

Exercise 18.7 ★★★

Let A(s)A(s) be a continuous family of antisymmetric 3×33 \times 3 matrices and F=FAF' = F A a matrix solution with F(s0)F(s_0) orthogonal. Show that F(s)F(s) is orthogonal for all ss. (Differentiate G=FFTG = F F^{\mathsf T} and use uniqueness in Cauchy–Lipschitz.) Explain the relevance to the Frenet system.

Solution

Solution of Exercise 18.7.

Let G(s)=F(s)F(s)TG(s) = F(s)F(s)^{\mathsf T}. Then, using F=FAF' = FA and (FT)=(F)T=ATFT(F^{\mathsf T})' = (F')^{\mathsf T} = A^{\mathsf T}F^{\mathsf T},

G=FFT+F(FT)=FAFT+FATFT=F(A+AT)FT=0G' = F'F^{\mathsf T} + F(F^{\mathsf T})' = FAF^{\mathsf T} + FA^{\mathsf T}F^{\mathsf T} = F(A + A^{\mathsf T})F^{\mathsf T} = 0

by antisymmetry. So GG is constant on the interval, equal to G(s0)=F(s0)F(s0)T=IG(s_0) = F(s_0)F(s_0)^{\mathsf T} = I: F(s)F(s) is orthogonal for every ss. (Alternatively, without computing GG' to zero: both GG and the constant II solve the linear system Y=YA+ATYY' = YA + A^{\mathsf T}Y with the same initial value, and Cauchy–Lipschitz uniqueness for linear systems, Chapter 16, forces GIG \equiv I.)

Relevance: the Frenet system (T,N,B)=(T,N,B)A(s)(T, N, B)' = (T, N, B)\,A(s) has the antisymmetric coefficient matrix

A=(0κ0κ0τ0τ0)A = \begin{pmatrix} 0 & -\kappa & 0\\ \kappa & 0 & -\tau\\ 0 & \tau & 0\end{pmatrix}

(columns expressing T,N,BT', N', B'). The computation above shows that a solution frame that starts orthonormal stays orthonormal — the key step in the fundamental theorem reconstructing a curve from (κ,τ)(\kappa, \tau).

Exercise 18.8 ★★★

(Total curvature of a closed convex curve) Let γ~\tilde\gamma be a unit-speed C2\mathcal{C}^2 closed plane arc of length LL (so γ~(s+L)=γ~(s)\tilde\gamma(s + L) = \tilde\gamma(s)), traversed once counterclockwise. Using the angle function φ\varphi with T=(cosφ,sinφ)T = (\cos\varphi, \sin\varphi) from Theorem 18.20, explain why φ(L)φ(0)\varphi(L) - \varphi(0) is a multiple of 2π2\pi, and show that 0Lκ(s) ⁣ds=φ(L)φ(0)\int_0^L \kappa(s)\,\dd s = \varphi(L) - \varphi(0). (For a circle of radius RR: κ=1R2πR=2π\int \kappa = \frac1R \cdot 2\pi R = 2\pi. The theorem of turning tangents asserts the value 2π2\pi for every simple closed curve; you are not asked to prove that.)

Solution

Solution of Exercise 18.8.

By Theorem 18.20 (uniqueness part), there is a C1\mathcal{C}^1 angle function φ\varphi with T(s)=(cosφ(s),sinφ(s))T(s) = (\cos\varphi(s), \sin\varphi(s)) and φ=κ\varphi' = \kappa. Hence

0Lκ(s) ⁣ds=φ(L)φ(0).\int_0^L \kappa(s)\,\dd s = \varphi(L) - \varphi(0) .

Since the arc is closed of period LL, T(L)=T(0)T(L) = T(0): (cosφ(L),sinφ(L))=(cosφ(0),sinφ(0))(\cos\varphi(L), \sin\varphi(L)) = (\cos\varphi(0), \sin\varphi(0)), so φ(L)φ(0)2πZ\varphi(L) - \varphi(0) \in 2\pi\Z. The total curvature of a closed curve is therefore always an integer multiple of 2π2\pi — the integer being the winding number of the tangent (the number of full turns TT makes). For the circle of radius RR: κ=1/R\kappa = 1/R and L=2πRL = 2\pi R, total curvature 2π2\pi, winding number 11; the theorem of turning tangents states this value holds for every simple closed curve.

Exercise 18.9 ★★★

Show that a biregular space curve with constant κ>0\kappa > 0 and τ=0\tau = 0 is (an arc of) a circle of radius 1/κ1/\kappa. (Use Proposition 18.27, then show the center γ+1κN\gamma + \frac1\kappa N is constant.)

Solution

Solution of Exercise 18.9.

Since τ0\tau \equiv 0, the curve lies in a plane (Proposition 18.27); work in that plane. Consider the candidate center

c(s)=γ(s)+1κN(s)(κ constant).c(s) = \gamma(s) + \frac{1}{\kappa}N(s) \qquad (\kappa \text{ constant}).

Differentiating with the Frenet formulas (N=κT+τB=κTN' = -\kappa T + \tau B = -\kappa T here):

c(s)=T+1κ(κT)=0,c'(s) = T + \frac1\kappa(-\kappa T) = 0 ,

so cc is a constant point Ω\Omega. Then γ(s)Ω=1κN(s)=1κ\norm{\gamma(s) - \Omega} = \norm{-\frac1\kappa N(s)} = \frac1\kappa for all ss: the curve lies on the circle of center Ω\Omega and radius 1/κ1/\kappa (in its plane), and being a nonconstant arc of it, it is an arc of that circle.

Exercise 18.10

Compute the arc length of the parabola y=x2/2y = x^2/2 over [0,a]\intcc0a and show that it equals

12(a1+a2+ln(a+1+a2)).\tfrac12\Bigl(a\sqrt{1 + a^2} + \ln\bigl(a + \sqrt{1 + a^2}\bigr)\Bigr).
Solution

Solution of Exercise 18.10.

For the graph γ(x)=(x,x2/2)\gamma(x) = (x, x^2/2), γ(x)=1+x2\norm{\gamma'(x)} = \sqrt{1 + x^2}, so L=0a1+x2 ⁣dxL = \int_0^a\sqrt{1+x^2}\,\dd x. Substituting x=sinhux = \sinh u ( ⁣dx=coshu ⁣du\dd x = \cosh u\,\dd u, uu from 00 to ua=ln(a+1+a2)u_a = \ln(a + \sqrt{1+a^2})):

L=0uacosh2u ⁣du=12[u+sinhucoshu]0ua=12(ln(a+1+a2)+a1+a2),L = \int_0^{u_a}\cosh^2 u\,\dd u = \frac12\bigl[u + \sinh u\cosh u\bigr]_0^{u_a} = \frac12\Bigl(\ln\bigl(a + \sqrt{1+a^2}\bigr) + a\sqrt{1+a^2}\Bigr),

using cosh2u=1+cosh2u2\cosh^2 u = \frac{1 + \cosh 2u}2 and sinhua=a\sinh u_a = a, coshua=1+a2\cosh u_a = \sqrt{1 + a^2}.

Exercise 18.11 ★★

Let γ\gamma be a regular C2\mathcal C^2 arc in Rn\R^n all of whose tangent lines pass through a fixed point PP. Prove that the trajectory of γ\gamma is contained in a straight line. (Parametrize by arc length, write γ(s)+λ(s)T(s)=P\gamma(s) + \lambda(s)T(s) = P and differentiate.)

Solution

Solution of Exercise 18.11.

Parametrize by arc length (Theorem 18.9) and set λ(s)=Pγ(s),T(s)\lambda(s) = \langle P - \gamma(s), T(s)\rangle, a C1\mathcal C^1 function; since PP lies on the tangent line at γ(s)\gamma(s), the vector Pγ(s)P - \gamma(s) is collinear with T(s)T(s), so P=γ(s)+λ(s)T(s)P = \gamma(s) + \lambda(s)T(s). Differentiating,

0=T(s)+λ(s)T(s)+λ(s)T(s)=(1+λ(s))T(s)+λ(s)T(s),0 = T(s) + \lambda'(s)T(s) + \lambda(s)T'(s) = \bigl(1 + \lambda'(s)\bigr)T(s) + \lambda(s)T'(s),

and T(s)T(s)T'(s) \perp T(s) (differentiate T2=1\norm T^2 = 1), so both components vanish: λ=1\lambda' = -1 and λT=0\lambda T' = 0. Then λ(s)=cs\lambda(s) = c - s vanishes at most once, so T=0T' = 0 on a dense set, hence everywhere by continuity: TT is a constant unit vector and γ(s)=γ(s0)+(ss0)T\gamma(s) = \gamma(s_0) + (s - s_0)T: a straight line (through PP, as it must be).

Exercise 18.12 ★★★

(Fundamental theorem for space curves) Let κ>0\kappa > 0 and τ\tau be continuous functions on an interval JJ. Carry out the program of the remark following Example 18.29: (a) show that the linear system F=FA(s)F' = FA(s), with A(s)A(s) the antisymmetric Frenet matrix built from κ,τ\kappa, \tau and F(s0)F(s_0) a direct orthonormal frame, has a unique global solution, which remains a direct orthonormal frame; (b) construct a unit-speed biregular curve with curvature κ\kappa and torsion τ\tau; (c) prove uniqueness up to a direct isometry of R3\R^3.

Solution

Solution of Exercise 18.12.

(a) The Frenet matrix

A(s)=(0κ0κ0τ0τ0)A(s) = \begin{pmatrix} 0 & -\kappa & 0\\ \kappa & 0 & -\tau\\ 0 & \tau & 0\end{pmatrix}

has continuous entries, so the linear system F=FA(s)F' = FA(s), F(s0)=F0F(s_0) = F_0 (a direct orthonormal matrix) has a unique solution on all of JJ (Theorem 16.4). By Exercise 18.7, F(s)F(s) is orthogonal for every ss; detF\det F is continuous with values in {±1}\{\pm1\} and equals 11 at s0s_0, so F(s)F(s) is direct for all ss.

(b) Read off the rows T,N,BT, N, B of FF (so that T=κNT' = \kappa N, N=κT+τBN' = -\kappa T + \tau B, B=τNB' = -\tau N) and set γ(s)=γ0+s0sT(u) ⁣du\gamma(s) = \gamma_0 + \int_{s_0}^s T(u)\,\dd u. Then γ=T\gamma' = T is a unit vector: unit speed; T=κNT' = \kappa N with κ>0\kappa > 0 and NN unit orthogonal to TT, so γ\gamma is biregular with curvature T=κ\norm{T'} = \kappa and principal normal NN; the binormal is TN=BT \wedge N = B (direct orthonormal frame), and B=τNB' = -\tau N identifies the torsion as τ\tau.

(c) Let γ1,γ2\gamma_1, \gamma_2 be unit-speed biregular curves with the same (κ,τ)(\kappa, \tau). There is a unique direct isometry Φ=ρ+w\Phi = \rho + w (ρSO(3)\rho \in SO(3)) sending γ1(s0)\gamma_1(s_0) to γ2(s0)\gamma_2(s_0) and the Frenet frame of γ1\gamma_1 at s0s_0 to that of γ2\gamma_2 at s0s_0. The curve Φγ1\Phi\circ\gamma_1 is unit-speed with the same invariants (its frame is ρ\rho applied to that of γ1\gamma_1, and ρ\rho preserves cross products, being direct). Now the frames of Φγ1\Phi\circ\gamma_1 and γ2\gamma_2 both solve F=FA(s)F' = FA(s) with the same initial value, so they coincide by uniqueness; in particular the tangents agree, and integrating from the common point s0s_0: Φγ1=γ2\Phi\circ\gamma_1 = \gamma_2.

18.7 Problem: envelopes — the astroid, two evolutes, and a caustic

A ladder of length 1 sliding down a wall (blue positions) never crosses the astroid x2/3 + y2/3 = 1 (red): the astroid is the envelope of the family of segments, tangent to every one of them.
A ladder of length 11 sliding down a wall (blue positions) never crosses the astroid x2/3+y2/3=1x^{2/3} + y^{2/3} = 1 (red): the astroid is the envelope of the family of segments, tangent to every one of them.

Problem 18.1

Weekend problem — the envelope machine and four classical curves

A one-parameter family of lines usually fails to cover the plane evenly: the lines pile up along a curve tangent to all of them, their envelope. Light rays make envelopes visible as caustics — the bright cusped curve in a mug of coffee. This problem builds the general envelope machine, then runs it four times: the sliding ladder (astroid), the normals of the parabola and of the cycloid (evolutes, with Huygens’ pendulum at the end), and the coffee-cup caustic (nephroid). Throughout, DtD_t denotes the line of equation a(t)x+b(t)y=c(t)a(t)\,x + b(t)\,y = c(t), where a,b,ca, b, c are C2\mathcal C^2 functions with (a(t),b(t))(0,0)(a(t), b(t)) \neq (0,0), and Δ(t)=a(t)b(t)a(t)b(t)\Delta(t) = a(t)b'(t) - a'(t)b(t).

Part I — The envelope machine.

  1. Suppose Δ(t)0\Delta(t) \neq 0. Show that the characteristic system

    {a(t)x+b(t)y=c(t)a(t)x+b(t)y=c(t)\begin{cases} a(t)\,x + b(t)\,y = c(t)\\ a'(t)\,x + b'(t)\,y = c'(t)\end{cases}

    has a unique solution E(t)=(x(t),y(t))E(t) = (x(t), y(t)), given by x=cbcbΔx = \dfrac{cb' - c'b}{\Delta}, y=acacΔy = \dfrac{ac' - a'c}{\Delta}.

  2. Assume moreover that EE is C1\mathcal C^1 near tt with E(t)0E'(t) \neq 0. Differentiating the first equation of the system, show a(t)x(t)+b(t)y(t)=0a(t)\,x'(t) + b(t)\,y'(t) = 0, and conclude that the curve EE passes through a point of DtD_t with the direction of DtD_t: the family is tangent to EE, which is called its envelope.
  3. Sanity check: the tangent lines of the parabola y=x2/2y = x^2/2 at the points (t,t2/2)(t, t^2/2) are txy=t2/2tx - y = t^2/2. Verify that the envelope machine returns the parabola itself.
  4. (Envelope of the normals) Let γ\gamma be unit-speed with κ(s)0\kappa(s) \neq 0. The normal line at γ(s)\gamma(s) is {M:Mγ(s),T(s)=0}\{M : \langle M - \gamma(s), T(s) \rangle = 0\}. Show that its characteristic system forces Mγ(s),N(s)=1/κ(s)\langle M - \gamma(s), N(s)\rangle = 1/\kappa(s), hence that the characteristic point is the center of curvature: the envelope of the normals is the evolute, recovering Exercise 18.6. Check Δ(s)=κ(s)\Delta(s) = \kappa(s).
  5. Two degenerations. For the pencil Dθ:xcosθ+ysinθ=0D_\theta : x\cos \theta + y\sin\theta = 0, show that the characteristic point is the origin for every θ\theta (the “envelope” collapses to a point, and E=0E' = 0: question 2 does not apply). For a family of parallel lines (a,ba, b constant), show Δ0\Delta \equiv 0 and that the characteristic system is in general inconsistent: no envelope.

Part II — The sliding ladder and the astroid. A segment of length 11 slides with one end Pt=(cost,0)P_t = (\cos t, 0) on the floor and the other Qt=(0,sint)Q_t = (0, \sin t) on the wall, t(0,π/2)t \in \intoo0{\pi/2}.

  1. Show that the line (PtQt)(P_tQ_t) has equation xsint+ycost=sintcostx\sin t + y\cos t = \sin t\cos t, and that the envelope machine gives the characteristic point

    E(t)=(cos3t, sin3t):E(t) = (\cos^3 t,\ \sin^3 t) :

    the astroid, of implicit equation x2/3+y2/3=1x^{2/3} + y^{2/3} = 1 (extended to the other quadrants by symmetry).

  2. Show that E(0)=0E'(0) = 0 and, using the local classification (Proposition 18.32), that the astroid has a cusp of the first kind at (1,0)(1, 0) — and likewise at its four axis points.
  3. Compute E(t)=32sin2t\norm{E'(t)} = \tfrac32\abs{\sin 2t} and deduce that the total length of the astroid is 66.
  4. Where does the ladder touch the astroid? Show E(t)=Pt+sin2t(QtPt)E(t) = P_t + \sin^2 t\,(Q_t - P_t): the contact point divides the ladder in the ratio sin2t:cos2t\sin^2 t : \cos^2 t, sweeping it from one end to the other as the ladder slides.
  5. Compute the area enclosed by the astroid: show that the first-quadrant area is 30π/2sin4tcos2t ⁣dt3\int_0^{\pi/2}\sin^4 t\cos^2 t\,\dd t, evaluate the integral by linearization (sin22t=1cos4t2\sin^2 2t = \tfrac{1 - \cos 4t}2), and conclude that the total area is 3π/83\pi/8.

Part III — The evolute of the parabola. Let γ(t)=(t,t2/2)\gamma(t) = (t, t^2/2).

  1. Show that the normal line at γ(t)\gamma(t) has equation x+ty=t+t3/2x + t\,y = t + t^3/2.
  2. Run the envelope machine: show that the envelope of the normals is

    E(t)=(t3, 1+32t2),E(t) = \Bigl(-t^3,\ 1 + \tfrac32 t^2\Bigr),

    with implicit equation x2=827(y1)3x^2 = \tfrac8{27}(y - 1)^3: a semicubical parabola.

  3. Cross-check with question 4: compute the center of curvature γ(t)+1κ(t)N(t)\gamma(t) + \frac1{\kappa(t)}N(t) from κ(t)=(1+t2)3/2\kappa(t) = (1 + t^2)^{-3/2} (Example 18.18) and recover the same point.
  4. Show that the evolute has a cusp of the first kind at (0,1)(0, 1), the center of curvature at the vertex — the point where κ\kappa is extremal, as predicted by the formula c=κκ2Nc' = -\frac{\kappa'}{\kappa^2}N of Exercise 18.6.
  5. How many normals of the parabola pass through a given point (x0,y0)(x_0, y_0)? Show the answer is governed by the cubic t32+(1y0)tx0=0\tfrac{t^3}2 + (1 - y_0)\,t - x_0 = 0; treat the axis case x0=0x_0 = 0 completely (one normal for y0<1y_0 < 1, three for y0>1y_0 > 1), and interpret the evolute as the transition curve.

Part IV — The coffee-cup caustic. Parallel rays of direction (1,0)(1, 0) strike the inside of the mirror circle x2+y2=1x^2 + y^2 = 1; the ray hitting Pθ=(cosθ,sinθ)P_\theta = (\cos\theta, \sin\theta) reflects according to the law of reflection.

  1. From the mirror symmetry in the normal (the radius), justify that the reflected direction is v=u2u,nnv = u - 2\langle u, n\rangle n with u=(1,0)u = (1,0), n=(cosθ,sinθ)n = (\cos\theta, \sin\theta), and compute v=(cos2θ,sin2θ)v = -(\cos2\theta, \sin2\theta).
  2. Show that the reflected ray lies on the line

    xsin2θycos2θ=sinθ.x\sin 2\theta - y\cos 2\theta = \sin\theta .
  3. Run the envelope machine (Δ=2\Delta = 2): show that the caustic is

    E(θ)=(3cosθcos3θ4, 3sinθsin3θ4),E(\theta) = \Bigl(\tfrac{3\cos\theta - \cos3\theta}4,\ \tfrac{3\sin\theta - \sin3\theta}4\Bigr),

    the nephroid.

  4. Compute E(θ)=32sinθ(cos2θ,sin2θ)E'(\theta) = \tfrac32\sin\theta\, (\cos2\theta, \sin2\theta); check that the tangent direction is the reflected-ray direction (question 16), locate the two cusps (±12,0)(\pm\tfrac12, 0), and show that the reflected ray crosses the axis y=0y = 0 at x=12cosθx = \frac1{2\cos\theta} — so nearly axial rays focus at x=12x = \tfrac12: the focal length R/2R/2 of a mirror of radius RR.
  5. Show that the nephroid has total length 66 and that near θ=0\theta = 0,

    E(θ)(12,0)=(34θ2+o(θ2), θ3+o(θ3)):E(\theta) - \bigl(\tfrac12, 0\bigr) = \bigl(\tfrac34\theta^2 + o(\theta^2),\ \theta^3 + o(\theta^3)\bigr) :

    a cusp of the first kind, pointing along the axis.

  6. Explain in one paragraph why the caustic is bright: through every point just outside the caustic pass two reflected rays, through every point on it the rays are “infinitely concentrated” (the map (θ,distance along ray)R2(\theta, \text{distance along ray}) \mapsto \R^2 has a critical point exactly on the envelope).

Part V — Huygens: the cycloid is its own evolute. Let γ(t)=(tsint, 1cost)\gamma(t) = (t - \sin t,\ 1 - \cos t), t(0,2π)t \in \intoo0{2\pi}, one arch of the cycloid.

  1. Compute κ(t)=14sin(t/2)\kappa(t) = -\dfrac1{4\sin(t/2)} and the center of curvature; show that the evolute is

    c(t)=(t+sint, cost1),c(t) = (t + \sin t,\ \cos t - 1),

    and that the substitution t=u+πt = u + \pi exhibits it as the original cycloid translated by (π,2)(\pi, -2): the evolute of a cycloid is a congruent cycloid (Huygens).

  2. Verify that the radius of curvature at the apex t=πt = \pi equals 44, which is half the length 88 of one arch (Exercise 18.1); locate the cusp of the evolute directly below the apex, at distance 44.
  3. (The string property) Let γ\gamma be unit-speed with κ>0\kappa > 0, κ\kappa of class C1\mathcal C^1 and R=1/κR = 1/\kappa strictly monotone. Using c=RNc' = R'N, show that the arc length of the evolute between c(s0)c(s_0) and c(s1)c(s_1) is R(s1)R(s0)\abs{R(s_1) - R(s_0)}. Interpret: a taut string unwound from the evolute, of length R(s0)R(s_0) at the start, has its free end trace the original curve — so a pendulum swinging between two cycloidal cheeks of Huygens’ clock describes a cycloid.
  4. Synthesis. The machine of Part I produced the astroid, a semicubical parabola, a nephroid and a cycloid. For each of the four families, state in one sentence where the hypotheses Δ0\Delta \neq 0 and E0E' \neq 0 held or failed, and what geometric event (cusp, focus, degeneracy) each failure of E0E' \neq 0 signalled. Where must extrema of curvature appear on the envelope of the normals, and why?
Solution

Solution of Problem 18.1.

1. The system is linear in (x,y)(x, y) with determinant Δ(t)=abab0\Delta(t) = a b' - a'b \neq 0: Cramer’s rule gives the unique solution

x=cbcbΔ,y=acacΔ.x = \frac{c b' - c' b}{\Delta}, \qquad y = \frac{a c' - a' c}{\Delta}.

2. Since a(t)x(t)+b(t)y(t)=c(t)a(t)x(t) + b(t)y(t) = c(t) identically, differentiating gives ax+by+ax+by=ca'x + b'y + ax' + by' = c'; the second characteristic equation kills ax+byca'x + b'y - c', so a(t)x(t)+b(t)y(t)=0a(t)x'(t) + b(t)y'(t) = 0: E(t)E'(t) is orthogonal to (a,b)(a, b), hence parallel to (b,a)(-b, a), the direction of DtD_t. As E(t)DtE(t) \in D_t (first equation) and E(t)0E'(t) \neq 0, the line DtD_t is exactly the tangent line of the curve EE at E(t)E(t).

3. Here (a,b,c)=(t,1,t2/2)(a, b, c) = (t, -1, t^2/2), so Δ=t01(1)=1\Delta = t\cdot0 - 1\cdot(-1) = 1 and

x=cbcbΔ=t220+t=t,y=acacΔ=ttt22=t22:x = \frac{c b' - c' b}{\Delta} = \tfrac{t^2}2\cdot 0 + t = t, \qquad y = \frac{a c' - a' c}{\Delta} = t\cdot t - \tfrac{t^2}2 = \tfrac{t^2}2 :

the envelope of the tangent lines of the parabola is the parabola, as it should be.

4. The normal line is M,T(s)=γ(s),T(s)\langle M, T(s)\rangle = \langle\gamma(s), T(s)\rangle: coefficients a=T1a = T_1, b=T2b = T_2, c=γ,Tc = \langle\gamma, T\rangle. Differentiating with Frenet (T=κNT' = \kappa N): a=κN1a' = \kappa N_1, b=κN2b' = \kappa N_2, and c=T,T+γ,κN=1+κγ,Nc' = \langle T, T\rangle + \langle\gamma, \kappa N\rangle = 1 + \kappa\langle\gamma, N\rangle. The second characteristic equation κM,N=1+κγ,N\kappa\langle M, N\rangle = 1 + \kappa\langle\gamma, N\rangle reads κMγ,N=1\kappa\langle M - \gamma, N\rangle = 1. The first says MγTM - \gamma \perp T, so Mγ=μNM - \gamma = \mu N with μ=1/κ\mu = 1/\kappa: the characteristic point is γ+1κN\gamma + \frac1\kappa N, the center of curvature, and the envelope of the normals is the evolute of Exercise 18.6. Finally Δ=T1κN2κN1T2=κdet(T,N)=κ0\Delta = T_1\kappa N_2 - \kappa N_1 T_2 = \kappa\det(T, N) = \kappa \neq 0.

5. Pencil: the system xcosθ+ysinθ=0x\cos\theta + y\sin\theta = 0, xsinθ+ycosθ=0-x\sin\theta + y\cos\theta = 0 has determinant 11 and solution (0,0)(0,0) for every θ\theta: E(0,0)E \equiv (0,0), E0E' \equiv 0, and there is no curve — merely the common point of all the lines. Parallel family: a=b=0a' = b' = 0 gives Δ0\Delta \equiv 0 and the second equation 0=c(t)0 = c'(t), which fails as soon as the family actually moves: no characteristic point, and indeed a family of parallel lines touches no curve along all its members.

6. The line through (cost,0)(\cos t, 0) and (0,sint)(0, \sin t) is xcost+ysint=1\frac x{\cos t} + \frac y{\sin t} = 1, i.e. xsint+ycost=sintcostx\sin t + y\cos t = \sin t\cos t. With (a,b,c)=(sint,cost,sintcost)(a, b, c) = (\sin t, \cos t, \sin t\cos t): a=costa' = \cos t, b=sintb' = -\sin t, c=cos2tc' = \cos 2t, Δ=sin2tcos2t=1\Delta = -\sin^2 t - \cos^2 t = -1. Cramer:

x=cbcb1=sin2tcost+cos2tcost=cost(sin2t+cos2tsin2t)=cos3t,y=acac1=sintcos2tsintcos2t=sint(cos2tcos2t+sin2t)=sin3t.\begin{align*} x &= \frac{cb' - c'b}{-1} = \sin^2 t\cos t + \cos 2t\cos t = \cos t\,(\sin^2 t + \cos^2 t - \sin^2 t) = \cos^3 t,\\ y &= \frac{ac' - a'c}{-1} = \sin t\cos^2 t - \sin t\cos 2t = \sin t\,(\cos^2 t - \cos^2 t + \sin^2 t) = \sin^3 t . \end{align*}

And (cos3t)2/3+(sin3t)2/3=1(\cos^3t)^{2/3} + (\sin^3t)^{2/3} = 1: the astroid.

7. E(t)=3(cos2tsint, sin2tcost)E'(t) = 3(-\cos^2 t\sin t,\ \sin^2 t\cos t) vanishes at t=0t = 0. There, E(0)=(3,0)0E''(0) = (-3, 0) \neq 0 gives p=2p = 2; the xx-component of EE is even in tt, so E(0)=(0,6)E'''(0) = (0, 6), not collinear: q=3q = 3. Even–odd: cusp of the first kind at (1,0)(1, 0) (Proposition 18.32), with tangent along the xx-axis. The symmetries xxx \mapsto -x, yyy \mapsto -y, (x,y)(y,x)(x, y) \mapsto (y, x) of the astroid transport the cusp to (1,0)(-1, 0) and (0,±1)(0, \pm1).

8. E(t)=3sintcost(cost,sint)E'(t) = 3\sin t\cos t\,(-\cos t, \sin t), so E(t)=3sintcost=32sin2t\norm{E'(t)} = 3\abs{\sin t\cos t} = \tfrac32\abs{\sin 2t}. One quadrant: 0π/232sin2t ⁣dt=32\int_0^{\pi/2}\tfrac32\sin 2t\,\dd t = \tfrac32, and by symmetry the total length is 432=64 \cdot \tfrac32 = 6.

9. E(t)Pt=(cos3tcost, sin3t)=sin2t(cost, sint)=sin2t(QtPt)E(t) - P_t = (\cos^3 t - \cos t,\ \sin^3 t) = \sin^2 t\,(-\cos t,\ \sin t) = \sin^2 t\,(Q_t - P_t). So the contact point is the barycenter of (Pt,cos2t)(P_t, \cos^2 t) and (Qt,sin2t)(Q_t, \sin^2 t): as tt runs from 00 to π/2\pi/2 it slides from the floor end to the wall end of the ladder.

10. In the first quadrant the region under the astroid has area 01y ⁣dx\int_0^1 y\,\dd x with x=cos3tx = \cos^3 t decreasing from 11 to 00 as tt goes from 00 to π/2\pi/2:

01y ⁣dx=π/20sin3t(3cos2tsint) ⁣dt=30π/2sin4tcos2t ⁣dt.\int_0^1 y\,\dd x = \int_{\pi/2}^{0}\sin^3 t\,(-3\cos^2 t\sin t)\,\dd t = 3\int_0^{\pi/2}\sin^4 t\cos^2 t\,\dd t .

Linearize: sin4tcos2t=(sintcost)2sin2t=18(sin22tsin22tcos2t)\sin^4 t\cos^2 t = (\sin t\cos t)^2\sin^2 t = \tfrac18\bigl(\sin^2 2t - \sin^2 2t\cos 2t\bigr), and 0π/2sin22t ⁣dt=π4\int_0^{\pi/2}\sin^2 2t\,\dd t = \tfrac\pi4 while 0π/2sin22tcos2t ⁣dt=[sin32t6]0π/2=0\int_0^{\pi/2}\sin^2 2t\cos 2t\,\dd t = \bigl[\tfrac{\sin^3 2t}6\bigr]_0^{\pi/2} = 0. So the integral is π32\tfrac\pi{32}, the quadrant area 3π32\tfrac{3\pi}{32}, and the enclosed area 43π32=3π84 \cdot \tfrac{3\pi}{32} = \tfrac{3\pi}8.

11. The tangent at γ(t)=(t,t2/2)\gamma(t) = (t, t^2/2) is directed by (1,t)(1, t), so the normal line is {(x,y):(xt)+t(yt2/2)=0}\{(x, y) : (x - t) + t\,(y - t^2/2) = 0\}, i.e. x+ty=t+t32x + t\,y = t + \tfrac{t^3}2.

12. (a,b,c)=(1,t,t+t3/2)(a, b, c) = (1, t, t + t^3/2): a=0a' = 0, b=1b' = 1, c=1+32t2c' = 1 + \tfrac32 t^2, Δ=1\Delta = 1. Then y=acac=1+32t2y = ac' - a'c = 1 + \tfrac32t^2 and

x=cbcb=t+t32t(1+32t2)=t3.x = cb' - c'b = t + \tfrac{t^3}2 - t\Bigl(1 + \tfrac32t^2\Bigr) = -t^3 .

Eliminating tt: t2=23(y1)t^2 = \tfrac23(y - 1) and x2=t6=827(y1)3x^2 = t^6 = \tfrac8{27}(y - 1)^3: a semicubical parabola with vertex (0,1)(0, 1).

13. T=(1,t)/1+t2T = (1, t)/\sqrt{1+t^2}, N=(t,1)/1+t2N = (-t, 1)/\sqrt{1+t^2}, and κ=(1+t2)3/2\kappa = (1+t^2)^{-3/2}, so

γ+1κN=(t,t22)+(1+t2)(t,1)=(t3, 1+32t2),\gamma + \frac1\kappa N = (t, \tfrac{t^2}2) + (1+t^2)\,(-t, 1) = \Bigl(-t^3,\ 1 + \tfrac32t^2\Bigr),

the same curve: the envelope of the normals is the locus of the centers of curvature, as question 4 promised.

14. E(t)=(3t2,3t)E'(t) = (-3t^2, 3t) vanishes at t=0t = 0; E(0)=(0,3)0E''(0) = (0, 3) \neq 0 gives p=2p = 2 and E(0)=(6,0)E'''(0) = (-6, 0) gives q=3q = 3: a cusp of the first kind at (0,1)(0, 1). The vertex is where κ=(1+t2)3/2\kappa = (1+t^2)^{-3/2} is maximal, so κ(0)=0\kappa'(0) = 0 and the evolute’s velocity κκ2N-\frac{\kappa'}{\kappa^2}N vanishes exactly there: cusps of the evolute sit at the extrema of curvature.

15. The normal at parameter tt passes through (x0,y0)(x_0, y_0) iff x0+ty0=t+t32x_0 + t\,y_0 = t + \tfrac{t^3}2, i.e.

t32+(1y0)tx0=0,\frac{t^3}2 + (1 - y_0)\,t - x_0 = 0 ,

a cubic in tt: one or three real roots (counted without multiplicity, for generic points). On the axis x0=0x_0 = 0 it factors as t(t22+1y0)=0t\bigl(\tfrac{t^2}2 + 1 - y_0\bigr) = 0: the root t=0t = 0 (the axis is the normal at the vertex), plus t=±2(y01)t = \pm\sqrt{2(y_0 - 1)} when y0>1y_0 > 1. So: one normal for y0<1y_0 < 1, three for y0>1y_0 > 1, and at y0=1y_0 = 1 the triple root marks the cusp of the evolute. In general a double root of the cubic means the point satisfies both the line equation and its tt-derivative — it lies on the envelope: the evolute is precisely the boundary between the one-normal and three-normal regions.

16. Reflection in the mirror reverses the normal component of the direction and keeps the tangential one: writing u=u,nn+utanu = \langle u, n\rangle n + u_{\mathrm{tan}}, the reflected direction is utanu,nn=u2u,nnu_{\mathrm{tan}} - \langle u, n\rangle n = u - 2\langle u, n\rangle n. Here u,n=cosθ\langle u, n\rangle = \cos\theta, so

v=(1,0)2cosθ(cosθ,sinθ)=(12cos2θ, 2sinθcosθ)=(cos2θ, sin2θ).v = (1, 0) - 2\cos\theta\,(\cos\theta, \sin\theta) = (1 - 2\cos^2\theta,\ -2\sin\theta\cos\theta) = -(\cos2\theta,\ \sin2\theta).

17. The reflected ray passes through Pθ=(cosθ,sinθ)P_\theta = (\cos\theta, \sin\theta) with direction (cos2θ,sin2θ)(\cos2\theta, \sin2\theta); a normal vector is (sin2θ,cos2θ)(-\sin2\theta, \cos2\theta), so the line is

sin2θ(xcosθ)+cos2θ(ysinθ)=0,-\sin2\theta\,(x - \cos\theta) + \cos2\theta\,(y - \sin\theta) = 0,

and the constant is sin2θcosθ+cos2θsinθ=sinθ-\sin2\theta\cos\theta + \cos2\theta\sin\theta = -\sin\theta: multiplying by 1-1, xsin2θycos2θ=sinθx\sin2\theta - y\cos2\theta = \sin\theta.

18. (a,b,c)=(sin2θ,cos2θ,sinθ)(a, b, c) = (\sin2\theta, -\cos2\theta, \sin\theta): a=2cos2θa' = 2\cos2\theta, b=2sin2θb' = 2\sin2\theta, c=cosθc' = \cos\theta, Δ=2sin22θ+2cos22θ=2\Delta = 2\sin^22\theta + 2\cos^22\theta = 2. Cramer, then product-to-sum formulas:

x=2sinθsin2θ+cosθcos2θ2=(cosθcos3θ)+12(cosθ+cos3θ)2=3cosθcos3θ4,y=sin2θcosθ2cos2θsinθ2=12(sin3θ+sinθ)(sin3θsinθ)2=3sinθsin3θ4:\begin{align*} x &= \frac{2\sin\theta\sin2\theta + \cos\theta\cos2\theta}{2} = \frac{(\cos\theta - \cos3\theta) + \frac12(\cos\theta + \cos3\theta)}{2} = \frac{3\cos\theta - \cos3\theta}4,\\ y &= \frac{\sin2\theta\cos\theta - 2\cos2\theta\sin\theta}2 = \frac{\frac12(\sin3\theta + \sin\theta) - (\sin3\theta - \sin\theta)}2 = \frac{3\sin\theta - \sin3\theta}4 : \end{align*}

the nephroid, a closed curve with two cusps.

19. Differentiating and factoring with sin3θsinθ=2cos2θsinθ\sin3\theta - \sin\theta = 2\cos2\theta\sin\theta, cosθcos3θ=2sin2θsinθ\cos\theta - \cos3\theta = 2\sin2\theta\sin\theta:

E(θ)=34(sin3θsinθ, cosθcos3θ)=32sinθ(cos2θ,sin2θ),E'(\theta) = \tfrac34\bigl(\sin3\theta - \sin\theta,\ \cos\theta - \cos3\theta\bigr) = \tfrac32\sin\theta\,(\cos2\theta, \sin2\theta),

parallel to the reflected direction of question 16: each reflected ray is tangent to the caustic, as the envelope property demands. E=0E' = 0 exactly at sinθ=0\sin\theta = 0: E(0)=(12,0)E(0) = (\tfrac12, 0) and E(π)=(12,0)E(\pi) = (-\tfrac12, 0), the two cusps. Setting y=0y = 0 in the line equation: xsin2θ=sinθx\sin2\theta = \sin\theta, so x=12cosθ12x = \frac1{2\cos\theta} \to \frac12 as θ0\theta \to 0: paraxial rays focus at distance R/2R/2 from the center — the focal length of the spherical mirror.

20. E(θ)=32sinθ\norm{E'(\theta)} = \tfrac32\abs{\sin\theta}, so the length is 3202πsinθ ⁣dθ=324=6\tfrac32\int_0^{2\pi}\abs{\sin\theta}\, \dd\theta = \tfrac32\cdot4 = 6. Near θ=0\theta = 0, with coskθ=1k2θ22+O(θ4)\cos k\theta = 1 - \tfrac{k^2\theta^2}2 + O(\theta^4) and sinkθ=kθk3θ36+O(θ5)\sin k\theta = k\theta - \tfrac{k^3\theta^3}6 + O(\theta^5):

x12=3θ2+O(θ4)4=34θ2+O(θ4),y=4θ3+O(θ5)4=θ3+O(θ5):x - \tfrac12 = \frac{3\theta^2 + O(\theta^4)}{4} = \tfrac34\theta^2 + O(\theta^4), \qquad y = \frac{4\theta^3 + O(\theta^5)}{4} = \theta^3 + O(\theta^5) :

p=2p = 2, q=3q = 3, a cusp of the first kind pointing along the axis — the bright point of the coffee-cup caustic.

21. Parametrize the illuminated points by Φ(θ,r)=Pθ+rvθ\Phi(\theta, r) = P_\theta + r\,v_\theta (position along each reflected ray). The Jacobian determinant det(θΦ,rΦ)=det(Pθ+rvθ,vθ)\det(\partial_\theta\Phi, \partial_r\Phi) = \det(P_\theta' + rv_\theta', v_\theta) is affine in rr and vanishes for exactly one r=r(θ)r = r_*(\theta) — and Φ(θ,r(θ))\Phi(\theta, r_*(\theta)) is the characteristic point, since there the ray direction and the variation of the family become dependent. Off the envelope the map is a local diffeomorphism, and a point just inside the caustic is hit by two nearby rays (two solutions θ\theta), a point outside by none from that part of the family; on the caustic the two merge. Light intensity is inversely proportional to the Jacobian’s absolute value, so it blows up along the envelope: the caustic is the bright curve, brightest of all at the cusp, where the degeneracy is worst.

22. x=1costx' = 1 - \cos t, y=sinty' = \sin t, x=sintx'' = \sin t, y=costy'' = \cos t, so xyyx=cost1x'y'' - y'x'' = \cos t - 1 and γ2=2(1cost)=4sin2t2\norm{\gamma'}^2 = 2(1 - \cos t) = 4\sin^2\tfrac t2; by Proposition 18.17,

κ(t)=(1cost)8sin3t2=2sin2t28sin3t2=14sint2(0<t<2π).\kappa(t) = \frac{-(1 - \cos t)}{8\sin^3\tfrac t2} = \frac{-2\sin^2\tfrac t2}{8\sin^3\tfrac t2} = -\frac1{4\sin\tfrac t2} \qquad (0 < t < 2\pi).

With T=(sint2,cost2)T = (\sin\tfrac t2, \cos\tfrac t2) (divide γ\gamma' by 2sint22\sin\tfrac t2) and N=(cost2,sint2)N = (-\cos\tfrac t2, \sin\tfrac t2):

γ+1κN=γ4sint2(cost2, sint2)=(tsint+2sint, 1cost2(1cost)),\gamma + \frac1\kappa N = \gamma - 4\sin\tfrac t2\,\Bigl(-\cos\tfrac t2,\ \sin\tfrac t2\Bigr) = \bigl(t - \sin t + 2\sin t,\ 1 - \cos t - 2(1 - \cos t)\bigr),

i.e. c(t)=(t+sint, cost1)c(t) = (t + \sin t,\ \cos t - 1). Substituting t=u+πt = u + \pi:

c=(u+πsinu, cosu1)=((usinu)+π, (1cosu)2):c = \bigl(u + \pi - \sin u,\ -\cos u - 1\bigr) = \bigl((u - \sin u) + \pi,\ (1 - \cos u) - 2\bigr) :

the cycloid γ(u)\gamma(u) translated by (π,2)(\pi, -2). The evolute of a cycloid is a congruent cycloid, hanging one level below.

23. R(t)=1/κ(t)=4sint2R(t) = 1/\abs{\kappa(t)} = 4\sin\tfrac t2, so R(π)=4R(\pi) = 4: half of the arch length 88 computed in Exercise 18.1. The evolute’s velocity c(t)=(1+cost,sint)c'(t) = (1 + \cos t, -\sin t) vanishes at t=πt = \pi: the cusp is c(π)=(π,2)c(\pi) = (\pi, -2), directly below the apex γ(π)=(π,2)\gamma(\pi) = (\pi, 2), at distance 4=R(π)4 = R(\pi), exactly the length of the osculating radius there.

24. From Exercise 18.6, c(s)=κ(s)κ(s)2N(s)=R(s)N(s)c'(s) = -\frac{\kappa'(s)}{\kappa(s)^2}N(s) = R'(s)\,N(s), so c(s)=R(s)\norm{c'(s)} = \abs{R'(s)} and, for monotone RR,

s0s1c(s) ⁣ds=s0s1R(s) ⁣ds=R(s1)R(s0).\int_{s_0}^{s_1}\norm{c'(s)}\,\dd s = \Bigl|\int_{s_0}^{s_1}R'(s)\,\dd s\Bigr| = \abs{R(s_1) - R(s_0)} .

Say RR decreases. A string laid along the evolute beyond c(s0)c(s_0) and prolonged by the segment from c(s0)c(s_0) to γ(s0)\gamma(s_0) (which is tangent to the evolute, by question 4) has, when peeled off up to c(s)c(s) and stretched taut, straight part of length R(s0)(R(s0)R(s))=R(s)R(s_0) - \bigl(R(s_0) - R(s)\bigr) = R(s) pointing from c(s)c(s) along the normal — landing exactly on γ(s)\gamma(s): the free end traces the original curve (“involute”). Huygens hung a pendulum between two cycloidal cheeks: the cord wraps on the evolute, so the bob describes a cycloid — the tautochrone, whose oscillation period does not depend on the amplitude.

25. Tangents of the parabola: Δ=1\Delta = 1 and E=(1,t)0E' = (1, t) \neq 0 everywhere — smooth envelope (the parabola itself). Sliding ladder: Δ=1\Delta = -1, but E=32sin2t(cost,sint)E' = \tfrac32\sin 2t\,(-\cos t, \sin t) vanishes at the quadrant ends — the four cusps of the astroid. Normals of the parabola and of the cycloid: Δ=κ0\Delta = \kappa \neq 0, and E=RNE' = R'N vanishes exactly where the curvature is extremal — cusps of the evolutes at (0,1)(0,1) and (π,2)(\pi, -2). Caustic: Δ=2\Delta = 2, and E=32sinθ(cos2θ,sin2θ)E' = \tfrac32\sin\theta\,(\cos2\theta, \sin2\theta) vanishes at θ=0,π\theta = 0, \pi — the two cusps of the nephroid, the focal points of the mirror. Extrema of curvature must produce cusps on the envelope of the normals, since the evolute’s velocity is RNR'N: this is why the evolute of the ellipse has four cusps (four vertices), and the degenerate families (pencil, parallels) are the cases where the machine outputs a point or nothing at all.