We study Einstein's equations with an isotropic but inhomogeneous metric in the cosmic rest frame. The equations are solved perturbatively in the late Universe. The leading plus next-to-leading order results agree with observations without using a cosmological constant or dark matter.

Considering radial geodesics in the Robertson-Walker metric leads us to abandon the co-moving coordinates. Instead we work in the cosmic rest frame. Since then the matter is in motion, the solution of Einstein's equations is more complicated. We calculate the first correction to standard cosmology which has an off-diagonal term b dt dr in the metric. It describes the late universe. We then solve Maxwell's equations in the new metric and discuss redshift and luminosities. We obtain the correct age of the universe T=14 Gyr= 1/H, without assuming a cosmological constant.

We continue to study the spherically symmetric vacuum solutions of the tensor-four-scalars theory which is a modification of general relativity. Our aim is to construct vacuum solutions with asymptotically constant circular velocity curves which can be used to represent the dark halo outside the baryonic matter. These solutions are obtained by means of a new non-standard gauge in the field equations. The metric is expressed in terms of the rotation curve $V(r)$ in closed form.

Massive gravity previously constructed as the spin-2 quantum gauge theory leads in the mass zero limit to a modification of general relativity. As a relic from the massive theory a vector field v survives which couples to the metric only. However, the coupling does not involve covariant derivatives so that v must be considered as 4 scalar fields. We analyse the static, spherically symmetric solutions of this theory. From the corresponding geodesics we find the circular velocity profile. Interpreting this as coming from a dark density profile, the theory predicts a flat density core for r going to 0. But the dark density profile is not universal.

Using Weierstrassian elliptic functions the exact geodesics in the Schwarzschild metric are expressed in a simple and most transparent form. The results are useful for analytical and numerical applications. For example we calculate the perihelion precession and the light deflection in the post-Einsteinian approximation. The bounded orbits are computed in the post-Newtonian order. As a topical application we calculate the gravitational red shift for a star moving in the Schwarzschild field.

Based on the geodesic equation in a static spherically symmetric metric we discuss the rotation curve and gravitational lensing. The rotation curve determines one function in the metric without assuming Einstein's equations. Then lensing is considered in the weak field approximation of general relativity. From the null geodesics we derive the lensing equation and corrections to it.