Alogrithmie sur matlab

[invite9333d8d3]

Bonjour
Dans le cadre d'un projet j'ai tenté de coder la méthode de brent pour trouver le zéro d'une fonction afin de pouvoir faire des tests sur les différents critères d'arrêts des algorithmes.
Toutefois j'ai un soucis car mon algo de converge pas vers la valeure escomptée j'aurais donc aimé savoir si vous pouviez m'aider.

Voici le code en question

 Cliquez pour afficher

% Calcul des racines d'une équation par la méthode de Brent
% INSA-Lyon - 12 Mai 2009
%
function [xs,n,lint] = brent(f,a,b,tol,nmax)
%
fa=f(a); fb=f(b);
if sign(fa) == sign(fb)
disp('Attention, problème de signe !'),return
end
%
if abs(fa) < abs(fb)
tmp=fa;
fa=fb;
fb=tmp;
end
%
c=a;
mflag=1;
n=0;
while (abs(b-a)>tol+eps*abs(a)) & (n<nmax) & (fb~=0)
n=n+1;
lint(n)=abs(b-a);
fc=f(c);
if (fa~=fc) & (fb~=fc)
xs(n)=(a*fb*fc)/((fa-fb)*(fa-fc))+(b*fa*fc)/((fb-fa)*(fb-fc))+(c*fa*fb)/((fc-fa)*(fc-fb));
else
xs(n)=b-fb*(b-a)/(fb-fa)
end

if (xs(n)<(3*a+b)/4) | (xs(n)>b) | (mflag & (abs(xs(n)-b)>=(abs(b-c)/2))) | (~mflag & (abs(xs(n)-b)>=(abs(c-d)/2)))
xs(n)=(a+b)/2;
mflag=1;
else
mflag=0;
end
fs=f(xs(n));
d=c;
c=b;
if fa*fs<0
b=xs(n);
else
a=xs(n);
end
fa=f(a);
fb=f(b);
if abs(fa)<abs(fb)
tmp=fa;
fa=fb;
fb=tmp;
end
end

Je vous remercie d'avance car je ne voit vraiment pas d'ou peut bien provenir mon erreur.
-----

[invitea6f35777]

Salut,

Pourquoi quand abs(fa)<abs(fb) tu échange seulement fa et fb et pas a et b ?