I am trying to prove that, if $\rho:S^n\times S^n\to \Bbb{R}$ is the distance induced by the Riemannian metric on the sphere $S^n\subset \Bbb{R}^{n+1}$, then $$\rho(p,q)=\arccos \langle p,q\rangle,\,\,\,\forall p,q\in S^n$$ with $\arccos$ defined from $[-1,1]$ to $[0,\pi]$.
If $p=q$, then the formula is trivial. We have then two cases: $p=-q$ and $p\neq- q$. I am having trouble even with the apparently easier case, $p=-q$.
If $p=-q$, then let $v\in S^n$ be any vector orthogonal to $p$. Then $\alpha:[0,\pi]\to S^n$, $\alpha(t)=(\cos t)p+(\sin t)v$ is a well defined differentiable path from $p$ to $q$ such that $\ell_0^\pi(\alpha)=\pi$ ("length of $\alpha$"). This guarantees that $\rho(p,q)\leq \pi (=\arccos \langle p,q\rangle$, in this case).
In order to show that $\rho(p,q)=\pi$, I must consider an arbitrary differentiable by parts path $\beta:[a,b]\to S^n$ from $p$ to $q$ and show that $$\pi\leq \ell_a^b(\beta)=\int_a^b|\beta'(t)|\,dt.$$
I've done some geometric observations and computations, but without success. How can I do this? I wish I could do it without using any facts about geodesics (they appear later in the book I'm studying with).
Let $p=(1,0,\dots,0)$ and $q=(-1,0,\dots,0)$. Then any path from $p$ to $q$ is of the form $\gamma(t)=(t,h(t))$, where $h(t)\in \Bbb{R}^n$. Then the length of the path is $$ L=\int_{-1}^{1}|\gamma'(t)|dt. $$ But $\gamma'(t)=(1,h'(t))$, and so $|\gamma'(t)|=\sqrt{1+|h'(t)|^2}$.
Now $t^2+\langle h(t),h(t)\rangle=\langle \gamma(t),\gamma(t)\rangle=1$, so we get $$ 2t+2\langle h(t),h'(t)\rangle=0 $$ by differentiating both sides with respect to t.
By Cauchy Schwartz we have $|h'(t)|^2|h(t)|^2\ge |\langle h(t),h'(t)\rangle|^2=|t|^2$, and so $$ 1+|h'(t)|^2\ge \frac{t^2}{|h(t)|^2}+1 $$
but $|h(t)|^2=1-t^2$, hence $$ \sqrt{1+|h'(t)|^2}\ge \sqrt{1+\frac{t^2}{1-t^2}}=\sqrt{\frac{1}{1-t^2}}, $$ and so $$ L=\int_{-1}^{1}\sqrt{1+|h'(t)|^2}dt\ge\int_{-1}^1 \sqrt{\frac{1}{1-t^2}}dt=\pi. $$
A similar reasoning applies for any other point: Assume $p=(1,0,\dots,0)$, then you can assume that the other point $q$ is of the form $(a,s,0,\dots,0)$. The path $\gamma(t)=(t,\sqrt{1-t^2},0,\dots,0)$ has length $$ L=\int_{a}^{1}|\gamma'(t)|dt=\int_{a}^1 \sqrt{\frac{1}{1-t^2}}dt =-arccos(t)|_{a}^1=arccos(a)=arccos(\langle p,q\rangle). $$ By a similar argument as above, any other path has the form $(t,h(t))$ with $\sqrt{1+|h'(t)|^2}\ge \sqrt{\frac{1}{1-t^2}}$, so the result follows.