Cordic et matlab

[Matlabo]

Bonjour;
J'ai lu la page Wikipédia de l'algorithme de CORDIC et j'ai voulu écrire un truc rapide sur Matlab pour le tester.

theta = pi/2.5; % Par exemple
y = 0;
x = 1;
tt = pi/4;
s = 1;
for i = 0:10
    t = power(2, -i);
    if s < 0
        x = x - t*y;
        y = y + t*x;
    else
        x = x + t*y;
        y = y - t*x;
    end
    
    s = theta - tt; 
    
    if s < 0
        tt = tt - atan(t);
    else
        tt = tt + atan(t);
    end
end

Cosx = x * 0.607
Sinx = y * 0.607

Le problème est que ça ne donne pas vraiment les bonnes valeurs ...
Où est-ce que j'ai bien pu me tromper ?

Merci
-----

[jiherve]

bonsoir,
il faut faire :
xn+1 = xn - t*yn;
yn+1= yn + t*xn;
donc passer par des variables intermédiaires.
JR