X=[1;0]; A=[3,8;1,3]
for k=1:10
    X=A*X ; 
    afficher ((X(2,1))^2)
end