Stack Exchange Network

Stack Exchange network consists of 175 Q&A communities including Stack Overflow, the largest, most trusted online community for developers to learn, share their knowledge, and build their careers.

Visit Stack Exchange

Questions tagged [astronomy]

Questions on the astronomical applications of Mathematica.

5
votes
1answer
53 views

Select Specific StarData

The help menu for StarData has many examples such as StarData[EntityClass[ "Star", {EntityProperty["Star", "DistanceFromEarth"] -> TakeSmallest[10]}]] ...
-1
votes
2answers
36 views

a ListAnimate problem

Thanks to @AlexTrounev, I learned how to use ListAnimate, but I still get stuck using it. Before ListAnimate everything is correct. I just cannot successfully use ListAnimate to make the plot. I have ...
0
votes
1answer
52 views
0
votes
1answer
92 views

Animate particle's trajectory

Does anyone know how to animate the trajectory of a particle, but make sure before the particle reaching a certain spot, there's nothing. I mean I want to make an animation that look like https://en....
1
vote
0answers
52 views

Solving the values of Christoffel symbols from a given metric tensor

Is there any easy mathematica package by which I can solve the Christoffel symbols from a given metric tensor?
3
votes
1answer
83 views

Creating the celestial sphere for astronomical demonstrations

Quick question: has anybody produced (and/or willing to produce and share) a nice celestial sphere with some stars and/or constellations? I am traveling , so I cannot do it myself right now, but I ...
1
vote
0answers
44 views

Calculating Christoffel Symbols [duplicate]

I am totally new in the astrophysics field and I need to calculate Christoffel symbols for a model. The line element for the FLRW metric and the geodesic equation is given below: ...
0
votes
2answers
88 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
1answer
78 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:...
12
votes
1answer
436 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 ...
6
votes
1answer
129 views

Does SolarEclipse only cover a limited date range?

If I use the example from the documentation for SolarEclipse and simply replace the years there with older ones, say from around Shakespeare's time, ...
13
votes
2answers
492 views

Processing images from NASA's Juno Jupiter Orbiter

Question Can you help me improve what I’ve done to create an enhanced-color picture from an image provided by NASA’s JunoCam? [My work is shown in the penultimate image of Jupiter in this post.] I’...
2
votes
1answer
74 views

InverseFunction's domain is not the range of the function that was inverted

I wrote a function to calculate how long of a period has past when a planet traverses an angle theta from perihelion in an orbit of a given eccentricity. ...
1
vote
2answers
72 views

How do I pull a multiple pieces of data from a while loop and put into a table

I am wanting to pull multiple pieces of data from a while loop and put it in a list. The example I have added below is a table with inputted data that I found. Im pretty sure the syntax is unrelated ...
1
vote
0answers
67 views

SatelliteData for a hypothetical satellite

SatelliteData appears to have some sort of orbit propagator behind it. It works for existing satellites, but is there any way to use it to predict the orbit of a ...
3
votes
1answer
100 views

Typeset the standard “Moon” symbol

As it says on the tin, I want a way to typeset the astronomical symbol for the moon, which is a crescent. I also want this in such a way that I could integrate it into a button as a label, like the ...
1
vote
1answer
99 views

list of dates for vernal equinox

I'm trying to get a list of dates and times of the vernal equinox for the last 30 years (or some other time interval). I'm flailing about and really have gotten nowhere. When I've tried typing ...
1
vote
1answer
77 views

TrueAnomaly of planet orbit as a function?

It is possible to obtain the true anomaly of a planet this way: ...
1
vote
1answer
780 views

How to convert astronomical RA and Dec from sexagesimal to degree?

Astronomers use right ascension (RA) and declination (Dec) to refer to the sky coordinates of celestial objects. RA and Dec can take be represented degrees or in sexagesimal (base 60). How can I ...
2
votes
1answer
203 views

Calculating with scientific notation

I want to calculate the gravitational force between the Earth and its Moon using Newton's gravitational formula: (G M1 M2)/d^2 and its appropriate physical units. I'd include the code I wrote for ...
4
votes
1answer
281 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 ...
2
votes
0answers
229 views

Plotting orbits of trajectories of bodies sharing the earth's orbit in coordinates centered at the barycenter of the solar system [closed]

I am attempting to try and plot the trajectories of a trio of satellites that are orbiting with the earth-moon system around the solar system's barycenter. Originally, I had the satellites in orbit ...
4
votes
1answer
385 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 ...
6
votes
1answer
143 views

TimeZone Anomaly

I ran into the following situation when I was trying to determine the time of sunrise (using the built-in Sunrise function) for a variable location on a variable ...
11
votes
1answer
379 views

Calculate moonrise

Mathematica has Sunrise[], but no Moonrise[]. But it has MoonPosition[], so I thought it ...
11
votes
3answers
449 views

How to speed up plotting with PlotPoints?

I'm plotting an ellipse 1/(1 - e Cos[x]) with e=0.9 using PolarPlot. Taking $x\in[0,2\pi)$, ...
10
votes
1answer
311 views

Plotting how sunrise times change over the year

