直線の式は OP↑ + t(PA↑ + PB↑)
A=(2, 7), B=(2, -2), P=(8, 1) は安産なので
(8, 1) + t((-6, 6) + (-6, -3)) が直線の式
T = 3t と置いて x = 8 - 4T, y = 1 + T
Tを消すと x = 8 - 4(y-1)
x + 4y = 12 が答え
ここから sympy / matplotlib 使う方法は知らん