It was shown in this answer that the locus of all points whose sum of distances from two points $A$ and $B$ is a constant $(2a)$, is an ellipsoid of revolution whose equation is
$ (p - C_e)^T ( a^2 I - U U^T ) (p - C_e) = a^2 (a^2 - c^2) $
where $p = [x,y,z]^T , U = \dfrac{1}{2} (A - B) $, and $ c = \| U \| $
and $C_e = \dfrac{1}{2} (A + B) $ is the center of the ellipsoid. At the minimum $a$ , the segment $CD$ passes through this ellipsoid at exactly one point.
Now, parametrically,
$ M(t) = C + t \ CD = C + t d $
Substitute $M(t)$ into the equation of the ellipsoid, you'll get an quadratic equation in $t$. Set its discriminant to zero and solve for $a$.
$ (C - C_e + t d )^T ( a^2 I - U U^T ) (C - C_e + t d) = a^2 (a^2 - c^2) $
Set $v = C - C_e $, then
$ (v + t d )^T ( a^2 I - U U^T ) (v + t d) = a^2 (a^2 - c^2) $
This is the quadratic equation
$ a_2 t^2 + a_1 t + a_0 = 0 $
where
$ a_2 = d^T (a^2 I - U U^T ) d $
$ a_1 = 2 d^T (a^2 I - U U^T ) v $
$ a_0 = v^T (a^2 I - U U^T ) v - a^2 (a^2 - c^2 ) $
Set the discriminant to zero, then
$ \Delta = a_1^2 - 4 a_2 a_0 = 0 $
i.e.
$ \bigg( a^2 d^T v - d^T {UU}^T v \bigg)^2 - \bigg( a^2 d^T d - d^T {UU}^T d \bigg) \bigg( a^2 v^T v - v^T {UU}^T v - a^2 (a^2 - c^2) \bigg) = 0 $
Now, let
$ w = a^2$
$ c_1 = d^T v , c_2 = d^T d , c_3 = v^T v , c_4 = d^T {UU}^T v , c_5 = d^T {UU}^T d , c_6 = v^T {UU}^T v $, then
$ ( c_1 w - c_4)^2 - ( c_2 w - c_5 )( c_3 w - c_6 - w (w - c^2) ) = 0 $
Expanding this cubic polynomial, we get
$ b_3 w^3 + b_2 w^2 + b_1 w + b_0 = 0 $
where
$ b_3 = c_2 $
$ b_2 = c_1^2 - c_2 ( c_3 + c^2 ) - c_5 $
$ b_1 = - 2 c_1 c_4 + c_2 c_6 + c_5 (c_3 + c^2) $
$ b_0 = c_4^2 - c_5 c_6 $
Note that $b_0$ is identically zero. Therefore, the cubic polynomial reduces to
$ w ( b_3 w^2 + b_2 w + b_1) = 0 $
Since $w \ne 0 $ (because $a$ cannot be zero), then our $w$ is one of the two roots of the quadratic
$ b_3 w^2 + b_2 w + b_1 = 0 $
To determine which of the two roots is the correct one, we need to keep in mind that $ w - c^2 = a^2 - c^2 $ must be positive.
I've applied the above formulas to the given points, and found that
$ a = 6.538573313 $
So that the minimum sum $MA + MB = 2 a = 13.07714663 $
To find the point $M$, recall that
$ M = C + t \ d $
The value of $t$ is determined from the quadratic
$a_2 t^2 + a_1 t + a_0 = 0 $
The value of $t$ at which this quadratic has its unique root is
$ t = - \dfrac{a_1}{2 a_2} $
Using the equations above for $a_2$ and $a_1$, one gets
$ t = \dfrac{1}{2} $
So that,
$ M = (1, 2, -1) + \dfrac{1}{2} (1, -4, 1) = ( \dfrac{3}{2} , 0 , 0 ) $