function dvdx = bumper(x, v)
global K m
dvdx = -(K*v^2*(x + 1)^3)/m;