ODE System Solution Maple: Plotting Earth's Orbit with Sun at Origin

Join the discussion
Ask a follow-up here, or get your own question answered by working scientists, mathematicians and engineers — people, not an autocomplete.
Real named experts · corrections over time · the nuance an AI answer skips
3 replies · 2K views
incredibe1999
Messages
2
Reaction score
0
Hey there! I was trying to plot a ODE solution, but am not getting what I should be. It is actually plotting the orbit of Earth with Sun at the origin.
My equations
upload_2018-4-5_18-53-13.png

The (-0.00011847) is GM.
The Initial Conditions:
upload_2018-4-5_18-54-10.png


The plot I get:
upload_2018-4-5_18-54-28.png

Should not I be getting a elliptic/circular plot as the orbital positions of the Earth? Also, here the period is t=1. The position is supposed to be same at t=0 and t=1. Help, please!
 

Attachments

  • upload_2018-4-5_18-53-13.png
    upload_2018-4-5_18-53-13.png
    4 KB · Views: 1,154
  • upload_2018-4-5_18-54-10.png
    upload_2018-4-5_18-54-10.png
    2.2 KB · Views: 1,097
  • upload_2018-4-5_18-54-28.png
    upload_2018-4-5_18-54-28.png
    3.4 KB · Views: 767
Physics news on Phys.org
scottdave said:
What is the code you used to generate this plot?
I used the ODE Analyzer.
 
Well, that is not very helpful for us in finding out what goes wrong.

If I use your picture (please post something one can copy as text next time -- it avoids typos) but take different values, even good old Euler in Excel gets me a plot I expect:

upload_2018-4-8_12-59-20.png


values:
Code:
dt 0,01
x_0 1
y_0 0
vx_0 0
vy_0 1
GM 1
so the expressions look OK (provided I reproduced them correctly)

Now: what are the units for your values (and for time) ?
 

Attachments

  • upload_2018-4-8_12-59-20.png
    upload_2018-4-8_12-59-20.png
    6.5 KB · Views: 707