Code complété :

import matplotlib.pyplot as plt
# Coller ci-dessous les valeurs obtenues par simulation
t = []
u = []
i = []
# Valeurs initiales
qc = [0]
uc = [0]
ic = [50e-3]
# Nombre de valeurs importées
N = len(t)
Dt = t[N-1]/(N-1)
# Calculs
for k in range(1, N):
qc.append(qc[k-1]+ic[k-1]*Dt)
uc.append(qc[k]/47e-6)
ic.append((6-uc[k])/120)
# Tracé des courbes de la tension
plt.subplot(1,2,1)
plt.plot(t, uc, color="blue", linewidth=2, label="calcul")
plt.plot(t, u, "r--", color="black", linewidth=2, label="simulation")
plt.xlabel("t (en s)")
plt.ylabel("u (en V)")
plt.title("Tension aux bornes du condensateur")
plt.legend()
# Tracé des courbes de l'intensité
plt.subplot(1,2,2)
plt.plot(t, ic, color="red", linewidth=2, label="calcul")
plt.plot(t, i, "r--", color="black", linewidth=2, label="simulation")
plt.xlabel("t (en s)")
plt.ylabel("i (en A)")
plt.title("Intensité du courant")
plt.legend()
# Affichage
plt.show()