| | 68 | |
| | 69 | I used mathematica to numerically solve the equation of motion for a uniform sphere as given in Carroll and Ostlie, and checked the solution was correct against a result they published in that book. For the 1/50 case, the function R(t) looks like, |
| | 70 | |
| | 71 | |
| | 72 | rho(t) looks like |
| | 73 | |
| | 74 | |
| | 75 | Here you see that the R(t) drops to zero at the freefall time. You see the density diverging at this point, given the unphysical condition r=0. The collapse accelerates in time, becomming faster and faster by the end. The density increases everywhere uniformly, and is only a function of time. |
| | 76 | |
| | 77 | For all the plots I made, the simulation time was less than the freefall time of the ambient. Additionally, R(t) never fell inside of the BE sphere. |
| | 78 | |
| | 79 | Here are the results, |