from __future__ import division, print_function
from vpython import *
scene.width = scene.height = 800

G = 6.7e-11
mEarth = 6e24
mcraft = 15e3
deltat = 60

Earth = sphere(pos=vector(0,0,0), radius=6.4e6, color=color.cyan)
craft = sphere(pos=vector(-10*Earth.radius, 0,0), radius=1e6,
               color=color.yellow, make_trail=True)
vcraft = vector(0,2e3,0)
pcraft = mcraft*vcraft
t = 0
scene.autoscale = False ##turn off automatic camera zoom

while t < 10*365*24*60*60:
    rate(100)     
    craft.pos = craft.pos + (pcraft/mcraft)*deltat
    t = t+deltat