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

AI Thread Summary
The discussion focuses on plotting the Earth's orbit around the Sun using an ODE solution in Maple, with the Sun positioned at the origin. The user is experiencing issues with the plot, expecting an elliptical or circular orbit but receiving an incorrect output. They mention using specific initial conditions and constants, including GM as -0.00011847, and seek clarification on the expected results and code used for plotting. Another participant suggests that using different values in Excel produces the expected plot, prompting a request for clarification on the units of the user's values and time. The conversation highlights the importance of correct parameterization and unit consistency in orbital simulations.
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,075
  • upload_2018-4-5_18-54-10.png
    upload_2018-4-5_18-54-10.png
    2.2 KB · Views: 1,024
  • upload_2018-4-5_18-54-28.png
    upload_2018-4-5_18-54-28.png
    3.4 KB · Views: 713
Physics news on Phys.org
What is the code you used to generate this plot?
 
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: 654

Similar threads

Replies
1
Views
2K
Replies
3
Views
2K
Replies
5
Views
2K
Replies
1
Views
3K
Replies
1
Views
3K
Replies
2
Views
3K
Replies
2
Views
2K
Replies
4
Views
3K
Replies
7
Views
3K
Back
Top