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

#Objects
Voyager = sphere(pos=vector(-5e10,-2e10,-3e9), radius=3e9, color=color.green, make_trail=True)
Asteroid = sphere(radius=5e9)

#Parameters and Initial Conditions
vVoyager = vector(1e5, 1e4, 1e4)

#Time and time step
dt = .25e6
t = 0
tf = 1e6

#MotionMap/Graph
trackVoyager = MotionMap(Voyager, tf, 5, markerScale=1e5)

#Calculation Loop
while t < tf:
	rate(10)

	Voyager.pos = Voyager.pos + vVoyager*dt
	
	trackVoyager.update(t,vVoyager)

	t = t + dt