D= a1*b2*c3*d4-a1*b2*d3*c4-a1*b3*c2*d4+a1*b3*d2*c4+a1*b4*c2*d3-a1*b4*d2*c3-a2*b1*c3*d4+a2*b1*d3*c4+a2*b3*c1*d4-a2*b3*d1*c4-a2*b4*c1*d3+a2*b4*d1*c3+a3*b1*c2*d4-a3*b1*d2*c4-a3*b2*c1*d4+a3*b2*d1*c4+a3*b4*c1*d2-a3*b4*d1*c2-a4*b1*c2*d3+a4*b1*d2*c3+a4*b2*c1*d3-a4*b2*d1*c3-a4*b3*c1*d2+a4*b3*d1*c2