Bonjour,

J'ai essayé de calculer la constante delta de feigenbaum avec ce lien:
Code:
% valeurs initiales
deltai=3.2; % pour i=1
ai_2=0;ai_1=1;bkmoins1=0;bprimekmoins1=0;
for i=2:10
    % premiere approximation ; pour j=0;
    aij=ai_1+(ai_1-ai_2)/deltai;
    % 2) Algorithme de Newton
    for j=1:100
        % calcul de la fraction
        for k=1:(2^i)
            % numerateur
            bk=aij-bkmoins1^2;
            % denominateur
            bprimek=1-2*bprimekmoins1*bkmoins1;
            %mise a jour
            bkmoins1=bk;
            bprimekmoins1=bprimek;
        end;
        aijplus1 = aij - bk/bprimek;
        % mise a jour
        aij=aijplus1;
    end;
    % calcul de delta
    deltai=(ai_1 -ai_2)/(aij-ai_1);
    % mise a jour
    ai_2=ai_1;
    ai_1=aij;
end;
display(deltai);
avec aij=; bkmoins1=; bprimekmoins1= et deltai=.
Mais apparemment il bug dès la deuxième itération et je ne sais pas pourquoi ...