For near , let be the unique geodesic with and . The exponential map is .
The constant initial velocity gives . Varying the initial velocity through yields , so
Thus is the identity. The inverse function theorem makes a diffeomorphism from a neighborhood of onto a neighborhood of . Coordinates from an orthonormal basis of are the geodesic normal coordinates.