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

#Window setup
scene.center = (0,0,0)
scene.range = 2000

#Objects
Station = cylinder(pos=(0,0,0), radius=500, axis=(0,-30,0), material=materials.silver)
DP = sphere(pos=vector(0,0,0), radius=30, color=color.green)
AP = sphere(pos=vector(900,900,0), radius=90, color=color.red)

#Parameters and Initial Conditions
mAP = 1

vAP = vector(0,-500,0)

pAP = mAP*vAP

#Time and time step
t = 0
dt = 0.001

#Graph

#Calculation Loop
while t >= 0:
	rate(500)

	if(mag(DP.pos-AP.pos)>DP.radius):

		AP.pos = AP.pos + (pAP/mAP)*dt

	else:

		print "collided!"

	t = t + dt