% ex2 du TP1: |
% ------------
% la fonction suite est definie dans suite.m

N=10;
A = [5/8 -1/4 1/8; 1/4 0 1/4; 1/8 -1/4 5/8]
b = [1; -1; 1];
u1= [5; 2; -4];
val1=10/3*ones(1,10);
val2=2/3*ones(1,10);
val3=10/3*ones(1,10);

u=suite(A,u1,b,10);

% Graphe:
hold on
plot(u(1,:),'g')
plot(val1,'r')
% plot(u(2,:),'r')
% plot(u(3,:),'b')
hold off

% u1b=[2; 1; 0];
% ub=suite(A,u1b,b,10);
% Graphe:
% hold on
% plot(ub(1,:),'g')
% plot(ub(2,:),'r')
% plot(ub(3,:),'b')
% hold off
% 
% Ab=[5 6 3; -1 5 -1; 1 2 0];
% u11=suite(Ab,u1,b,10);
% Graphe:
% hold on
% plot(u11(1,:),'g') 
% plot(u11(2,:),'r')
% plot(u11(3,:),'b')
% hold off
% 
% u12=suite(Ab,u1b,b,10);
% hold on
% plot(u12(1,:),'g')
% plot(u12(2,:),'r')
% plot(u12(3,:),'b')
% hold off

