boucle for sous matlab

[inviteccca3812]

bonsoir a tous,
voila je voudrai savoir s'il été possible sous matlab d'utiliser la constante d’incrémentation i d'une boucle for sous matlab.
en effet j'ai essayé de réaliser ces instruction et je me heurte toujours a une erreur.

load test.txt

for i=1:2
test()
R1=[test(:,i)*test(:,i)]
R2=[test(:,2*i)*test(:,2*i)];
R3=[test(:,3*i)*test(:,3*i)];
R4=[test(:,4*i)*test(:,4*i)];
R5=[test(:,5*i)*test(:,5*i)];
R6=[test(:,6*i)*test(:,6*i)];
R7=[test(:,7*i)*test(:,7*i)];
R8=[test(:,8*i)*test(:,8*i)];
R9=[test(:,9*i)*test(:,9*i)];
%calcul des coordonnée du repere2
%calcul des coordonnée du point D
Xd=(R2/6 - R1/3 + R3/6 + (- (2*R1^2)/9 + (2*R1*R2)/9 + (2*R1*R3)/9 + (4*R1)/9 - (2*R2^2)/9 + (2*R2*R3)/9 + (4*R2)/9 - (2*R3^2)/9 + (4*R3)/9 - 8/9)^(1/2)/2 + 1/3);
Yd=(R1/6 - R2/3 + R3/6 + (- (2*R1^2)/9 + (2*R1*R2)/9 + (2*R1*R3)/9 + (4*R1)/9 - (2*R2^2)/9 + (2*R2*R3)/9 + (4*R2)/9 - (2*R3^2)/9 + (4*R3)/9 - 8/9)^(1/2)/2 + 1/3);
Zd=(R1/6 + R2/6 - R3/3 + (- (2*R1^2)/9 + (2*R1*R2)/9 + (2*R1*R3)/9 + (4*R1)/9 - (2*R2^2)/9 + (2*R2*R3)/9 + (4*R2)/9 - (2*R3^2)/9 + (4*R3)/9 - 8/9)^(1/2)/2 + 1/3);
D=[Xd,Yd,Zd]
%calcul des coordonnée du point E
Xe=(R5/6 - R4/3 + R6/6 + (- (2*R4^2)/9 + (2*R4*R5)/9 + (2*R4*R6)/9 + (4*R4)/9 - (2*R5^2)/9 + (2*R5*R6)/9 + (4*R5)/9 - (2*R6^2)/9 + (4*R6)/9 - 8/9)^(1/2)/2 + 1/3);
Ye=(R4/6 - R5/3 + R6/6 + (- (2*R4^2)/9 + (2*R4*R5)/9 + (2*R4*R6)/9 + (4*R4)/9 - (2*R5^2)/9 + (2*R5*R6)/9 + (4*R5)/9 - (2*R6^2)/9 + (4*R6)/9 - 8/9)^(1/2)/2 + 1/3);
Ze=(R4/6 + R5/6 - R6/3 + (- (2*R4^2)/9 + (2*R4*R5)/9 + (2*R4*R6)/9 + (4*R4)/9 - (2*R5^2)/9 + (2*R5*R6)/9 + (4*R5)/9 - (2*R6^2)/9 + (4*R6)/9 - 8/9)^(1/2)/2 + 1/3);
E=[Xe,Ye,Ze]
%calcul des coordonnée du point F
Xf =(R8/6 - R7/3 + R9/6 + (- (2*R7^2)/9 + (2*R7*R8)/9 + (2*R7*R9)/9 + (4*R7)/9 - (2*R8^2)/9 + (2*R8*R9)/9 + (4*R8)/9 - (2*R9^2)/9 + (4*R9)/9 - 8/9)^(1/2)/2 + 1/3);
Yf=(R7/6 - R8/3 + R9/6 + (- (2*R7^2)/9 + (2*R7*R8)/9 + (2*R7*R9)/9 + (4*R7)/9 - (2*R8^2)/9 + (2*R8*R9)/9 + (4*R8)/9 - (2*R9^2)/9 + (4*R9)/9 - 8/9)^(1/2)/2 + 1/3); 
Zf=(R7/6 + R8/6 - R9/3 + (- (2*R7^2)/9 + (2*R7*R8)/9 + (2*R7*R9)/9 + (4*R7)/9 - (2*R8^2)/9 + (2*R8*R9)/9 + (4*R8)/9 - (2*R9^2)/9 + (4*R9)/9 - 8/9)^(1/2)/2 + 1/3);
F=[Xf,Yf,Zf]
%reconstruction de l'origine du repere 2 dans 1
O2=[(Xd+1),(Ye-1),(Zf-1)]
end

si quelqu’un pouvait m’aidai , merci.
-----

[inviteb9f49292]

Je ne sais pas si c'est à cause de ça mais "i" c'est l'imaginaire pur sous MATLAB