from __future__ import division
from visual import *
from visual.graph import *


funct1 = gcurve(color=color.cyan)

t=0
dt=0.001
x=10
v=0
tf=10
k=5
b=0

funct1.plot(pos=(t,x)) # Plot the initial point

while t<tf:
    a=-k*x - b*v    # Calculate a using the current x & v
    v=v+a*dt        # Use the acceleration to update the velocity
    x=x+v*dt        # Use the updated velocity to update the position
    t=t+dt          # Advance the time by a step dt
    funct1.plot(pos=(t,x)) # Plot the new point