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