#Représentation des vecteurs vitesse du système
N=14 #Nombre de positions
dt=0.040
echelle=0.02 #coefficient pour ajuster la longueur des vecteurs
Vx=[] # crée une liste vide nommée Vx
Vy=[] # crée une liste vide nommée Vx

Vx.append((x[1]- x[0])/dt) # ajoute la valeur entre parenthèse à la fin de la liste Vx
Vy.append((y[1]- y[0])/dt) # ajoute la valeur entre parenthèse à la fin de la liste Vy
plt.arrow(x[0],y[0],Vx[0]*echelle,Vy[0]*echelle, head_width=0.05, head_length=0.07, fc='lightblue', ec='blue')
plt.text(x[0]-0.05,y[0]+0.05,"M"+str(0)) #Création et placement du nom du point

#Affichage
plt.axis('scaled') #repère orthonormé
plt.show()