I have now tested the simulation with various initial conditions, and the terminal velocity seem to be in accordance with the equations  :)