import matplotlib.pyplot as plt
import numpy as np


E=5

N=1000
fe=50000
Te=1/fe
durée=N*Te
f=100
R=20E3
C=159E-9

t=np.linspace(0,0.04,N)
e=0
for i in range(len(t)):
    e+=((4*E/np.pi)*np.sin((2*i+1)*2*np.pi*f*t))/(2*i+1)

tau=R*C
s=[0]*(len(e))
s[0]=0
for i in range(1,len(e)):
    s[i]=(Te/tau)*(e[i-1]-s[i-1])+s[i-1]
    
plt.plot(t,e,label="signal d'entrée e(t)",linestyle='-',color='b')
plt.plot(t,s,label="signal de sortie s(t)",linestyle='-',color='r')
plt.grid(True)
plt.legend()
plt.show()