### Relevant code and results for yesterday's post on 'Estimating the value of \pi'

[]

In [58]:
import numpy as np
import matplotlib.pyplot as plt
%matplotlib inline


In [77]:
a = 100
alist = np.arange(0,a+1,1)
anotherlist = np.zeros(iterations)
count = 0.
iterations = 100
for i in xrange(iterations):
tempx = np.random.choice(alist)
tempy = np.random.choice(alist)
if tempx**2 + tempy**2 < (a**2)+1:
count += 1.
anotherlist[i] = count/(i+1)*4.


In [78]:
print anotherlist[-1]
plt.plot(anotherlist)


3.0



Out[78]:

[<matplotlib.lines.Line2D at 0xa8b9f58c>]


In [85]:
a = 100
alist = np.arange(0,a+1,1)
andanotherlist = []
for iterations in [10+1,10**2+1,10**3+1,10**4+1,10**5+1,10**6+1]:
#anotherlist = np.zeros(iterations)
count = 0.
#iterations = 100
for i in xrange(iterations):
tempx = np.random.choice(alist)
tempy = np.random.choice(alist)
if tempx**2 + tempy**2 < (a**2)+1:
count += 1.
#anotherlist[i] = count/(i+1)*4.
andanotherlist.append(count/iterations*4.)


In [86]:
plt.plot(andanotherlist,'o')


Out[86]:

[<matplotlib.lines.Line2D at 0xa8a019ec>]


In [87]:
a = 1000
alist = np.arange(0,a+1,1)
andanotherlist = []
for iterations in [10+1,10**2+1,10**3+1,10**4+1,10**5+1,10**6+1]:
#anotherlist = np.zeros(iterations)
count = 0.
#iterations = 100
for i in xrange(iterations):
tempx = np.random.choice(alist)
tempy = np.random.choice(alist)
if tempx**2 + tempy**2 < (a**2)+1:
count += 1.
#anotherlist[i] = count/(i+1)*4.
andanotherlist.append(count/iterations*4.)


In [88]:
plt.plot(andanotherlist,'o')


Out[88]:

[<matplotlib.lines.Line2D at 0xa89cffec>]


In []:

### 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…

### on MOOCs.

For those of you who don't know, MOOC stands for Massively Open Online Course.

The internet is an awesome thing. It's making education free for all. Well, mostly free. But it's surprising at the width and depth of courses being offered online. And it looks like they are also having an impact on students, especially those from universities that are not top ranked. Students in all parts of the world can now get a first class education experience, thanks to courses offered by Stanford, MIT, Caltech, etc.

I'm talking about MOOCs because one of my new year resolutions is to take online courses, atleast 2 per semester (6 months). And I've chosen the following two courses on edX - Analyzing Big Data with Microsoft R Server and Data Science Essentials for now. I looked at courses on Coursera but I couldn't find any which was worthy and free. There are a lot more MOOC providers out there but let's start here. And I feel like the two courses are relevant to where I …

### On programmers.

I just watched this brilliant keynote today. It's a commentary on Programmers and the software development industry/ecosystem as a whole.

I am not going to give you a tl;dr version of the talk because it is a talk that I believe everyone should watch, that everyone should learn from. Instead, I am going to give my own parallel-ish views on programmers and programming.
As pointed out in the talk, there are mythical creatures in the software development industry who are revered as gods. Guido Van Rossum, the creator of Python, was given the title Benevolent Dictator For Life (BDFL). People flock around the creators of popular languages or libraries. They are god-like to most programmers and are treated like gods. By which, I mean to say, we assume they don't have flaws. That they are infallible. That they are perfect.
And alongside this belief in the infallibility of these Gods, we believe that they were born programmers. That programming is something that people are born wit…