My two favorite ways:
Alternative A. Assuming that you are familiar with the Taylor series of both sine and cosine (some purists think that you should not rely on pictures when defining these functions, and power series are the way to go).
Consider the matrix
$$
A=\left(\begin{array}{cc}0&-1\\1&0\end{array}\right).
$$
We can use this as an infinitesimal generator of rotations. IOW it is then easy see, using the Taylor series of $e^x$, that we have
$$
e^{xA}=\sum_{n=0}^\infty\frac{x^nA^n}{n!}=\left(\begin{array}{cc}\cos x&-\sin x\\\sin x&\cos x\end{array}\right)
$$
that we recognize as the familiar rotation matrix $R(x)$.
Differentiating this gives, using the chain rule and the known derivative of $e^x$,
$$
\frac d{dx}e^{Ax}=Ae^{Ax}.
$$
Using the rotation matrix representation allows us to rewrite this in the form
$$
\frac d{dx} R(x)=A R(x)=\left(\begin{array}{cc}-\sin x&-\cos x\\\cos x&-\sin x\end{array}\right).
$$
The derivatives of the matrix entries can be read from this. We get $D\cos x$ and $D\sin x$ from two entries actually.
Alternative B. Here I'm using the unit circle and arcs along it in the usual way for the definitions of sine and cosine.
Assume that we have already travelled the distance $\alpha$ along the unit circle in ccw direction, starting from $(1,0)$ obviously. Assume that we have reached the point $P=(x,y)$. We know that $x=\cos\alpha$ and
$y=\sin\alpha$. The question I want you to ask yourself is: "which direction are we heading now?" A moment's reflection reveals that we are heading in the direction $\vec{t}$ of the tangent of the unit circle drawn at $P$. That vector $\vec{t}$ is thus orthogonal to the vector $(x,y)$. This allows us to conclude that $\vec{t}$ is parallel to $(-y,x)$. Furthermore, given that $P$ was assumed to be on the unit circle we also know that $||(-y,x)||=\sqrt{x^2+y^2}=1$. By drawing a picture it is clear that the vector $(-y,x)$, drawn starting from $P$, points in the direction of increasing $\alpha$.
Hence we can use $\vec{t}=(-y,x)$. So if we move along the circle an infinitesimal distance $h$ we move from point $P=(x,y)$ to the point
$$P'=(x,y)+h\vec{t}=(x-hy,y+hx)=(\cos(\alpha+h),\sin(\alpha+h)).$$
(Ok, this is really just a good approximation. See the last paragraph for a better way of justifying this.)
Hence
$$
\lim_{h\to0}\frac{\cos(\alpha+h)-\cos\alpha}h=\lim_{h\to0}\frac{-hy}h=-y=-\sin\alpha,
$$
and similarly
$$
\lim_{h\to0}\frac{\sin(\alpha+h)-\sin\alpha}h=\lim_{h\to0}\frac{hx}h=x=\cos\alpha.
$$
Alternative B probably also feels a bit hand-wavy. It is equivalent to the text book derivation of the formula for the tangent vector of a parametrized plane curve (this time the curve parametrized by $(x,y)=(\cos\alpha,\sin\alpha)$. Using the natural parameter (=arc length $s$), we always get a tangent vector that has unit length. And this time I won't differentiate the component functions but I use the known (unit!) tangent direction along a circle instead :-)