Bonjour,
j'aimerais savoir comment écrire un système d'équations différentielles du deuxième ordre pour l'intégrer ?
Je sais écrire une fonction support d'une équation différentielle:

def rhs(y, t, w0, gamma): #second membre de l'equa diff
    return (y[1], -w0**2*sin(y[0]) - gamma*y[1])


t = np.linspace(0, 1500, 10000)
y = odeint(rhs, (1, 1.1), t, args=(0.5, 0.01))

Voilà comment faire pour travailler sur un système du second ordre?

Merci