In [1]:

    

import csv
from PyGravity import *
base = PyGravity()
base.time_interval = 60*60


    

In [2]:

    

#load paricles from CSV file
base.read_file('./example_data.csv')
for item in base.particle_list:
    print item


    

    




a: Position: [  1.00000000e+10   0.00000000e+00   0.00000000e+00], Velocity: [ 0.      -5.61464  0.     ], Mass: 3.55e+21
b: Position: [ -1.00000000e+10   0.00000000e+00   0.00000000e+00], Velocity: [ 0.       5.61464  0.     ], Mass: 3.55e+21
c: Position: [  0.00000000e+00  -1.00000000e+10   0.00000000e+00], Velocity: [-5.61464  0.       0.     ], Mass: 3.55e+21
d: Position: [  0.00000000e+00   1.00000000e+10   0.00000000e+00], Velocity: [ 5.61464  0.       0.     ], Mass: 3.55e+21

In [3]:

    

for i in range(20):
    base.step_all()
    if i % 5 == 0:
        print (base.currant_time*base.time_interval)/60, 'minutes'
        for item in base.particle_list:
            print item


    

    




0 minutes
a: Position: [  1.00000000e+10  -2.02127040e+04   0.00000000e+00], Velocity: [ -8.16332407e-06  -5.61464000e+00   0.00000000e+00], Mass: 3.55e+21
b: Position: [ -1.00000000e+10   2.02127040e+04   0.00000000e+00], Velocity: [  8.16332407e-06   5.61464000e+00   0.00000000e+00], Mass: 3.55e+21
c: Position: [ -2.02127040e+04  -1.00000000e+10   0.00000000e+00], Velocity: [ -5.61464000e+00   8.16332407e-06   0.00000000e+00], Mass: 3.55e+21
d: Position: [  2.02127040e+04   1.00000000e+10   0.00000000e+00], Velocity: [  5.61464000e+00  -8.16332407e-06   0.00000000e+00], Mass: 3.55e+21
0 minutes
a: Position: [  1.00000000e+10  -1.21276224e+05   0.00000000e+00], Velocity: [ -4.89799444e-05  -5.61464000e+00   0.00000000e+00], Mass: 3.55e+21
b: Position: [ -1.00000000e+10   1.21276224e+05   0.00000000e+00], Velocity: [  4.89799444e-05   5.61464000e+00   0.00000000e+00], Mass: 3.55e+21
c: Position: [ -1.21276224e+05  -1.00000000e+10   0.00000000e+00], Velocity: [ -5.61464000e+00   4.89799444e-05   0.00000000e+00], Mass: 3.55e+21
d: Position: [  1.21276224e+05   1.00000000e+10   0.00000000e+00], Velocity: [  5.61464000e+00  -4.89799444e-05   0.00000000e+00], Mass: 3.55e+21
0 minutes
a: Position: [  1.00000000e+10  -2.22339744e+05   0.00000000e+00], Velocity: [ -8.97965648e-05  -5.61464000e+00   0.00000000e+00], Mass: 3.55e+21
b: Position: [ -1.00000000e+10   2.22339744e+05   0.00000000e+00], Velocity: [  8.97965648e-05   5.61464000e+00   0.00000000e+00], Mass: 3.55e+21
c: Position: [ -2.22339744e+05  -1.00000000e+10   0.00000000e+00], Velocity: [ -5.61464000e+00   8.97965648e-05   0.00000000e+00], Mass: 3.55e+21
d: Position: [  2.22339744e+05   1.00000000e+10   0.00000000e+00], Velocity: [  5.61464000e+00  -8.97965648e-05   0.00000000e+00], Mass: 3.55e+21
0 minutes
a: Position: [  1.00000000e+10  -3.23403264e+05   0.00000000e+00], Velocity: [ -1.30613185e-04  -5.61464000e+00   0.00000000e+00], Mass: 3.55e+21
b: Position: [ -1.00000000e+10   3.23403264e+05   0.00000000e+00], Velocity: [  1.30613185e-04   5.61464000e+00   0.00000000e+00], Mass: 3.55e+21
c: Position: [ -3.23403264e+05  -1.00000000e+10   0.00000000e+00], Velocity: [ -5.61464000e+00   1.30613185e-04   0.00000000e+00], Mass: 3.55e+21
d: Position: [  3.23403264e+05   1.00000000e+10   0.00000000e+00], Velocity: [  5.61464000e+00  -1.30613185e-04   0.00000000e+00], Mass: 3.55e+21

In [3]:

    




    

In [ ]: