Kerr metric geodesics using EL equations

Let there exist a Kerr metric $$ds^2$$ in Boyer–Lindquist coordinates such that $$ds^2=-\left(1-\frac{r_sr}{\Sigma}\right)dt^2+\left(\frac{\Sigma}{\Delta}\right)dr^2+\Sigma d \theta^2+\left(r^2+a^2+\frac{r_sra^2}{\Sigma}\sin^2(\theta)\right)\sin^2(\theta)d \phi^2-\left(\frac{2r_sra\sin^2(\theta)}{\Sigma}\right)dt d \phi,$$ such that $$a=\frac{J}{Mc}$$ $$\Sigma=r^2+a^2\cos^2(\theta)$$ $$\Delta=r^2-r_sr+a^2$$ $$r_s=2M$$ $$L=\frac{1}{2}\left[-\left(1-\frac{r_sr}{\Sigma}\right)t'^2+\left(\frac{\Sigma}{\Delta}\right)r'^2+\Sigma \theta'^2+\left(r^2+a^2+\frac{r_sra^2}{\Sigma}\sin^2(\theta)\right)\sin^2(\theta)\phi'^2-\left(\frac{2r_sra\sin^2(\theta)}{\Sigma}\right)t' \phi'\right]$$ where we are letting $$c=G=1$$hence the equation of $$ds^2$$. To find the geodesics of the Kerr metric $$ds^2$$ we need to apply the Euler-Lagrange equations $$\frac{\partial L}{\partial x^\mu}-\frac{d}{d \tau}\left(\frac{\partial L}{\partial x'^\mu}\right)=0$$ where $$x^\mu=(x^0,x^1,x^2,x^3)=(t,r,\theta,\phi)$$ and $$\forall x'^\mu=\frac{d x^\mu}{d \tau}.$$ In order to obtain the system of ODE geodesic equations for the Kerr metric will I just solve $$\frac{\partial L}{\partial t},\frac{\partial L}{\partial t'},\frac{\partial L}{\partial r},\frac{\partial L}{\partial r'},\frac{\partial L}{\partial \theta},\frac{\partial L}{\partial \theta'},\frac{\partial L}{\partial \phi}, \frac{\partial L}{\partial \phi'}$$ then just plug them into the equations $$\frac{\partial L}{\partial t}-\frac{d}{d \tau}\left(\frac{\partial L}{\partial t'}\right)=0$$ $$\frac{\partial L}{\partial r}-\frac{d}{d \tau}\left(\frac{\partial L}{\partial r'}\right)=0$$ $$\frac{\partial L}{\partial \theta}-\frac{d}{d \tau}\left(\frac{\partial L}{\partial \theta'}\right)=0$$ $$\frac{\partial L}{\partial \phi}-\frac{d}{d \tau}\left(\frac{\partial L}{\partial \phi'}\right)=0$$ and numerically solve or is there something I am missing?

• @m4r35n357 Wouldnt we arive at the trajectory equations by just applying the Euler-Lagrange equations on the Lagrangian $L$ and then using constants to manipulate the results?