I wonder if there is a way to reproduce the 4-force generalization for Newton's equation for a free particle i.e.
$$ m\frac{d^2x^\mu}{d \tau^2} = 0, \qquad \text{ for } \, \mu =0,1,2,3, \tag{1} $$
with just calculus of variations? If I start out with the action
$$ S = -mc^2\int\frac{dt}{\gamma}\tag{2} $$
with
$$ L =- \frac{mc^2}{\gamma},\qquad \gamma:= \frac{1}{\sqrt{1-\frac{{\bf v}^2}{c^2}}},\tag{3} $$
and put this in Euler-Lagrange equation I obtain
$$ \frac{d}{dt}\left( \frac{\partial L}{\partial v^i} \right) - \left(\frac{\partial L}{\partial x^i}\right) = 0 $$ $$\qquad\Leftrightarrow\qquad \frac{d}{dt}\left( \gamma m v^i\right)= \frac{d}{dt}p^i= 0 \qquad \text{ for } \, i =1,2,3. \tag{4} $$
I believe this is done right? But what about the case
$$ m\frac{d^2x^0}{d \tau^2} = 0,\tag{5} $$
i.e. the case for the temporal index $\mu = 0$?