Numpy

import numpy
 
incomes = numpy.random.normal(27000, 15000, 10000)
numpy.mean(incomes)
 
# %matplotlib inline
import matplotlib.pyplot as pyplot
pyplot.hist(incomes, 50)
pyplot.show()
 
print numpy.mean(incomes) # 27066
print numpy.median(incomes) # 26820
 
incomes = numpy.append(incomes, [1000000000])
 
print numpy.mean(incomes) # 127053
print numpy.median(incomes) # 26822