Skip to main content

All 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,...
user avatar
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; ...
isaac5122's user avatar
  • 847
1 vote
0 answers
303 views

Numerical integration for equations of motion

...
Isaac's user avatar
  • 11
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 ...
isaac5122's user avatar
  • 847
0 votes
1 answer
130 views

Integrating the Planck Function over certain filter ranges [closed]

ClearAll["Global`*"] ...
Michael's user avatar
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 (...
user61882's user avatar
  • 359
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 ...
Markus Roellig's user avatar
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 ...
Reb.Cabin's user avatar
  • 8,693
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 ...
Emlie's user avatar
  • 73
1 vote
0 answers
228 views

Orbit recognition error in planet motion [closed]

...
Narasimham's user avatar
  • 3,198
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 ...
user's user avatar
  • 1,897
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 ...
indigoblue's user avatar
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: ...
RedRover's user avatar
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 ...
dr.blochwave's user avatar
  • 8,798
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}}$...
InquisitiveInquirer's user avatar

15 30 50 per page