Pb génération en matlab

[tpscience]

Bonsoir,

Je suis en train d'établir un programme en Matlab, or dès le début je suis confronté à un problème que je n'arrive pas à résoudre.
Il me semble comprendre d'où ça vient mais je n'arrive pas à trouver une solution.
Voici ce que je définis : un simple vecteur dont la taille n'est pas importante, fixons par exemple une dimension de 5. Or, mes composantes sont définies par rapport à un irrationnel, d'où le manque de précision il me semble.

Voici en résumé les parties importantes :

...

a = 1.46;
b = 1.33;
c = a/b;
d = 5e-5;
e = c*d/(2*(c-1));
f = 2*e;

for k=1:5
    if k==1
        g(k) = f/2;
    else
        g(k) = g(k-1) + f;
    end
end

...

Or, si je demande g(k)-g(k-1)-2*e pour quel que soit la valeur de k, le résultat obtenu n'est pas 0...?!!!

Merci de votre aide.
-----

[lou_ibmix_xi]

Une première remarque, utilisé "e" comme nom de variable n'est pas malin puisque "e" est censé contenir la base du logarithme néperien.

Ton problème vient, à mon avis, que tu n'as pas compris comment fonctionne la représentation en nombres flottants, cet article expose bien la problématique.
une illustration rapide, sous MATLAB essaie le bout de code qui suit:

146/133 - 1.46/1.33

tu t'attendais à avoir 0? Il ya donc un bug dans MATLAB ?
Voilà pourquoi il ne faut pas comparer deux flottants pour l'égalité, ou pourquoi il y a des systèmes de calculs symboliques / formels type maxima...