[step:Implication $(5) \implies (1)$ via extending an inextensible geodesic]Assume $(M, d)$ is complete. Let $\gamma : [0, b) \to M$ be a maximal geodesic with $b < \infty$; we derive a contradiction. Without loss of generality $\gamma$ has constant speed $|\dot\gamma|_g \equiv c \geq 0$ by [Geodesics Have Constant Speed](/theorems/2709). If $c = 0$, then $\gamma$ is constant: the geodesic equation $\nabla_{\dot\gamma}\dot\gamma = 0$ with $\dot\gamma \equiv 0$ admits the constant curve as a global solution defined on all of $\mathbb{R}$, so the maximal domain is $\mathbb{R}$ already, contradicting $b < \infty$. Hence assume $c > 0$.
Pick any sequence $t_n \uparrow b$ with $t_n \in [0, b)$. For $m, n$ with $t_m, t_n$ close to $b$,
\begin{align*}
d(\gamma(t_m), \gamma(t_n)) \leq \ell(\gamma|_{[\min(t_m, t_n), \max(t_m, t_n)]}) = c \cdot |t_m - t_n|.
\end{align*}
Since $(t_n)$ is Cauchy in $[0, b)$ (it converges to $b$), $(\gamma(t_n))$ is Cauchy in $M$. By $(5)$, there exists $q \in M$ with $\gamma(t_n) \to q$.
The limit $q$ does not depend on the sequence $(t_n)$: if $(t_n')$ is another sequence with $t_n' \uparrow b$, form the merged sequence $(s_k)$ defined by $s_{2k-1} := t_k$ and $s_{2k} := t_k'$. Both $(t_n)$ and $(t_n')$ converge to $b$, so $(s_k) \to b$, and in particular $(s_k)$ is Cauchy in $\mathbb{R}$. By the Lipschitz estimate from constant speed, $d(\gamma(s_j), \gamma(s_k)) \leq c\, |s_j - s_k|$, so $(\gamma(s_k))$ is Cauchy in $(M, d)$ and converges by (5) to some single limit $q^* \in M$. Both subsequences $(\gamma(t_n))$ and $(\gamma(t_n'))$ inherit this limit, forcing the two limits to coincide. Hence $\gamma$ extends continuously to $\bar\gamma : [0, b] \to M$ with $\bar\gamma(b) = q$.
We invoke the following uniform local existence statement for the geodesic flow.
[claim:Uniform local existence on a neighbourhood of $q$]
There exist $\rho, \varepsilon > 0$ such that for every $r \in M$ with $d(r, q) < \rho$ and every $w \in T_r M$ with $|w|_g \leq 1$, the geodesic $\gamma_w$ with $\gamma_w(0) = r$ and $\dot\gamma_w(0) = w$ is defined on the entire interval $[-\varepsilon, \varepsilon]$.
[/claim]
[proof]
The geodesic spray $G$ is a smooth vector field on $TM$: the geodesic equation $\nabla_{\dot\gamma}\dot\gamma = 0$, second-order on $M$, is the first-order ODE $\dot{(r,w)} = G(r,w)$ on $TM$. By the standard theorem on smooth flows of smooth vector fields, the maximal flow domain $\mathcal{D} \subseteq \mathbb{R} \times TM$ of $G$ is open in $\mathbb{R} \times TM$, and the flow $\Phi : \mathcal{D} \to TM$ is smooth.
Fix $\rho_0 > 0$ small enough that $\overline{B}_d(q, \rho_0)$ is compact (such $\rho_0$ exists because, by [Exponential Map as a Local Diffeomorphism](/theorems/2712) applied at $q$, there is a normal neighbourhood of $q$ inside which sufficiently small closed metric balls have compact closure). The closed unit-disk bundle restricted to this ball,
\begin{align*}
K := \{(r, w) \in TM : r \in \overline{B}_d(q, \rho_0),\ |w|_g \leq 1\},
\end{align*}
is compact in $TM$ (a fibre bundle with compact fibre $\overline{B}(0, 1) \subseteq \mathbb{R}^n$ over the compact base).
The set $\{0\} \times K \subseteq \mathcal{D}$ since the flow is defined at time $0$ from any initial condition. By openness of $\mathcal{D}$, for each $(r, w) \in K$ there exist $\varepsilon_{(r,w)} > 0$ and an open neighbourhood $U_{(r,w)} \subseteq TM$ of $(r,w)$ such that $(-\varepsilon_{(r,w)}, \varepsilon_{(r,w)}) \times U_{(r,w)} \subseteq \mathcal{D}$. The collection $\{U_{(r,w)} : (r,w) \in K\}$ is an open cover of the compact set $K$; extract a finite subcover indexed by $(r_1, w_1), \dots, (r_N, w_N)$ and set $\varepsilon := \min_{1 \leq i \leq N} \varepsilon_{(r_i, w_i)} > 0$. Then $(-\varepsilon, \varepsilon) \times K \subseteq \mathcal{D}$. Choosing any $\rho \in (0, \rho_0)$ gives the claim.
[/proof]
Apply the claim with the unit vector $w := \dot\gamma(t_0) / c$ (recall $c = |\dot\gamma|_g > 0$): for any $r$ with $d(r, q) < \rho$, the geodesic at $r$ in any unit direction exists for time at least $\varepsilon$. By [Geodesic Rescaling](/theorems/2710), the geodesic at $r$ with initial velocity of norm $c$ exists for time at least $\varepsilon / c$.
Choose $t_0 \in [0, b)$ close enough to $b$ that $d(\gamma(t_0), q) < \rho$ and $b - t_0 < \varepsilon / c$. Then the geodesic with initial conditions $(\gamma(t_0), \dot\gamma(t_0))$ — which is $\gamma$ itself by uniqueness — extends for time at least $\varepsilon / c > b - t_0$ past $t_0$. This means $\gamma$ extends past $b$, contradicting maximality of the domain $[0, b)$.
By the symmetric argument with $\gamma : (a, 0] \to M$ (or equivalently reversing the parametrisation), the geodesic also extends past any finite left endpoint. Hence every maximal geodesic has domain $\mathbb{R}$, proving $(1)$.[/step]