Das ist eigentlich eine Aufgabe, die man standardmäßig mit dem Pythagoras löst unter Verwendung des Skalarproduktes.
Setze $$C(-1,2,1),\: D(5,4,1),\: P(2,4,6)$$Richtungsvektor der Geraden: \(r=D-C = (6,2,0)\)
Hypotenuse: \(c=|P-C| = |(3,2,5)|\)
Kathete (Projektion von c auf r): \(a=\left|\left(c\cdot \frac r{|r|}\right)\frac r{|r|}\right| = \frac{|c\cdot r|}{|r|}\)
Gesuchte Distanz per Pythagoras.
$$b^2 = c^2 -a^2 = |(3,2,5)|^2 - \frac{|(3,2,5)\cdot (6,2,0)|^2}{|(6,2,0)|^2} = \frac{259}{10}$$$$\Rightarrow b=\sqrt{25.9}\approx 5.09$$
Rechnung hier.
Probe hier.