Let X be a smooth manifold and $v, w \in \Gamma(TX)$ be two vector fields. Under the natural correspondence of vector fields and derivations on $C^\infty(X)$ let $\delta, \epsilon: C^\infty(X) \rightarrow C^\infty(X)$ be the corresponding derivations.
Then the commutator $[\delta, \epsilon] = \delta \circ \epsilon - \epsilon \circ \delta$ is another derivation. Define the Lie bracket $[v,w]$ as the vector field corresponding to this commutator. Using local coordinates on X, we can express the action of the two derivations on some $a \in C^\infty(X)$ as $$\delta(a) = \sum v_i \frac{\partial a}{\partial x_i} \ \text{ and } \ \epsilon(a) = \sum w_i \frac{\partial a}{\partial x_i}$$
Using this how can $[v,w]$ be expressed using the same local coordinates? I suspect the answer is the following, but I don't know how to derive it.
$$[v,w] = \sum_{i,j} \bigg(v_i \frac{\partial w_j}{\partial x_i}- w_i \frac{\partial v_j}{\partial x_i}\bigg) \frac{\partial}{\partial x_j}$$
Associating the vector field $v_i$ with the differential operator $v:=v_i\partial^i$ gives $vwf=v_i\partial^i(w_j\partial^jf)$ for a sufficiently nice function $f$. The product rule, plus the symmetry of derivatives to delete the second-order terms, gives $[v,\,w]f=Rf$, where $R$ is the expression for $[v,\,w]$ you desire. The rest is an exercise.