Utilisation de bibliothèques Numpy/Scipy

Question

Un certain nombre d'algorithmes existent pour résoudre une équation, la plupart d'entre eux sont déjà implantées dans le module Scipy.

On a utiliser \(\tt scipy.integrate.odeint(F3,conditions initiales, t)\)

Il faudra modifier la fonction \(\tt F3(thet,thetp,t)\) pour quel soit utilisable par odeint, en effet les conditions initiales doivent être regroupées dans une liste [angle, dérivé_angle] et retourner la grandeur et sa dérivée dans une liste [angle, dérivé_angle] .

Indice

Solution

1
def F3Scipy(V,t):
2
    if V[1] >=0:
3
        theta=V[1]
4
        thetap=-g/L1*sin(V[0] )-mu/m1*V[1] -0.05/(L1*m1)
5
    else:
6
        theta=V[1]
7
        thetap=-g/L1*sin(V[0])-mu/m1*V[1] +0.05/(L1*m1)
8
    return [theta,thetap]
9
10
sol=odeint(F3Scipy,[q10,0],t)
11
angle= sol[:, 0]