points avec vitesse.sagews

g = 9.81
l = 1
theta = pi/2
actualTheta = theta
T = 2*pi #*sqrt(l/g)#*(1+(theta^2)/16)
vitesse = ((2*theta)/ (T/2)) #rad/s^2*

liste =[(0,0)  for i in range(0,(T/2)*10,1)] #liste pour stocker les positions du pendules 
for time in range(0, (T/2)*10, 1): #on ne peut que utiliser des entiers
    liste[time] = (sin(actualTheta)*l,(1-cos(actualTheta))*l)             #on stocke la position actuel du pendule (x,y)
    actualTheta -= vitesse*0.1    #on ajoute l'angle parcouru en 0.1sec à 'langle acutel
point(liste)                        #affichage