import numpy as np
from matplotlib import pyplot as plt

y = [10, 20, 30, 40, 50, 52, 60, 70, 200]
x = np.random.uniform(low=1, high=200, size=(len(y),))

min_sum, Y = sum(y), 0
for cur_y in y:
    cur_sum = 0
    for j in y:
        # print('%d += abs(%d - %d)' % (cur_sum, cur_y, j))
        cur_sum += abs(cur_y - j)
    print('for y [%d] sum is [%d]' % (cur_y, cur_sum))
    if cur_sum < min_sum:
        min_sum = cur_sum
        Y = cur_y
print(min_sum, Y)

plt.axhline(y=Y)
plt.scatter(x, y)

low_middle = int((len(y) - 1) / 2)  # floor
median = (y[low_middle] + y[low_middle + 1]) / 2
print(median)
plt.axhline(y=median, color='r')

plt.show()

'''
רק אחרי שמימשתי הבנתי שהסכום המינימלי הוא ממוצע בין האמצעיים, חציון
מזכיר לי חישוב של רגרסיה
'''