si potrebbe fare anche con Sin e Cos? Oppure trovando sempre il coefficiente angolare m, trovare m'=-1/m

e poi prendere
X= Px +/- m' * d
Y= Py +/- m' * d
dopo d è la distanza dal punto Px,Py.

Così si trovano 4 punti, da qui elimino il più distante e il più vicino....
Da tener presente che nel caso in cui P1x=P2x m=inf quindi non funzionerebbe...