Bonjour,
Pour obtenir ce tracé il faut utiliser le calcul des racines de f(x) par la méthode d'itération du point fixe : Xn+1 = Xn +f(Xn)
Seulement cet algorithme ne converge pas toujours comme c'est le cas avec ton équation :
1 2 3 4 5 6 7 8 9 10 11 12
| x=0:0.01:3;
f=inline('x-exp(sin(x))');
plot(x,-f(x)+x,'r',x,x,'k')
axis equal
hold on
x0=2.5;
for n=1:50
xn=-f(x0)+x0;
line([x0 x0 xn x0],[xn x0 xn xn]);
x0=xn;
end
x0 |
Si tu veux voir la convergence tu peux prendre comme fonction :
f=inline('x-exp(0.8*sin(x))');
Partager