if you are only every gonna use critical damping, you can do something like: function crit_response(dt,pos,vel,rate) local dissipation = math.exp(-d

if you are only every gonna use critical damping, you can do something like: function crit_response(dt,pos,vel,rate) local dissipation = math.exp(-dt*rate) local disspatePos = pos*dissiptation local disspateVel = vel*dissiptation local posCarried = 1 + dt*rate local velCarried = 1 - dt*rate local posFromVel = dt local velFromPos = -dt*rate*rate local newpos = posCarried*dissipatePos + posFromVel*dissipateVel local newvel = velFromPos*dissipatePos + velCarried*dissipateVel return newpos , newvel end (intentionally verbose for self-didactic purposes to show how you can actually split the dissipation step into its own loop out from the elastic calculation, meaning that you can use this as a more accurate initial guess in numeric simulations compared to just constant-acceleration approximations)