Questions préliminaires Script : def subdivision1(tmax:float, n:int) -> list: Dt = tmax/n # Pas de temps return [k*Dt for k in range(n+1)] print(subdivision1(1., 5)) Renvoie : [0.0, 0.2, 0.4, 0.6, 0.8, 1.0] Script : def subdivision2(tmax:float, Dt:float) -> list: n = int(tmax/Dt) return [k*Dt for k in range(n+1)] print(subdivision2(0.3, 0.07)) Renvoie : [0.0, 0.07, 0.14, 0.21000000000000002, 0.28] (Ignorez l’erreur de python) def resolution(e, U:float, tau:float, tmax:float, Dt:float) -> (list, list): t = subdivision2(tmax, Dt) u = [U] uk = U r = Dt/tau for k in range(len(t)-1): uk = uk + r * (e(t[k])-uk) u.
Plus…