Homework #3
ASTR 3800
Due February 2, 2016

1. Download the python modules from ThinkStats. As a test that you can run somebody else's modules and packages, run the code of chap06soln.py.  Take the  part of the main that reads:
sample = np.power(10, log_sample)

mean, median = density.Summarize(sample)

cdf = thinkstats2.Cdf(sample)

print('cdf[mean]', cdf[mean])

pdf = thinkstats2.EstimatedPdf(sample)

thinkplot.Pdf(pdf)

thinkplot.Show(xlabel='household income',

ylabel='PDF')

and  comment it out.

If you get a blank plot, like I did in Canopy, edit the hinkplot.py module.  Go to line 627&628 and comment out

the Clf lines as shown below.

def Show(**options):

"""Shows the plot.

For options, see Config.

options: keyword args used to invoke various pyplot functions

"""

clf = options.pop('clf', True)

Config(**options)

# if clf:

# Clf()

pyplot.show(options)

Print out your plot and turn it in.

2.  Make a plot of a normalized Gaussian distribution.  Center it at 100 and give it a 15 point standard deviation as for IQ.
Sum up the bins.  How many bins does it take to get the total to be within 10^-6 of unity?

Use your python code to integrate between values.
What IQ has a probability of 10billion to one of being above?  (That's the smartest person in history.)
What is the probability of having an IQ below zero?  (yes, yes, I know its zero.  But mathematically it's not.)