And there is another method using the non - commutative algebra of the 2-vectors made of the base-vectors of R5:
1st method was
b = {0, 0, 0, 1, 2}
a = {1, 2, 1, 3, 2}
c = {2, 0, 1, 1, 1}
Norm[b - a] Norm[(c - a) - ((c - a).(b - a)/(b - a).(b - a)) (b - a)]
the 2nd method then is
x1 = b - a
x2 = c - a
jj = Subsets[Range[5], {2}]
vv = Table[
Det[
{Part[x1, jj[[k]]], Part[x2, jj[[k]]]}
],
{k, 1, Length[jj]}
]
Sqrt[vv.vv]