Conditions Initiales
Pour chaque simulation, les conditions initiales auront les mêmes résultats peu importe la force choisie.
Simulation 1
Pour la première simulation, nous imposons une vitesse initiale de rayon Vo :
Nous choissisons arbitrairement:

Nous posons G=1.0 et M correspond à la masse totale soit :

Or précédemment, nous avons definit dans la théorie :

Ce qui nous mène à la vitesse vitesse initiale suivante:

En ce qui concerne les positons, nous fixons les conditions initiales à
Figure 1 Conditions initiales pour N = 8000 particules
Figure 2 Conditions initiales pour N = 1000 particules
Simulation 2
Pour la seconde simulaton, la vitesse initiale de rayon

Nous reprendrons les mêmes conditions sur G et M que pour la simulation 1.
Ce qui implique, la vitesse initiale est: Vo=1.0
Figure 3 Conditions initiales pour N = 1000 particules