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

grandeurs=["L","C","RP","RN","R","Usat","V0","derivV0","Tmax"]
unites=["H","F","Ω","Ω","Ω","V","V","V/s","s"]

simu1=[0.04,1e-7,1050,708,1416,15,0.01,0.1,0.005]
simu2=[0.04,1e-7,708,708,1416,15,0.01,0.1,0.005]
simu3=[0.04,1e-7,750,708,1416,15,0.1,0.1,0.02]

def simulation(valeurs,nombrePas=2000):
    print("***Paramètres imposés***")
    for i in range(len(valeurs)):
        print(grandeurs[i],"=",valeurs[i],unites[i])
    print("***Paramètres calculés***")
    #COMPLETER LE CODE POUR LES AFFICHAGES
    def ED(vectV,t):
        #CODE A COMPLETER EQUATION d'EVOLUTION
    T=#à compléter....
    Y0=#à compléter....
    #choisir la méthode voulue (mettre # pour désactiver l'autre)
    rep=integr.odeint(ED,Y0, T)
    #rep=Euler(ED,Y0,T)
    #COMPLETER LE CODE POUR LES TRACES
    
def Euler(ED,Y0,T):
    pass #à faire dans un deuxime temps