Carte d'attracteur

[invitee495456e]

Bonjour tout le monde! Je n'arrive pas à mettre en graphique les autres que les points où j = 0. Le graphique ne semble afficher que les points sortis par j=0, d'où une liasse de points autour de x=0.

Voici mon code source:

#include <stdio.h>
#include <plplot.h>
#include <math.h>
int main (void)
{
  float a, da, xj, xk, xkp1, xmin, xmax, ymin, ymax;
  double xp[1], yp[1], xpb[1], ypb[1];
  int niter, k, j;
  niter=150;
  xmin=0., xmax=15., ymin=0., ymax=2.;
  plinit();
  plenv(xmin,xmax,ymin,ymax,0,1);
/* Executable*/
  a=2.1;
  da=0.1;
  for (j=0; j<=10;j++)
    {
    xj=0.1*j;
    while (a <= 15.)
      {
      xk=xj;
      for (k=0; k<niter;k++)
         {
         xkp1=a*(pow(xk,2))*(1.-xk);
         xk=xkp1;
//       printf("%f\n", xk);
         }
//      return;
      xp[0]=a;
      yp[0]=xk;
      a+=da;
    printf("%lf %lf\n", xp[0], yp[0]);
      if(xk<0.000001)
      {
         plpoin(1,xp,yp,1);
      }
      else
      {
         plpoin(1,xp,yp,2);
      }
    }
  }
  plend();
}

Qu'est-ce qui ne fonctionne pas là-dedans?
-----

[invite4492c379]

Hello,

Mmmm pas simple à te comprendre ... ton algo est pas bon à mon avis, surtout si le comportement du programme te semble incohérent.

Je suppose que l'erreur provient de la boucle while, quand tu en ressors, a>15 (ça arrive quand j vaut 0 forcément). Pour toutes les autres itérations de j, a restera supérieur à 15 donc plus de points.

Revois ton algo.

Edit : Et un debuger est souvent utile dans ces cas, même si une trace manuelle te donne rapidement la solution

[invite2d7144a7]

Bonjour,

À quoi j'ajoute : pourquoi ce mélange de float et double ?

Ça ne sert à rien, sauf à amener des conversions fréquentes, qui demandent du temps, et augmentent les problèmes liés à la précision de la représentation des réels dans nos chers ordis.

[jiherve]

bonsoir,
Bonne réponse de Whoami.
Cependant une remarque les flottants appartiennent bien sur aux réels mais ne sont pas les réels c'est un petit ensemble de fractionnaires.
je ne résiste encore pas à fournir ce lien:
http://download.oracle.com/docs/cd/E..._goldberg.html
JR

[A voir en vidéo sur Futura]

[invite2d7144a7]

Bonjour,

En info, on parle de réels, tout en sachant très bien que ce n'est qu'un petit ensemble très limité (et c'est un point auquel la plupart des débutants se heurtent, surtout ceux qui sont autodidactes).

[inviteccac9361]

Bonjour,

Bonjour tout le monde! Je n'arrive pas à mettre en graphique les autres que les points où j = 0.

C'est normal, ta boucle for (j=0; j<=10;j++) ne sert à rien puisque tu utilise la boucle while (a <= 15.)
Si tu voulais l'utiliser, il te faudrait réinitialiser ta valeur de "a" dans la boucle de "j"

Comme ceci :

  a=2.1; // inutile donc 
  da=0.1;
  for (j=0; j<=10;j++)
    {
    a=2.1;
    xj=0.1*j;
    while (a <= 15.)

sinon, une fois que "a" a dépassé 15, à la prochaine iteration de "j" soit j=1 la boucle des "a" ne sera jamais vérifiée "while (a <= 15.)"