All Questions
Tagged with astronomy differential-equations
30
questions
2
votes
1
answer
132
views
Solving a ODE (Peeble's Equation) with initial condition a the maximum input range
I have a differential equation which is,
$$
\frac{d X_e}{dz} = \frac{C_r(z)}{H(z)}\times \frac{1}{(1+z)}\Big[\frac{38\eta T^3}{25\pi^2}\alpha (T)X_e^2 - b(T)(1 - X_e) \Big]
$$
Where
$$
\alpha(T) = \...
1
vote
0
answers
107
views
Path of moon orbiting the Earth [closed]
I'm trying to make an animation showing the path of the moon orbiting around the earth, but I'm having some trouble plotting it. For some reason, Mathematica doesn't give me any error messages but ...
5
votes
0
answers
207
views
Fast N-Body Simulation with CUDA [closed]
There is a Mathematica code for 3 body simulation
...
3
votes
1
answer
500
views
2 body orbits around a barycenter
I'm attempting to plot a barycentric orbit of the sun and Jupiter.
Initial Conditions
...
6
votes
1
answer
789
views
gravitational two body problem for the orbit of the sun and jupiter around their barycenter
This is a simulation of the sun and jupiter orbiting the respective barycenter.
Initial conditions
...
4
votes
1
answer
750
views
N body simulation
I'm trying to solve a 3 body problem consisting of the earth moon and sun.
I'm wanting to show the earth orbiting the sun and the moon orbiting the earth as a result of the gravitational forces.
...
2
votes
1
answer
141
views
NDSolve Maxsteps Reached Error
I'm trying to solve a system of Equations. The code is commented below, I am mainly conserved with the last few lines, specifically the NDSolve line. No matter what I do to the Maxsteps option, the ...
7
votes
1
answer
562
views
Can NDSolve address spherical gravitational collapse?
As far as I know spherical gravitational collapse -of central importance to theoretical physics- is thought to be out of the scope of Wolfram Mathematica.
However there are cases, as in this article,...
0
votes
2
answers
240
views
Plotting a parametric solution from NDSolve
Im trying to animate multiple points around multiple parametric plots which depicts orbital motion of body around a planet. Initial conditons are;
...
1
vote
0
answers
303
views
0
votes
1
answer
861
views
Plotting elliptical orbits using Verlet method
I'm looking to solve a 2 body system consisting of the Earth and a body of mass in orbit.
I'm using the Verlet method to numerically integrate Newton's equation of motion in order to plot an ...
8
votes
1
answer
3k
views
Interpreting Mathematica code on black holes
I am trying to understand the code written down on page 7 of this document (code is in Mathematica) I understand pretty much all of the code on the previous page needed to setup the page 7 code (...
1
vote
2
answers
258
views
Using WhenEvent to calculate the period of an orbit
I'm having trouble trying to have Mathematica tell me when the period of the orbit of my system is. Namely, I'm trying to use Mathematica to find the orbital period of the Earth around the Sun. Here's ...
0
votes
1
answer
217
views
Plotting the orbit of the Earth around the Sun. Error: step size is effectively zero; singularity or stiff system suspected
I've been asked to plot the orbit of the Earth around the sun over a period of three years. I've copied my exact code below. When I try to run it I get a series of errors, the first being:
NDSolve:...
15
votes
1
answer
2k
views
Speed up NDSolve compared to Python (calls to LSODA)
I migrated a numerical model code from Python to Mathematica and am surprised how much faster the Python version runs. Profiling of the Python version tells me that it is about 100 times faster (120 ...