Python orbit Beispiele

Beispiel #1

# -*- coding: utf-8 -*-

import os

# compare Euler's method and the Euler-Cromer method for eccentric orbits
#
# M. Zingale (2013-02-19)

import numpy
import pylab
import math
from orbit import *

# circular orbit
o = orbit(1.0, 0.6)  # eccentricity = 0.6

# period
P = o.keplerPeriod()

histEuler = o.intEuler(0.0125, P)
histEC = o.intEulerCromer(0.0125, P)
#histRK2 = o.intRK2(0.0125, P)

pylab.plot(histEuler.x, histEuler.y, label="Euler's method", color="k")

# mark the Sun
pylab.scatter([0], [0], s=250, marker=(20, 1), color="k")
pylab.scatter([0], [0], s=200, marker=(20, 1), color="y")

# draw a vertical line that the semi-major axis should fall on
yy = numpy.linspace(-2.0, 2.0, 100)

Beispiel #2

# -*- coding: utf-8 -*-

import os

# compare different ODE methods on the circular orbit problem
#
# M. Zingale (2013-02-19)

import numpy
import pylab
import math
from orbit import *

# circular orbit
o = orbit(1.0, 0.0)  # eccentricity = 0

# orbital period
P = o.keplerPeriod()

histEuler = o.intEuler(0.05, P)
histEulerCromer = o.intEulerCromer(0.05, P)
histRK2 = o.intRK2(0.05, P)
histRK4 = o.intRK4(0.05, P)

pylab.plot(histEuler.x, histEuler.y, label="Euler's method", color="k")

# mark the Sun
pylab.scatter([0], [0], s=250, marker=(5, 1), color="k")
pylab.scatter([0], [0], s=200, marker=(5, 1), color="y")

pylab.plot(histEulerCromer.x,

Beispiel #3

Datei: circular_orbits_compare.py Projekt: NicovincX2/Python-3.5

# -*- coding: utf-8 -*-

import os

# compare different ODE methods on the circular orbit problem
#
# M. Zingale (2013-02-19)

import numpy
import pylab
import math
from orbit import *

# circular orbit
o = orbit(1.0, 0.0)   # eccentricity = 0

# orbital period
P = o.keplerPeriod()

histEuler = o.intEuler(0.05, P)
histEulerCromer = o.intEulerCromer(0.05, P)
histRK2 = o.intRK2(0.05, P)
histRK4 = o.intRK4(0.05, P)

pylab.plot(histEuler.x, histEuler.y, label="Euler's method", color="k")

# mark the Sun
pylab.scatter([0], [0], s=250, marker=(5, 1), color="k")
pylab.scatter([0], [0], s=200, marker=(5, 1), color="y")

Beispiel #4

Datei: orbit-eccentric.py Projekt: bt3gl/Numerical-Methods-for-Physics

# compare Euler's method and the Euler-Cromer method for eccentric orbits
#
# M. Zingale (2013-02-19)

import numpy
import pylab
import math
from orbit import *

# circular orbit
o = orbit(1.0, 0.6)  # eccentricity = 0.6

# period
P = o.keplerPeriod()

histEuler = o.intEuler(0.0125, P)
histEC = o.intEulerCromer(0.0125, P)
#histRK2 = o.intRK2(0.0125, P)


pylab.plot(histEuler.x, histEuler.y, label="Euler's method", color="k")

# mark the Sun
pylab.scatter([0],[0],s=250,marker=(20,1),color="k")
pylab.scatter([0],[0],s=200,marker=(20,1),color="y")

# draw a vertical line that the semi-major axis should fall on
yy = numpy.linspace(-2.0, 2.0, 100)
pylab.plot(0.0*yy, yy, ls=":", color="0.5")