Hier werden die Unterschiede zwischen zwei Versionen gezeigt.
Beide Seiten der vorigen Revision Vorhergehende Überarbeitung Nächste Überarbeitung | Vorhergehende Überarbeitung | ||
ss20:himmelskoerper.comments [2020/09/13 13:55] icks |
ss20:himmelskoerper.comments [2020/09/18 21:38] (aktuell) icks |
||
---|---|---|---|
Zeile 1: | Zeile 1: | ||
[[Himmelsmechanik|← Zurück zur Projektseite]] | [[Himmelsmechanik|← Zurück zur Projektseite]] | ||
- | ===== Code-Kommentierung ===== | + | ===== Code ===== |
+ | ====Kepler-Methode==== | ||
<file python keplermethode.py> | <file python keplermethode.py> | ||
import vpython as vpy | import vpython as vpy | ||
import math | import math | ||
+ | |||
+ | vpy.scene.width = 1920 | ||
+ | vpy.scene.height = 1080 | ||
class Body: #planet/trabant | class Body: #planet/trabant | ||
Zeile 88: | Zeile 92: | ||
</file> | </file> | ||
+ | |||
+ | |||
+ | ====Newton-Methode==== | ||
<file python projektnewton.py> | <file python projektnewton.py> | ||
import vpython as vpy | import vpython as vpy | ||
import math | import math | ||
+ | |||
+ | vpy.scene.width = 1920 | ||
+ | vpy.scene.height = 1080 | ||
class Body: #planet/trabant | class Body: #planet/trabant | ||
Zeile 147: | Zeile 157: | ||
a += G * M[i2]*(X[i2]-X[i])*1/(((X[i2]-X[i]).mag**3)) | a += G * M[i2]*(X[i2]-X[i])*1/(((X[i2]-X[i]).mag**3)) | ||
- | print(X[i], self.himmelskoerper_liste[i].v) | + | #print(X[i], self.himmelskoerper_liste[i].v) |
- | print(leapfrog(X[i], self.himmelskoerper_liste[i].v, a)) | + | #print(leapfrog(X[i], self.himmelskoerper_liste[i].v, a)) |
self.himmelskoerper_liste[i].ob.pos = leapfrog(X[i], self.himmelskoerper_liste[i].v, a)[0] | self.himmelskoerper_liste[i].ob.pos = leapfrog(X[i], self.himmelskoerper_liste[i].v, a)[0] | ||
self.himmelskoerper_liste[i].v = leapfrog(X[i], self.himmelskoerper_liste[i].v, a)[1] | self.himmelskoerper_liste[i].v = leapfrog(X[i], self.himmelskoerper_liste[i].v, a)[1] |