from visual import * track = box(pos=vector(0,-.075, 0), size=(1.0, 0.05, .10)) cart = sphere(pos=vector(0,1.5,0), radius=0.01, color=color.green) cart.m = 0.05 cart.p = cart.m*vector(0, 0, 0) deltat = 1.0/30 t = 0 while cart.pos.y>0: t = t + deltat rate(30) cart.pos = cart.pos + (cart.p/cart.m)*deltat ball= sphere(pos=cart.pos, radius=0.05, color=color.green) Fnet = vector(0, -9.8*cart.m, 0) cart.p = cart.p + Fnet*deltat print 't=', t, 'y=', cart.pos.y