Bonjour,

Je veux étudier une suite récurrente de matrice définie ainsi:

U(n+1) = 1/3 * M * U(n) + M

avec M=[2 -3/4 ; 2 -1/2] et U(0) = M

et calculer U(100) et U(1000).

Mon code est le suivant :

Code : Sélectionner tout - Visualiser dans une fenêtre à part
1
2
3
4
5
6
7
M=[2 -3/4;2 -1/2]
U(0)=M;
for i=0:1000
  U(i+1)=(1/3)*M*U(i)+M;
end
U(100)
U(1000)
L'erreur est la suivante :

Code : Sélectionner tout - Visualiser dans une fenêtre à part
1
2
3
U(0)=M;
       !--error 21
invalid index
Je pense que le code que j'ai fait est bon pour des scalaires, mais alors on fait comment avec des matrices?

Merci