#Imports
from __future__ import division
from visual import *
from physutil import *

#Objects
axel = cylinder(pos=vector(0,-5,0), axis=vector(0,0.5,0), radius=0.1, color=color.white)
platform = cylinder(pos=vector(0,-5,0), axis=vector(0,0.1,0), radius=10, material=materials.marble)

#Parameters and Initial Conditions
I = 1
L = 0

#Time and time step
t=0
tf=10
dt=0.01

#MotionMap/Graph
#Larrow = arrow(color=color.orange)

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

	taunet = 0

	platform.rotate(angle=(L/I)*dt, axis=vector(0,1,0), origin=axel.pos)

#	Larrow.pos = axel.pos
#	Larrow.axis = vector(0,0,0)

	t = t + dt