I am given the variational principle: $$\delta \int_{\lambda 1}^{\lambda_2} \left(\frac{ds}{d\lambda}\right)^2d\lambda=0$$ where $ds^2=-dt^2+e^{2\Psi}dr^2+r^2(d\theta^2+\sin^2(\theta)d\phi^2)$ with $\Psi$ only depending on $r$.
I want to use this to write the geodesic equation for $r$ with respect to $\lambda$.
My first problem is interpreting $$\left(\frac{ds}{d\lambda}\right)^2$$ So $ds^2$ is the line element, where $ds^2=g_{\mu\nu}dx^\mu dx^\nu$ and since we are, I imagine, working with timelike separated paths, we then have: $$ds = \sqrt{-g_{\mu\nu}dx^\mu dx^\nu}$$ and hence I have: $$\left(\frac{\sqrt{-g_{\mu\nu}dx^\mu dx^\nu}}{d\lambda}\right)^2,$$ but what does this expression even mean?
Nothing seem to depend on $\lambda$? Furthermore I don't know what this has to do with the geodesic equation.