// Ce script permet d'estimer g^(k)(r) lorsque l'ordre est k.
//
// cas g3 : ordre = 3 (s'arrter ˆ n=3, avec x0 = 1 par exple)  g'''(r)=0.09 ???
// cas g2 : ordre = 2 (n=4 ; x0=2)  g''(r)=0.87 ???
// cas g1 : ordre = 1 (n=10 ; x0=2)  g'(r)=0.5 ???
//
ieee(1);
exec tp2fonctions.sce;
exec tp2PtFixe.sci;
n= 10; // Changer n et rŽessayer autant que cela vous plaira. 
r3=sqrt(3);
ordreConv=1;
x0 = 2;
xk=iter(g1,x0,n); 
yk = abs(xk(2:$)-r3) ./ (abs(xk(1:$-1)-r3).^ordreConv) // doit tendre vers g'(r) si diff de 0.
nk = (1:length(yk));
clf; plot2d(nk,yk);