KHELIFI Ahmed Aziz ahmedazizkhelifi
🍉
ahmedazizkhelifi
/ vp2.py
Last active
August 19, 2020 00:54
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| planet = vp.sphere(pos=vp.vector(1,0,0), radius=0.05, color=vp.color.green, | |
| mass = 1, momentum=vp.vector(0,30,0), make_trail=True ) | |
| star = vp.sphere(pos=vp.vector(0,0,0), radius=0.2, color=vp.color.yellow, | |
| mass = 2.0*1000, momentum=vp.vector(0,0,0), make_trail=True) |
ahmedazizkhelifi
/ vp3.py
Created
August 19, 2020 00:18
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| def gravitationalForce(p1,p2): | |
| G = 1 #real-world value is : G = 6.67e-11 | |
| rVector = p1.pos - p2.pos | |
| rMagnitude = vp.mag(rVector) | |
| rHat = rVector / rMagnitude | |
| F = - rHat * G * p1.mass * p2.mass /rMagnitude**2 | |
| return F |
ahmedazizkhelifi
/ vp4.py
Created
August 19, 2020 00:31
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| while True: | |
| vp.rate(500) | |
| #calculte the force using gravitationalForce function | |
| star.force = gravitationalForce(star,planet) | |
| planet.force = gravitationalForce(planet,star) | |
| #Update momentum, position and time | |
| star.momentum = star.momentum + star.force*dt | |
| planet.momentum = planet.momentum + planet.force*dt | |
| star.pos = star.pos + star.momentum/star.mass*dt | |
| planet.pos = planet.pos + planet.momentum/planet.mass*dt |
ahmedazizkhelifi
/ vp5.py
Last active
August 19, 2020 00:59
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| star = vp.sphere(pos=vp.vector(0,0,0), radius=0.2, color=vp.color.yellow, | |
| mass = 1000, momentum=vp.vector(0,0,0), make_trail=True) | |
| planet1 = vp.sphere(pos=vp.vector(1,0,0), radius=0.05, color=vp.color.green, | |
| mass = 1, momentum=vp.vector(0,30,0), make_trail=True) | |
| planet2 = vp.sphere(pos=vp.vector(0,3,0), radius=0.075, color=vp.vector(0.0,0.82,0.33),#RGB color | |
| mass = 2, momentum=vp.vector(-35,0,0), make_trail=True) | |
| planet3 = vp.sphere(pos=vp.vector(0,-4,0), radius=0.1, color=vp.vector(0.58,0.153,0.68), |
ahmedazizkhelifi
/ vp5.py
Created
August 19, 2020 01:13
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| while (True): | |
| vp.rate(500) | |
| #Calculte the force using gravitationalForce function | |
| star.force = gravitationalForce(star,planet1)+gravitationalForce(star,planet2)+gravitationalForce(star,planet3) | |
| planet1.force = gravitationalForce(planet1,star)+gravitationalForce(planet1,planet2)+gravitationalForce(planet1,planet3) | |
| planet2.force = gravitationalForce(planet2,star)+gravitationalForce(planet2,planet1)+gravitationalForce(planet2,planet3) | |
| planet3.force = gravitationalForce(planet3,star)+gravitationalForce(planet3,planet1)+gravitationalForce(planet3,planet2) | |
| #Update momentum, position and time |
OlderNewer