import numpy as np
import matplotlib.pyplot as plt

# 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 et pas de temps en USI
n = 10000
tf = 10*tau
dt = tf/(n-1)

# Listes des instants et tensions d’entree associees
t = [i*dt for i in range(n)]
e = [E0*np.cos(2*np.pi*f*ti) for ti in t]
# Création d’une liste pour les tensions du condensateur
u = [0 for i in range(n)]

# Initialisation
u[0] = u0
# Resolution par la methode d’Euler => ligne 27 à compléter !
for i in range(n-1):
    u[i+1] = 

# Representation graphique => lignes 30 et 35 à compléter

plt.xlabel("t (s)")
plt.ylabel("u (V)")
plt.title("Euler ED1 tension d'un condensateur RC série entrée sinuoïdale")
plt.grid()


# Affichage des differentes durees caracteristiques du probleme
print("tau en s:",tau)
print("T en s:",1/f)
print("dt en s:",dt)