Solution (source code)

= Solution

For an <embedded submanifold> $X\subseteq\mathbb R^N$ and $p\in X$, the <tangent space> is
$$
\boxed{
T_pX=\{\gamma'(0):
\gamma:(-\varepsilon,\varepsilon)\to X
\text{ is smooth and }\gamma(0)=p\}.
}
$$
Equivalently, in a local parametrization $\varphi:U\subseteq\mathbb R^d\to X$ with $\varphi(u)=p$, one has $T_pX=\operatorname{im}D\varphi_u$.