I will provide my proof to this method if you want to, but here it is :
Consider a tetrahedron with vertices labelled 1,2,3,4. Let the sides opposite to each vertex be labelled the same number as that vertex.
Note that if two vectors are (anti)parallel one of them is just a scalar multiple of the other (non-null vectors). This is what you must exploit.
If you pick one vertex and do this, you will get a parametric family of solutions (each coordinate will be of the form at+b where t is a parameter)
Then pick a different vertex and perform the same procedure, but instead plug in these parametric coordinates to find the value of your parameter.
Depending upon the situation, use a suitable form (the first or second expression, both mean the same thing).
This is extremely long though.
EDIT : The summation is badly expressed, sorry about that. What I mean is,
If I label my vertices 1,2,3,4 and i let k=2
Then, due to my labeling, cyclically choose two out of p1, p3 and p4 and find the cross product. I.e. the LHS is equal to
t(p1×p4 + p4×p3 + p3×p1). Always stick to one cycle.
Also the inradius is given by √(|t|).