import numpy as np
import matplotlib.pyplot as plt
import scipy.integrate as scin

# Parametres du circuit et condition initiale en USI
R = 1e3
C = 1e-6
tau = R*C
u0 = 0
# Amplitude et frequence de la fem du GBF en USI
E0 = 2
f = 1e3
# Nombre d’iterations, instant final en USI
n = 1000
tf = 10*tau

# Listes des instants
time = np.linspace(0,tf,n)

# Définition de la fonction F(u,t)
def F(u,t):
    return(1/tau*(E0*np.cos(2*np.pi*f*t)-u))

# Resolution par la fonction odeint
u=scin.odeint(F,u0,time)

# Representation graphique
plt.plot(time,u)
plt.xlabel("t (s)")
plt.ylabel("u (V)")
plt.title("Odeint ED1 tension d'un condensateur RC série entrée sinuoïdale")
plt.grid()
plt.show()