Rocketsim [patched] -

[ m(t) \cdot \fracd\vecvdt = \vecF thrust + \vecF drag + \vecF gravity + \vecF lift ]