Bonjour,

J'ai la boucle suivante:

Code : Sélectionner tout - Visualiser dans une fenêtre à part
1
2
3
4
 
for n=2:N-1
    xPoint(n)=a(n-1)/b(n)*max(x(n-1)-x(n),0)^(3/2);
end
Je l'ai vectorisée comme ceci:
Code : Sélectionner tout - Visualiser dans une fenêtre à part
1
2
3
4
5
6
 
a2=a(1:end-1);
b2=b(2:end-1);
x2=x(1:end-2);
x3=x(2:end-1);
xPoint=a./b.*max(x2-x3,0).^(3/2);
Est-ce la meilleur façon de faire ? Par que c'est certes plus rapide que la boucle for mais ça gache de la mémoire avec les matrices intermédiaires.


Merci.