# capacité numérique charge d'un condensateur

# importation des bibliothèque
from scipy.integrate import odeint
import numpy as np
import matplotlib.pyplot as plt

# conditions initiales
E=10
u0=0
tau=0.1 #constante de temps du circuit
t0=0
tmax=0.5
nt=100#nombre de points de discrétisation
Te=(tmax-t0)/(nt-1)# pas de discrétisation

t=[t0]
s=[u0]

si=0

for i in range(1,nt):
    ti=i*Te
    si=si+Te/tau*(E-si)
    t.append(ti)
    s.append(si)


plt.xlabel("temps t en s ")
plt.ylabel("tension en V ")
plt.plot(t,s,color='red',label=' avec euler')
plt.grid(True)
plt.legend(loc="best")
plt.title("charge d'un condensateur")
plt.show()