2nd order DE for planet's density

  • Thread starter Thread starter dinospamoni
  • Start date Start date
  • Tags Tags
    2nd order Density
dinospamoni
Messages
90
Reaction score
0

Homework Statement



I attached the problem because it's easier

Homework Equations





The Attempt at a Solution



The main problem I have with this problem is trying to find the density as a function of radius.
I have been thinking for hours but can't come up with anything.

What I have for the n=0 question:
for clarity, I'm just using x

θ''+(2/x)θ' +1 = 0

is my equation and I'm using mathematica to numerically solve it.
 

Attachments

  • Screen capture00000000000000000 2013-04-17 at 9.39.43 PM.JPG
    Screen capture00000000000000000 2013-04-17 at 9.39.43 PM.JPG
    85.3 KB · Views: 401
  • Screen capture00000000000000000 2013-04-17 at 9.39.50 PM.JPG
    Screen capture00000000000000000 2013-04-17 at 9.39.50 PM.JPG
    49.6 KB · Views: 372
Physics news on Phys.org
The variable theta is defined as the density of the body as a function of radius. Theta is not given explicitly in terms of r, but the values of theta are determined by solving the Lane-Emden ODE subject to the initial conditions specified in the text. As the text describes, the initial conditions themselves are not explicit, and the Lane-Emden equation must be solved in an iterative fashion while manipulating the initial conditions. Basically, follow the recipe given in the text and check the values of theta to see if they give a maximum value when r = 0 and a minimum value when r = R (the surface of the body).
 
There are two things I don't understand about this problem. First, when finding the nth root of a number, there should in theory be n solutions. However, the formula produces n+1 roots. Here is how. The first root is simply ##\left(r\right)^{\left(\frac{1}{n}\right)}##. Then you multiply this first root by n additional expressions given by the formula, as you go through k=0,1,...n-1. So you end up with n+1 roots, which cannot be correct. Let me illustrate what I mean. For this...
Back
Top