Solution (source code)

= Solution

Taking the first marginal of the normal vector in part (i) gives
$$
\boxed{X_1-Y_1\sim N\bigl(0,\,2(a-c)\bigr).}
$$
Equivalently its <variance> is $a+a-2c$; the <covariance> between the two sons must be retained when taking their difference.