Here $\omega$ is a smooth form on a manifold, $X, Y_1, \ldots, Y_p$ are smooth vector fields, and $L_X$ represents the Lie derivative.
I am having trouble proving this standard identity. All the books that I have checked in describe it as a trivial exercise, but I would appreciate a hint or a sketch of an argument. Even in the reduction to checking for $\omega = f \, dx_1 \vee dx_2$ in a local coordinate system I don't see an obvious way forward.
$$L_X( \omega (Y_1,\ldots, Y_n)) = (L_X \omega) (Y_1, \ldots, Y_n) + \sum_{i = 1}^n \omega (Y_1, \ldots, Y_{i-1}, L_X Y_i, \ldots, Y_n)$$
It is "trivial" - work in a local coordinate system and apply the usual tricks for computing the derivative of a product. Take advantage of the multilinearity of forms and the existence of a local rectification for $X$ (at nonzero points) to use the canonical identification between the tangent spaces on $R^n$.