y[i] = (vBD[j] * x[k]) - (vBD[k] * x[j]);
}
+ std::cout << "v12 = " << vBA[0] << ", " << vBA[1] << ", " << vBA[2] << std::endl;
+ std::cout << "v43 = " << vDC[0] << ", " << vDC[1] << ", " << vDC[2] << std::endl;
+ std::cout << "z = " << vBD[0] << ", " << vBD[1] << ", " << vBD[2] << std::endl;
+ std::cout << "p = " << p[0] << ", " << p[1] << ", " << p[2] << std::endl;
+ std::cout << "x = " << x[0] << ", " << x[1] << ", " << x[2] << std::endl;
+ std::cout << "y = " << y[0] << ", " << y[1] << ", " << y[2] << std::endl;
+
u = (x[XX] * x[XX]) + (x[YY] * x[YY]) + (x[ZZ] * x[ZZ]);
v = (y[XX] * y[XX]) + (y[YY] * y[YY]) + (y[ZZ] * y[ZZ]);