  Let $r$ be a positive integer, $N$ be a nonnegative integer and $\Omega
\subset \mathbb{R}^{r}$ be a domain. Further, for all multi-indices $\alpha \in
\mathbb{N}^{r}$, $|\alpha|\leq N$, let us consider the partial differential
operator $D^{\alpha}$ defined by \[
  D^{\alpha}= \frac{\partial^{|\alpha|}}{\partial x_{1}^{\alpha_{1}}\cdots
\partial x_{r}^{\alpha_{r}}}, \] where $\alpha= (\alpha_{1}, \ldots,
\alpha_{r})$. Here by definition we mean $D^{0}\equiv \mathrm{id}$. An easy
computation shows that if $f, g\in \mathscr{C}^{N}(\Omega)$ and $\alpha \in
\mathbb{N}^{r}, |\alpha|\leq N$, then we have \[ \tag{$\ast$} D^{\alpha}(f\cdot
g) = \sum_{\beta\leq \alpha}\binom{\alpha}{\beta}D^{\beta}(f)\cdot D^{\alpha -
\beta}(g). \] This paper is devoted to the study of identity $(\ast)$ in the
space $\mathscr{C}(\Omega)$. More precisely, if $r$ is a positive integer, $N$
is a nonnegative integer and $\Omega \subset \mathbb{R}^{r}$ is a domain, then
we describe those mappings $T_{\alpha} \colon \mathscr{C}(\Omega)\to
\mathscr{C}(\Omega)$, $\alpha \in \mathbb{N}^{r}, |\alpha|\leq N$ that satisfy
identity $(\ast)$ for all possible multi-indices $\alpha\in \mathbb{N}^{r}$,
$|\alpha|\leq N$. Our main result says that if the domain is
$\mathscr{C}(\Omega)$, then the mappings $T_{\alpha}$ are of a rather special
form. Related results in the space $\mathscr{C}^{N}(\Omega)$ are also
presented.
