ad[i, j] = np.sqrt((mx[i] - mx[j])**2 + (my[i] - my[j])**2)