Planet path using netwon's law with forward method

1 vue (au cours des 30 derniers jours)
BabyMilo
BabyMilo le 18 Fév 2012
Modifié(e) : Michael le 16 Oct 2013
Given the initial position (ie x1,y1) and the velocity (ie vx,vy) of the planet.
How to plot a graph of the path of the planet around the center point (0,0) assuming the center is Earth, using time step of 10s for 500000s?
formula for forward method:
x(t+dt)=x(t)+v(t)*dt
v(t+dt)=v(t)+a(t)*dt
a=GM/r^2
G = 6.67e-11, M = 5.97e24
Many thanks.
  3 commentaires
Afficher 1 commentaire plus ancien
Walter Roberson
Walter Roberson le 18 Fév 2012
Where is "r" defined?
If I recall correctly, this kind of forward prediction is unable to handle the orbit of Mercury, which has notable influence from drag of the rotational gravitational frame.
James Tursa
James Tursa le 19 Fév 2012
@Walter: Well, OP has Euler's method listed which isn't going to do all that great with *any* orbit.

Connectez-vous pour commenter.

Connectez-vous pour répondre à cette question.

Réponses (1)

Gregory Vernon
Gregory Vernon le 19 Fév 2012
You should look into velocity verlet methods.

Catégories

En savoir plus sur Earth and Planetary Science dans Help Center et File Exchange

Tags

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by