import numpy as np
import matplotlib.pyplot as plt #Permet de tracer des graphes
from scipy.integrate import odeint
from scipy.integrate import cumtrapz
# on pose toute les constantes
g=9.81
m=0.175
V0=14
x0=0
y0=1
a=10
# on génére la ligne de temps, l'intervalle de temps et la ligne des x pour le grapgique
durée=15
N=1000
dt=durée/N
t=np.linspace(0,durée,N)
def egalitéV(A,t):
    dAdt=np.array([-0.033*(0.08+2.72*(a+4)**2)*(1/m)*(A[0]**2+A[1]**2),(-g+0.033*(1.5+1.4*a))*(1/m)*(A[0]**2+A[1]**2)])
    return dAdt

V=np.ones((N,2))
normeV=np.ones((N,5))
A0=np.array[(14,0)]
V=odeint(egalitéV,A0,t)
M=np.ones((N,2))
x_calc = cumtrapz(V[:,0], t, initial=x0)
y_calc = cumtrapz(V[:,1], t, initial=y0)
for k in range (N):
     M[k,0]=x_calc[k]
     M[k,1]=y_calc[k]
plt.clf()
plt.plot(M[:,0],M[:,1],linewidth=0.5)
plt.xlim(0,30)
plt.ylim(0,12)
plt.grid()
plt.show()