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


#objects
funct1 = gcurve(color=color.cyan)

#constants and initial conditions
t=0
dt=0.001
x=10
v=0
tf=10
k=5


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

#execution
while t<tf:
    a=-k*x          # 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