# Median Absolute Deviation, MAD, is available in SciPy
from scipy.stats import median_abs_deviation

my_array = list(range(0, 10)) + [1000]
mad = median_abs_deviation(my_array)
print(mad) # prints 3.0