I created a somewhat lengthy Python code to show how the time of sunrise varies over a year, resulting in the following figure (daylight savings were disregarded): The dataset can be found here. The ...
1
vote
0answers
197 views
1
vote
0answers
58 views

Is MoonPosition working as expected?

I have noticed some surprising behavior when executing the MoonPosition function. When I provide a range of dates, I get a TemporalData answer: ...
7
votes
1answer
181 views

Mass distribution in NDSolve

I'm trying to model the solar system with NDSolve, basically a numerical n-body solution. I model the initial positions of bodies from a data set thusly: ...
4
votes
0answers
174 views

Predicting Manhattanhenge

At the time of writing, Manhattanhenge is tonight in New York City (ref 1), "an event during which the setting sun is aligned with the east–west streets of the main street grid of Manhattan, New York ...
4
votes
1answer
113 views

How To Improve Mathematica StarData Performance When Compared To AstronomicalData

My tests of using AstronomicalData versus StarData show that StarData requests are 15 times ...
2
votes
1answer
88 views

How To Replace xAxis and yAxis Number Values With An Associated Name

This application works great. The thing is, the associations: {{xAxis, 6}, {4, 5, 6, 7, 8}}, {{yAxis, 7}, {5, 6, 7, 8}}, while functional, are not very intuitive. How can I replace xAxis and yAxis ...
8
votes
0answers
84 views

Extracting, Formatting, and Exporting Mathematica AstronomicalData - Is There an Easier Way?

If any of you have tried to use Mathematica AstronomicalData and its newest replacements like StarData, you may already be familiar with the large data issues, and the almost impossible download times....
1
vote
2answers
89 views

How to Refresh ChartLabels when Manipulate BarChart Take Data Positions Change

See attached sample code at the bottom of this post. Under Manipulate, as I step through BarChart data using Take with different start/end positions, the chart of the data is refreshed, but the ...
0
votes
2answers
151 views

How to Replace Plot Slot Numbers With Variable Names

The format of my data looks like: ...
3
votes
1answer
154 views

Working with StarData Versus AstronomicalData

Using AstronomicalData I could easily do something like: ...
13
votes
5answers
4k views

Beautiful dance of the Earth and Venus around the sun

I have seen very nice dynamic graphics created by users of Ma. I would like to get help whether the following graphics can be created using Ma. I was not able to create the dynamic picture here. ...
6
votes
1answer
383 views

GeoGraphics - mapping the Moon

The link here shows a nice way to depict the map of the Moon in Mathematica. But, is there a way to center the Mollweide projection onto the far side of the Moon?
1
vote
1answer
725 views

How to add moon to orbit Earth that Orbit Sun

I want to add moon(Blue) that orbiting earth(Red). while Earth is orbiting Sun. this is my code so far that consist of earth, sun and moon. how to write a code so that Moon is orbiting earth instead ...
26
votes
1answer
855 views

Interstellar Image Effects

I'd like to know how to create this type of gravitational lensing effect in an image (here are some star fields from the hubble) "I wrote down the equations and I tested them in Mathematica by ...
4
votes
1answer
171 views

Plotting the Equation of Time: SiderealTime and SunPosition

The equation of time describes the discrepancy between two kinds of solar time. The word equation is used in the medieval sense of reconcile a difference. The two times that differ are the apparent ...
17
votes
2answers
4k views

What is the closest star to Earth other than the sun?

How to find the Sun's nearest neighbor star using Mathematica? I tried but it didn't work
5
votes
2answers
296 views

How to find the Andromeda Galaxy

I would like to find the Andromeda Galaxy ($M31$) with the unaided eye. One way is to use the following image and graphic but I can’t find it just by looking in a dark sky. It is hard to detect ...
12
votes
1answer
664 views

Can I create a radar chart of sunrise/sunset times for a year

Can I create a radar chart (or just a simple curve) of sunrise/sunset times, for a given location thougout a year? It would be something like this, but with a smooth curve: (source: hassanhameed.com)...
8
votes
2answers
571 views

Extracting and displaying data from GW150914 hdf5 files

Has anyone looked at the LIGO data using Mathematica or know how to do this? I'm not familiar with the hdf5 format but I've been able to import the data, consisting primarily of 131,072 floating ...
10
votes
1answer
311 views

3D orbits of moons around their respective planets

I have seen a lot of solar system models (both animated and not), but I am trying to replicate the same thing with moons. Either just one planet and its lunar system or the our whole solar system with ...
3
votes
1answer
101 views

Lack of precision in solar declination leads to odd Fourier result

Here's the sun's declination for a 10 year period (roughly 2000-2010): ...
3
votes
1answer
128 views

Convert star name into entity value

I have a list of star names: stan = Select[#, Last[#] < Quantity[153.5, "LightYears"] &] &@ StarData[StarData[], {"Name", "DistanceFromEarth"}] ...
6
votes
2answers
570 views

Series solution of the Lane-Emden equation

I have the following problem: I'd like to show how the Lane-Emden equation would look like if someone would solve it as a series. Here's the code: ...