Planet path using netwon's law with forward method
显示 更早的评论
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 个评论
James Tursa
2012-2-18
What code do you have so far?
Walter Roberson
2012-2-18
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
2012-2-19
@Walter: Well, OP has Euler's method listed which isn't going to do all that great with *any* orbit.
回答（1 个）
Gregory Vernon
2012-2-19
0 个投票
You should look into velocity verlet methods.
类别
在 帮助中心 和 File Exchange 中查找有关 Earth and Planetary Science 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!
Translated by 
