import numpy as np
import matplotlib.pyplot as plt
from scipy.integrate import odeint


## Variables globales utiles pour l'ensemble du script


A = 3.1
B0 = 1/3
w0 = 1/(10e3*22e-9)
Q = 1/3


# Definition de l'equation differentielle : dY/dt = F(Y,t)
def F(Y,t):
    u2 = 0    # A MODIFIER
    v2 = 0    # A MODIFIER
    return [0,0]    # A MODIFIER

# Condition initiale
Y0 = [0.01,0]

# Resolution numerique
temps = np.linspace(0,0.03,5000)
Y_calcule = odeint(F,Y0,temps)

# Recuperation de la liste des valeurs de u2 au cours du temps
U2 = Y_calcule[:,0]

# Trace
plt.figure()
plt.plot(temps,U2,label='u2(t)')
plt.xlabel('t (s)')
plt.ylabel('tension (V)')
plt.legend()
plt.show()