// Permet de dŽterminer l'ordre de convergence : // 1 si g'(r) =/= 0, 2 si g'(r) = 0 et g''(r) =/= 0, ainsi de suite ... // cas g3 : ordre = 3 (s'arrter ˆ n=4, avec x0 = 1 par exple) // cas g2 : ordre = 2 (n=10 ; x0=2) // cas g1 : ordre = 1 (n=10 ; x0=2) clf; ieee(1); exec tp2fonctions.sce; exec tp2PtFixe.sci; n= 10; // Changer n et rŽessayer autant que cela vous plaira. r3=sqrt(3); x0 = 2; xk=iter(g1,x0,n); xxk= log10(abs(xk(1:$-1)-r3)); yk = log10(abs(xk(2:$)-r3)); nk = (1:length(yk)); clf; plot2d(xxk,yk);