| b c d | | a c d | | a b d | | a b c |
A = | f g h | B = -| e g h | C = | e f h | D = -| e f g |
| j k l | | i k l | | i j l | | i j k |
Where |...| indicates the determinant.
This follows from Cramer's rule, and the fact that the plane equation is
just a 4D dot-product (Ax+By+Cz+Dw=0).
Also: provided that D is non-zero, you can set D=1 and A,B,C will be the
result of inverting the left-hand 3x3 submatrix and multiplying by
-(d,h,l)^T. If D=0, the 3x3 submatrix won't have an inverse (its
determinant is -D).