Approximating distribution in an array.

approx_distribution.py

import numpy as np

def approx_distribution(array):
    array = np.array(array)
    array.sort()
    radii = np.std(array)
    array_size = len(array)
    max_val = np.max(array)
    min_val = np.min(array)
    container = np.zeros(array_size)
    count = 0
    for i in np.linspace(min_val + radii, max_val - radii, array_size):
        container[count] = len(array[np.where((array > i - radii) & (array < i + radii))])
        count += 1
    return container