Чего здесь не хватает, так это того, что вам нужно учитывать отклонения от начальной точки, а не x = 0 (а также учитывать также знак направления, что правильно указывают другие). То есть, если ваша начальная точка x0, ваше уравнение должно быть больше похоже на:
x += (x-x0)/5
Вот рисунок для движения в положительном и отрицательном направлениях (обратите внимание, что положение находится на вертикальной оси, а время на горизонтальной)
А вот и код Python. (Обратите внимание, что я добавил термин dt, так как слишком странно проводить динамическое моделирование без явного времени.)
from pylab import *
x0, b, dt = 11.5, 5, .1
xmotion, times = [], []
for direction in (+1, -1):
x, t = x0+direction*dt/b, 0 # give system an initial kick in the direction it should move
for i in range(360):
x += dt*(x-x0)/b
t += dt
xmotion.append(x)
times.append(t)
plot(times, xmotion, '.')
xlabel('time (seconds)')
ylabel('x-position')
show()