D'Alembert equation reads, for the considered case,
$$\mu\frac{\partial^2 y}{\partial t^2} = T_0 \frac{\partial^2 y}{\partial x^2}\:.\tag{1}$$
It is nothing but $F=ma$ along the vertical direction ($y$). Here $y$ denotes the small deformation of the string along the vertical direction from the stationary (horizontal) configuration. We disregard horizontal deformations. $T_0 \frac{\partial^2 y}{\partial x^2}$ is the $x$-derivative of vertical component of force acting on a infinitesimal segment of string $dx$, when the deformation is small, assuming that the absolute value of the tension $T_0$ is constant. In other words
$$f_y = T_0 \frac{\partial y}{\partial x}\tag{2}$$
is the $y$ component of either the tension or the force acting at the endpoints.
From (1) and the defintion of $e(t,x)$
$$e(t,x) = \frac{1}{2}\left[T_0\left(\frac{\partial y}{\partial x}\right)^2 + \mu\left(\frac{\partial y}{\partial t}\right)^2\right]$$
one easily sees (just computing derivatives) that
$$\frac{\partial e}{\partial t} = \frac{\partial}{\partial x} \left(\frac{\partial y}{\partial t} T_0\frac{\partial y}{\partial x}\right)\:.\tag{3}$$
Integrating this identity along the string, from the point (maybe endpoint) $x_1$ to the point (maybe endpoint) endpoint $x_2$, we have
$$\frac{d}{dt}\int_{x_1}^{x_2} e(t,x) dx = \left.\left(\frac{\partial y}{\partial t}T_0\frac{\partial y}{\partial x}\right)\right|_{x_2} -
\left.\left(\frac{\partial y}{\partial t}T_0\frac{\partial y}{\partial x}\right)\right|_{x_1}\tag{4}$$
As $\frac{\partial y}{\partial t}$ is the vertical velocity of the endpoint, taking (2) into account, we conclude that $\frac{\partial y}{\partial t}T_0\frac{\partial y}{\partial x}$ is the power of the force acting on the endpoint (or at a generic point of the string).
Equations (3) and (4) establish relations, respectively local and global, between the density of energy and the power of the forces acting along the string.