% fonction suite de l'Ex2:
% Calcul nbIt fois: u=A*u1+b

function u= suite(A,u1,b,nbIt)
u=zeros(3,nbIt+1);
u(:,1)=u1;
for k = 1:nbIt+1
    u(:,k+1) = A*u(:,k)+b; 
end