The center between both points is determined
((x1+x2)/2, (y1+y2)/2)
while (abnormal) normal
(-(y2-y1), (x2-x1))
If we normalize this vector, we get
(-(y2-y1), (x2-x1)) / sqrt((x2-x1)^2+(y2-y1)^2)
and if we combine both, we will find two points
((x1+x2)/2, (y1+y2)/2) +- D * (-(y2-y1), (x2-x1)) / sqrt((x2-x1)^2+(y2-y1)^2)
.