\frac{d}{dt}(\Delta m_1\vec{v_1}+\Delta m_2\vec{v_2})=(\Delta m_1+\Delta m_2)\vec{f}+\oint_{\Delta S_1+\Delta S_2}\vec{P_n}dS