If $X$ is a smooth vector field on $U$ and $\omega$ is a smooth differential form, then $i_{X}\omega$ is smooth
The above is an Exercise 14.22 in Lee's introduction to smooth manifolds.
Recall the definition of smooth differential form: A smooth differential form on $U$ is a map $\omega:U\to \Lambda^k(T^*U)$ such that $\omega_p\in \Lambda^k(T^*_pU)$ and $\omega_p = \sum_{I}'\omega_I(p)dx^I|_p$ where $w_I:U\to \Bbb R$ are $\textbf{smooth}$ functions.
Also, if a smooth vector field on $U$ is a map $X:U\to TU$ such that $X_p\in T_pU$ and $X_p = \sum_{i=1}^n X^i(p)\frac{\partial}{\partial x^i}|_p$ where $X^i:U\to\Bbb R$ are smooth functions.
Here, $(i_{X}\omega)_p:= i_{X_p}\omega_p$. So I need to transform this $i_{X_p}\omega_p$ to the above form to check whether the coefficient function is smooth or not. But how can I transform? Any help?
Edit: For simplicity, I changed my question to an open subset $U\subset\Bbb R^n$ not $M$.
Edit: If $X = \sum_{i=1}^n X^i \frac{\partial}{\partial x^i}$ and $\omega_p = \sum_{I}' \omega_I dx^I$ then by multilinearity, $i_Xw = w(X,...) = w(\sum_{j}X^j\frac{\partial}{\partial x^j},...) = \sum_{j}X^j\omega(\frac{\partial}{\partial x^j},...) = \sum_jX^j(\sum_I'\omega_I dx^I)(\frac{\partial}{\partial x^j},...) = \sum_j\sum_I'X^jw_I dx^I(\frac{\partial}{\partial x^j},...) =\sum_j\sum_I'X^jw_I (\frac{\partial}{\partial x^j}\lrcorner dx^I)$