/////// TP1 Exercice 1 ///////


N=input('Nombre d''itérations ? ')
x0=input('Première coordonnée initiale x0 ? ')
y0=input('Seconde coordonnée initale y0 ? ')

X=iterations(N,x0,y0);
I=0:N;
clf()
figure
subplot(1,2,1)
plot(I,X(1,:),'*')
xlabel('Itération')
ylabel('Première coordonnée')

subplot(1,2,2)
Y=sqrt(X(1,:).^2+X(2,:).^2);
Z=log(Y);
plot(I,Z,'*')
xlabel('Itération')
ylabel('Log de la taille')

N0=input('Première itération du régime linéaire ? ')
[a b]=reglin(N0:N,Z(N0+1:N+1));
plot(I,Z,'*',I,a*I+b,'r')
xlabel('Itération')
ylabel('Log de la taille')

//Pour mettre en défaut choisir x0=1 y0=(1-sqrt(5))/2
//Comparer avec cette donnée ce que l'on voit avec N=20, puis avec N=80