Ich würde zunächst den Normalenvektor der Ebene bestimmen
[1, -3, 1] ⨯ [0, 2, -1] = [1, 1, 2]
Nun brauch ich dazu nur einen Vektor der Senkrecht ist und die gleiche Länge hat.
[x, y, z] * [1, 1, 2] = 0
x + y + 2·z = 0
x = - y - 2·z
x^2 + y^2 + z^2 = 1^2 + 1^2 + 2^2
x^2 + y^2 + z^2 = 6
(-y - 2·z)^2 + y^2 + z^2 = 6
2·y^2 + 4·y·z + 5·z^2 = 6
y = √(3 - 3/2·z^2) - z
x = - (√(3 - 3/2·z^2) - z) - 2·z = - √(3 - 3/2·z^2) - z
Damit haben die Richtungsvektoren die Form
[- √(3 - 3/2·z^2) - z, √(3 - 3/2·z^2) - z, z] + [1, 1, 2]
[1 - √(3 - 3/2·z^2) - z, 1 + √(3 - 3/2·z^2) - z, 2 + z]
Der Betrag von z darf hier nicht über √2 = 1.4 sein. Ich setzte mal für z = -1, 0 und 1 ein
[2 - √6/2, √6/2 + 2, 1], [1 - √3, √3 + 1, 2], [- √6/2, √6/2, 3]
Diese 3 Vektoren sollten also die Bedingung erfüllen. Du solltest also mit P und einem dieser Richtungsvektoren die Gerade aufstellen können. Geprüft habe ich das jetzt aber nicht.
Ich schau mal ob der Winkel mit der Ebene 45 Grad beträgt.
arccos([1, 1, 2] * [- √6/2, √6/2, 3] / (|[1, 1, 2]| * |[- √6/2, √6/2, 3]|)) = 45 Grad
Das sieht also soweit recht gut aus.