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 ] - [[RY ]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: {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 01/28/2014 16:39:29 -CET-)



[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 AVirtualSpaceTimeTravelMachine [Retour à AVirtualSpaceTimeTravelMachine]]
[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 (c) Jean-François Colonna, 2000-2014.
Copyright (c) France Telecom R&D and CMAP (Centre de Mathématiques APpliquées) UMR CNRS 7641 / Ecole Polytechnique, 2000-2014.