Skip to main content

Estimating the radius of the Sun - The medieval way

I still don't quite understand the whole process but I'll talk about what I do! But if you don't want to listen to me talk, you can go ahead and read my reference. And before we measure the radius of the Sun, we need to measure the distance to the Sun. Or atleast that's how it was done in the olden days.

So, to start, Halley envisioned a way to estimate the distance to the Sun, using the transit of Venus across the Sun. Depending on what latitude one observes from, the time taken by Venus to cross the Sun is different. If one records the time taken and has an estimate of the drift velocity of Venus, in radians per sec, one can estimate the distance between the two paths Venus takes across the Sun's surface. How to estimate the drift velocity of Venus you ask? Well, we know that it takes the Earth 365 days to complete 360 degrees and that it takes Venus 242 days to do the same. We can therefore estimate the relative velocity with which Venus drifts away or comes towards us. And using this, one can estimate the distance between Earth and Venus (this is where I am stuck!).

Coming to the distance between the Sun and the Earth, it can be written as the sum of the distance between Sun and Venus and Venus and Earth. The distance between Sun and Venus was estimated to be 0.72 times that of the distance between the Sun and the Earth. How you ask? Kepler's third law of course, which relates the time taken for a celestial body to complete one full revolution to it's distance from the center of revolution. We know the time it takes for the Earth to go around the Sun and astronomers were also able to measure how long it takes Venus to go around the Sun (one more thing I need to look up!). Using this, they were able to relate the distance between the Sun and the Earth and the Earth and Venus. Okay, I need to be more quantitative than qualitative. Sigh...

But well, you understand how I can estimate the radius of the Sun once we have an estimate of the distance to the Sun, because we know that the angular size of the Sun is approximately half a degree.

On a serious note, this is getting interesting. I need to make proper notes of this, after I make sense of it. While I was searching for how astronomers estimated the distance to the Sun, I came across similar estimates and methods used for the rest of the planets. It's all geometry, trigonometry and sheer brilliance! I never realized that the lower rungs of the cosmic distance ladder were this interesting. Actually, I thought the lowest rung was stellar parallax, but then even to reach that step, we need the value of One Astronomical Unit i.e the distance between the Sun and the Earth!

Apart from the link I provided at the beginning of this post, there were a few other articles I found online that were of (a bit of) help. Some of them are 1, 2, 3, 4, 5, 6.

Popular posts from this blog

Animation using GNUPlot

Animation using GNUPlotI've been trying to create an animation depicting a quasar spectrum moving across the 5 SDSS pass bands with respect to redshift. It is important to visualise what emission lines are moving in and out of bands to be able to understand the color-redshift plots and the changes in it.
I've tried doing this using the animate function in matplotlib, python but i wasn't able to make it work - meaning i worked on it for a couple of days and then i gave up, not having found solutions for my problems on the internet.
And then i came across this site, where the gunn-peterson trough and the lyman alpha forest have been depicted - in a beautiful manner. And this got me interested in using js and d3 to do the animations and make it dynamic - using sliders etc.
In the meanwhile, i thought i'd look up and see if there was a way to create animations in gnuplot and whoopdedoo, what do i find but nirvana!

In the image, you see 5 static curves and one dynam…

Pandas download statistics, PyPI and Google BigQuery - Daily downloads and downloads by latest version

Inspired by this blog post :, I wanted to play around with Google BigQuery myself. And the blog post is pretty awesome because it has sample queries. I mix and matched the examples mentioned on the blog post, intent on answering two questions - 
1. How many people download the Pandas library on a daily basis? Actually, if you think about it, it's more of a question of how many times was the pandas library downloaded in a single day, because the same person could've downloaded multiple times. Or a bot could've.
This was just a fun first query/question.
2. What is the adoption rate of different versions of the Pandas library? You might have come across similar graphs which show the adoption rate of various versions of Windows.
Answering this question is actually important because the developers should have an idea of what the most popular versions are, see whether or not users are adopting new features/changes they provide…

Adaptive step size Runge-Kutta method

I am still trying to implement an adaptive step size RK routine. So far, I've been able to implement the step-halving method but not the RK-Fehlberg. I am not able to figure out how to increase the step size after reducing it initially.

To give some background on the topic, Runge-Kutta methods are used to solve ordinary differential equations, of any order. For example, in a first order differential equation, it uses the derivative of the function to predict what the function value at the next step should be. Euler's method is a rudimentary implementation of RK. Adaptive step size RK is changing the step size depending on how fastly or slowly the function is changing. If a function is rapidly rising or falling, it is in a region that we should sample carefully and therefore, we reduce the step size and if the rate of change of the function is small, we can increase the step size. I've been able to implement a way to reduce the step size depending on the rate of change of …