F=m\frac{dv}{dt}+v_R\frac{dm}{dt}