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

Click For Summary
SUMMARY

The forum discussion focuses on plotting the Earth's orbit using an ODE solution in Maple, with the Sun positioned at the origin. The user reports issues with the generated plot not reflecting the expected elliptical or circular trajectory, despite using the gravitational constant GM = -0.00011847 and initial conditions. The conversation highlights the need for correct parameter values and units to achieve accurate orbital representations, emphasizing the importance of using tools like the ODE Analyzer and Excel for verification.

PREREQUISITES
  • Understanding of Ordinary Differential Equations (ODEs)
  • Familiarity with Maple software for mathematical modeling
  • Knowledge of gravitational constants and their applications in orbital mechanics
  • Basic skills in data visualization and plotting techniques
NEXT STEPS
  • Research how to set initial conditions for ODEs in Maple
  • Learn about the gravitational constant and its significance in orbital calculations
  • Explore the use of Excel for simulating ODE solutions and plotting trajectories
  • Investigate the role of time units in ODE modeling and their impact on results
USEFUL FOR

Students and professionals in physics, mathematicians working with ODEs, and anyone interested in simulating celestial mechanics and orbital dynamics.

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,131
  • upload_2018-4-5_18-54-10.png
    upload_2018-4-5_18-54-10.png
    2.2 KB · Views: 1,075
  • upload_2018-4-5_18-54-28.png
    upload_2018-4-5_18-54-28.png
    3.4 KB · Views: 744
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: 681

Similar threads

  • · Replies 1 ·
Replies
1
Views
2K
  • · Replies 3 ·
Replies
3
Views
2K
  • · Replies 5 ·
Replies
5
Views
2K
  • · Replies 1 ·
Replies
1
Views
3K
  • · Replies 1 ·
Replies
1
Views
3K
  • · Replies 2 ·
Replies
2
Views
4K
  • · Replies 2 ·
Replies
2
Views
3K
  • · Replies 3 ·
Replies
3
Views
2K
  • · Replies 4 ·
Replies
4
Views
3K
  • · Replies 7 ·
Replies
7
Views
3K