Itérer $X_{n+1}=MX_n$

Partie 6 · ★★

Le TP calcule $X_1=MX_0$ , $X_2=MX_1$ … Écris evolution(M, x0, n) qui renvoie la liste $[X_0, X_1, \dots, X_n]$ (donc $n+1$ vecteurs). matvec est fourni. Exemple du TP : $M=\begin{pmatrix}1{,}2&-0{,}3\\-0{,}2&1{,}1\end{pmatrix}$ , $X_0=(100,80)$ donne $X_1=(96,68)$ et $X_2=(94{,}8,\ 55{,}6)$ .

def matvec(A, x):
    return [sum(a * b for a, b in zip(row, x)) for row in A]

def evolution(M, x0, n):
    xs = [x0]
    for _ in range(___):
        xs.append(matvec(M, xs[-1]))
    return xs

Ctrl/⌘+Entrée

Exercice suivant → Quand le modèle devient irréaliste