• Support PF! Buy your school textbooks, materials and every day products Here!

Kepler Problem

  • Thread starter ice109
  • Start date
  • #1
1,707
5

Homework Statement



plot the orbit of two masses using the equation for an ellipse and coordinate system at one of the foci

Masses Initial Positions Initial Velocities
i mi x1 x2 x3 v_1 v_ 2 v_ 3
1 1 0.651 0.585 -0.238 -0.755 -0.828 -0.865 -0.726
2 0.931 -0.096 0.000 0.357 -0.209 0.107 -0.660


Homework Equations



[tex] \mathbf{r} = \mathbf{r_1} - \mathbf{r_2} [/tex]

[tex]r=\frac{a\cdot(1-\varepsilon^{2})}{1+\varepsilon\cdot\cos\theta};\,\![/tex]

[tex]\varepsilon = \sqrt{1 + \frac{2EL^{2}}{k^{2}\mu}}[/tex]

[tex]L=\left|\mathbf{L} \right|=\left| \mathbf{r} \times \mathbf{p}\right| [/tex]

[tex] \mu = \frac{1}{\frac{1}{m_{1}} + \frac{1}{m_{2}}} = \frac{m_{1}m_{2}}{m_{1} + m_{2}} [/tex]

[tex] E =\frac{\mu \dot{r}^2}{2} -\frac{k}{r} + \frac{L^2}{2\mu r^2} [/tex]

[tex] k = - G m_1 m_2 [/tex]



The Attempt at a Solution



so all i have to do is plot r a function of theta which seems simple enough. the equation for the conic section only has two parameters i don't know E and L. so from what i understand both are constant of motion so i was so i can just evaluate them at the initial conditions plug them into, the formula and plot? am i right?

to check i calculated the eccentricity: 1.11843. could someone check to make that's one should get for the relative orbit of the light one to the heavier one?
 
Last edited:

Answers and Replies

  • #2
1,707
5
anyone?
 

Related Threads on Kepler Problem

  • Last Post
Replies
8
Views
2K
  • Last Post
Replies
1
Views
2K
Replies
4
Views
6K
Replies
19
Views
11K
Replies
2
Views
582
  • Last Post
Replies
4
Views
6K
Replies
1
Views
994
Replies
1
Views
1K
Replies
0
Views
776
  • Last Post
Replies
4
Views
952
Top