# x,y - наши массивы координат
mi, mj = np.meshgrid(np.arange(len(x)), np.arange(len(y)))
mx=x[mi]-x[mj]
my = y[mi] - y[mj]
sqr = mx**2+my**2
mi[sqr == sqr.max()]
# x,y - наши массивы координат
mi, mj = np.meshgrid(np.arange(len(x)), np.arange(len(y)))
mx=x[mi]-x[mj]
my = y[mi] - y[mj]
sqr = mx**2+my**2
mi[sqr == sqr.max()]