Let (M, g)$(M, g)$ be a smooth Riemannian manifold, p∊M$p \in M$, and exp-P$\exp_P$ the exponential map at the point P$P$: expP: T(P)⟶M$\exp_P: T(P) \to M$
It seems clear to me that exp-P$\exp_P$ is smooth on U/{0}$U \setminus \{0\}$, where U$U$ is a neighborhood of the origin 0∊ T(P)$0 \in T(P)$, because exp-P$\exp_P$ is defined from the geodesics, which are solutions of a system of ordinary differential equations.
But is it true that exp-P$\exp_P$ is also always smooth at the origin 0∊T(P)$0 \in T(P)$?
I am not interested for the moment on the ways to calculate these high order derivatives, but just on the theoretical question about their existence.
Thank you!