The theory of ordinary differential equations shows that if is smooth then the differential equation with initial condition has a unique solution for sufficiently small and the solution depends smoothly on and .
There is a standard technique (see for example ) for computing the change of variables to normal coordinates , at a point as a formal Taylor series expansion.