  Let $\mathcal{H}$ be the class of all analytic self-maps of the open unit
disk $\mathbb{D}$. Denote by $H^n f(z)$ the $n$-th order hyperbolic derivative
of $f\in \mathcal H$ at $z\in \mathbb{D}$. For $z_0\in \mathbb{D}$ and $\gamma
= (\gamma_0, \gamma_1 , \ldots , \gamma_{n-1}) \in {\mathbb D}^{n}$, let
${\mathcal H} (\gamma) = \{f \in {\mathcal H} : f (z_0) = \gamma_0,H^1f (z_0) =
\gamma_1,\ldots ,H^{n-1}f (z_0) = \gamma_{n-1} \}$. In this paper, we determine
the variability region $V(z_0, \gamma ) = \{ f^{(n)}(z_0) : f \in {\mathcal H}
(\gamma) \}$, which can be called ``the generalized Schwarz-Pick Lemma of
$n$-th derivative". We then apply the generalized Schwarz-Pick Lemma to
establish a $n$-th order Dieudonn\'e's Lemma, which provides an explicit
description of the variability region $\{h^{(n)}(z_0): h\in \mathcal{H},
h(0)=0,h(z_0) =w_0, h'(z_0)=w_1,\ldots, h^{(n-1)}(z_0)=w_{n-1}\}$ for given
$z_0$, $w_0$, $w_1,\dots,w_{n-1}$. Moreover, we determine the form of all
extremal functions.
