All Questions
Tagged with astronomy numerical-integration
17
questions
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 ...
0
votes
1
answer
130
views
Integrating the Planck Function over certain filter ranges [closed]
ClearAll["Global`*"]
...
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 (...
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 ...
5
votes
1
answer
578
views
Lagrangian form of Earth-Sun Kepler blows up NDSolve?
I've looked at this 1-D SE question. That question seems to have difficulties around diverging potential functions, and I don't have that. I've looked at this 2-D 3-body demonstration, but that starts ...
4
votes
1
answer
764
views
Orbit followed by a particle around Schwarzschild Black Hole
The following is a equation which describes various possible orbits of a particle around the Schwarzschild black hole spacetime in general relativity. I want to solve it from ...
1
vote
0
answers
228
views
7
votes
2
answers
653
views
Solving the path of Earth around Sun
This maybe isn't a universally helpful question. Maybe a little of a code dump. But here goes.
I'm trying to solve the path Earth moves around Sun from Earths mass, Suns mass, Earths initial velocity ...
5
votes
1
answer
715
views
Using ParametricNDSolve and NMinimize to find accurate interplanetary trajectory arrival periapse radius
I have an interplanetary trajectory simulation that calculates the required velocity (using a Lambert solver) and escape angle for a spacecraft to travel from Earth to Mars. This model does not use a ...
7
votes
1
answer
5k
views
Combining Gravity Turn and Orbit Models
I have a mathematical model for the motion of an orbiting spacecraft about Earth:
...
2
votes
2
answers
412
views
Problem evaluating a complicated integral
Apologies for stating my problem poorly in the first instance, thank you for the help in narrowing down the issue.
Problem
For those interested, it's for the exact unbiased inverse of the ...
9
votes
3
answers
2k
views
RK4 Gravity Simulator
I have the following RK4 solver which splits the two 2nd order ODEs, used to calculate x and y positions under the influence of a gravitating body where $$x''(t)=\frac{G m x(t)}{(x(t)^2+y(t)^2)^{3/2}}$...