# -*- coding: utf-8 -*-
"""
Created on Mon Aug 29 16:20:13 2022

@author: benja
"""

import matplotlib.pyplot as plt # pour le tracé graphique
from math import floor                # pour la partie entiere : floor()

## paramètres du circuit
R = 10**3                  # en Ohm
C = 3*10**-7               # en Farad
tau = R*C                  # pratique pour les calculs

## définition de l'intervalle d'étude
t0 = 0                     # instant initial
tf = 0.003                 # instant final
n = 1000                   # nombre d'itérations souhaitées 
h =               # à compléter : pas de calculs


## définition de la tension d'excitation
# cas d'une dent de scie periodique avec e(0)=-a et e(T)=a
def e(t):
    a = 2                         # amplitude en V
    T = 0.001                     # periode en s
    t_translat =            # a completer : pour se ramener dans [0,T]
    return                  # a completer : valeur du signal
def e_libre(t):
   return 0



## initialisations
u =                       # à compléter : initialisation de la suite des (ui) avec la CI
Lu = []                   # à compléter : initialisation de la liste des tensions
Lt =                      # à compléter : liste des n+1 instants

## boucle methode Euler
for i in range(n):   #il reste n instants à déterminer (i allant de 0 à n-1)
    u =    # à compléter : calcul de la valeur suivante
                   # à compléter : ajout u à la liste des tensions


Le = [e(t) for t in Lt]
plt.plot(Lt,Lu, label = 'u(t)')            # représentation de u(t)
plt.plot(Lt,Le, label = 'e(t)')            # représentation de e(t)
plt.title('Reponse du circuit en fonction du temps')
plt.xlabel('t(s)')
plt.ylabel('tension(V)')
plt.legend()
plt.show()