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