click to view the MPEG movie (cliquez pour voir le film MPEG)

Bidimensional display of the rounding-off errors when computing the Verhulst dynamics [Visualisation bidimensionnelle des erreurs d'arrondi dans le calcul de la dynamique de Verhulst].




This picture displays the two following Verhulst dynamics (for 1000 iterations):
                    X  = [[R+1]*X   ] - [R*[X   *X   ]]
                     n           n-1         n-1  n-1

Y = [[R+1]*Y ] - [[R*Y ]*Y ] n n-1 n-1 n-1
computed with the same initial conditions:
                    X  = 0.5
                     0

Y = 0.5 0
They are strictly equivalent mathematically speaking and thus the point set {Xn,Yn} should be on a straight line with equation:
                    Y = X
in the bidimensional space. Unfortunately, when the R parameter increases (from 2.0 -picture 01- to 3.0 -picture 16-) the rounding-off errors become huge and give birth to two different numerical Verhulst dynamics. Then, the {Xn,Yn} points appear "anywhere" inside a square.

Each point {Xn,Yn} is displayed as a small sphere. Its color is chosen according to a color scale of increasing luminance Color palette 'cercle.35' : {Blue,Red,Magenta,Green,Cyan,Yellow,White} when 'n' vary from 1 to 1000. Points with a small 'n' value (most of the blue ones) are always on the straight line with equation:
                    Y = X
for:
                    Y  ~ X
                     n    n
during the first iterations:

This double precision (64 bits) computation was made on a SGI O200 under IRIX 6.5 and C compiler release 7.2.1 [see the same computation on a PC under LINUX 2.2 and C compiler release 2.95.2].


[More information (in english/en anglais)]
[Plus d'informations (en français/in french)]


(CMAP28 WWW site: this page was created on 12/18/2000 and last updated on 02/28/2022 11:11:56 -CET-)



[See all related pictures (including this one) [Voir toutes les images associées (incluant celle-ci)]]

[See the following comment(s): Verhulst dynamics [Voir le(s) commentaire(s) suivant(s): dynamique de Verhulst]]
[Please visit the related DeterministicChaos picture gallery [Visitez la galerie d'images DeterministicChaos associée]]
[Please visit the related DeterministicFractalGeometry picture gallery [Visitez la galerie d'images DeterministicFractalGeometry associée]]
[Please visit the related SensitivityToRoundingOffErrors picture gallery [Visitez la galerie d'images SensitivityToRoundingOffErrors associée]]

[Go back to AVirtualMachineForExploringSpaceTimeAndBeyond [Retour à AVirtualMachineForExploringSpaceTimeAndBeyond]]

[The Y2K Bug [Le bug de l'an 2000]]

[Site Map, Help and Search [Plan du Site, Aide et Recherche]]
[Mail [Courrier]]
[About Pictures and Animations [A Propos des Images et des Animations]]


Copyright © Jean-François Colonna, 2000-2022.
Copyright © France Telecom R&D and CMAP (Centre de Mathématiques APpliquées) UMR CNRS 7641 / Ecole Polytechnique, 2000-2022.