[step:Establish the chain homotopy identity $dK + Kd = \mathrm{id}_M^* - (i\circ p)^*$ via Cartan's Magic Formula]Fix $k \ge 1$ and $\omega \in \Omega^k(M)$. We claim
\begin{align*}
d(K\omega) + K(d\omega) = \mathrm{id}_M^* \omega - (i\circ p)^* \omega.
\end{align*}
Since exterior differentiation commutes with pullback and with integration against the parameter $t$ (the integration is over a parameter independent of $M$, and differentiation under the integral sign is justified by the smoothness of $H^*\omega$ on $M \times [0,1]$),
\begin{align*}
d(K\omega) = \int_0^1 d\bigl(j_t^*(\iota_{\partial_t} H^*\omega)\bigr) d\mathcal{L}^1(t) = \int_0^1 j_t^*\bigl(d\, \iota_{\partial_t} H^*\omega\bigr) d\mathcal{L}^1(t).
\end{align*}
Similarly, because $d \circ H^* = H^* \circ d$,
\begin{align*}
K(d\omega) = \int_0^1 j_t^*\bigl(\iota_{\partial_t} H^* d\omega\bigr) d\mathcal{L}^1(t) = \int_0^1 j_t^*\bigl(\iota_{\partial_t}\, d\, H^*\omega\bigr) d\mathcal{L}^1(t).
\end{align*}
Adding these and applying [Cartan's Magic Formula](/theorems/1535) — which states $\mathcal{L}_X = d\iota_X + \iota_X d$ for any smooth vector field $X$ on any smooth manifold — to $X = \partial_t$ on $M \times [0,1]$,
\begin{align*}
d(K\omega) + K(d\omega) = \int_0^1 j_t^*\bigl(\mathcal{L}_{\partial_t} H^*\omega\bigr) d\mathcal{L}^1(t).
\end{align*}
We now identify the integrand with $\frac{d}{dt}\bigl(j_t^* H^*\omega\bigr)$. For any $\eta \in \Omega^k(M \times [0,1])$ and any $x \in M$, $v_1, \dots, v_k \in T_xM$, the time-slice pullback $j_t^*\eta$ satisfies
\begin{align*}
\frac{d}{dt}\bigl(j_t^*\eta\bigr)_x(v_1, \dots, v_k) = \bigl(j_t^*\mathcal{L}_{\partial_t}\eta\bigr)_x(v_1, \dots, v_k),
\end{align*}
because the flow of $\partial_t$ on $M \times [0,1]$ is $\Phi_s(x, t) = (x, t+s)$ (defined on the appropriate open subset of $M \times \mathbb{R}$ containing the orbits we use), so $j_t \circ \Phi_s|_{M \times \{0\}} = j_{t+s}$ in a neighbourhood of any $t \in (0,1)$, and the Lie derivative is by definition $\mathcal{L}_{\partial_t}\eta = \frac{d}{ds}\big|_{s=0} \Phi_s^*\eta$, hence $j_t^*\mathcal{L}_{\partial_t}\eta = \frac{d}{ds}\big|_{s=0} j_t^*\Phi_s^*\eta = \frac{d}{ds}\big|_{s=0} j_{t+s}^*\eta$. At the endpoints $t = 0$ and $t = 1$, one-sided derivatives suffice and the formula extends by continuity since the integrand is smooth in $t$.
Applying the [fundamental theorem of calculus](/theorems/632) to the smooth $t$-family of $k$-forms $t \mapsto j_t^* H^*\omega$,
\begin{align*}
d(K\omega) + K(d\omega) = \int_0^1 \frac{d}{dt}\bigl(j_t^* H^*\omega\bigr) d\mathcal{L}^1(t) = j_1^* H^*\omega - j_0^* H^*\omega.
\end{align*}
By the contravariance of pullback applied to the identities $H \circ j_1 = \mathrm{id}_M$ and $H \circ j_0 = i \circ p$,
\begin{align*}
j_1^* H^*\omega = (H \circ j_1)^*\omega = \mathrm{id}_M^*\omega = \omega, \qquad j_0^* H^*\omega = (H \circ j_0)^*\omega = (i\circ p)^*\omega.
\end{align*}
Combining,
\begin{align*}
d(K\omega) + K(d\omega) = \omega - (i \circ p)^*\omega.
\end{align*}[/step]