[step:Solve the tangent ODE inside $G$ and compare it with the matrix exponential]
Let $TG:=\bigsqcup_{A\in G}T_A G$ denote the tangent bundle of $G$. Define the vector field
\begin{align*}
W:G\to TG,
\end{align*}
by sending each $A\in G$ to the tangent vector $AV\in T_A G$. Since $G$ is an embedded submanifold of the finite-dimensional real vector space $M(n,\mathbb C)$, $G$ is a finite-dimensional smooth manifold. The map $A\mapsto AV$ is a smooth real-linear map from $M(n,\mathbb C)$ to $M(n,\mathbb C)$, and its restriction to $G$ is smooth. The previous step shows that this restricted map has value in $T_A G$ at each $A\in G$, so $W$ is a smooth tangent vector field on the smooth manifold $G$. By the local existence theorem for smooth ordinary differential equations on finite-dimensional smooth manifolds, applied to the smooth vector field $W$ with initial point $I\in G$, there exist $\varepsilon>0$ and a smooth curve
\begin{align*}
\beta:(-\varepsilon,\varepsilon)\to G
\end{align*}
such that
\begin{align*}
\beta(0)=I
\end{align*}
and
\begin{align*}
\beta'(t)=\beta(t)V
\end{align*}
for all $t\in(-\varepsilon,\varepsilon)$.
Viewed as a curve in the ambient vector space $M(n,\mathbb C)$, the same curve solves the matrix [initial value problem](/page/Initial%20Value%20Problem)
\begin{align*}
Y'(t)=Y(t)V,\qquad Y(0)=I.
\end{align*}
The curve
\begin{align*}
\alpha:\mathbb R\to M(n,\mathbb C),\quad t\mapsto \exp(tV)
\end{align*}
also solves this initial value problem. Indeed, by [citetheorem:8770], $\alpha$ is a smooth one-parameter subgroup of $GL(n,\mathbb C)$. The expansion at $0$ obtained from the absolutely convergent matrix exponential series [citetheorem:8777] gives $\alpha'(0)=V$. For each $t\in\mathbb R$, the homomorphism identity gives $\alpha(t+h)=\alpha(t)\alpha(h)$, so differentiating with respect to $h$ at $h=0$ yields
\begin{align*}
\alpha'(t)=\alpha(t)V=\exp(tV)V.
\end{align*}
Also $\alpha(0)=I$. The right-hand side $Y\mapsto YV$ is a linear, hence smooth and locally Lipschitz, vector field on the finite-dimensional real vector space $M(n,\mathbb C)$. Therefore uniqueness for linear ordinary differential equations in finite-dimensional real vector spaces applies, and
\begin{align*}
\beta(t)=\exp(tV)
\end{align*}
for every $t\in(-\varepsilon,\varepsilon)$. Since $\beta(t)\in G$, we have
\begin{align*}
\exp(tV)\in G
\end{align*}
for every $t\in(-\varepsilon,\varepsilon)$.
[/step]