+ double a[3];
+
+ /*First Vector*/
+ a[0] = bbGetInputInX()[1]-bbGetInputInX()[0];
+ a[1] = bbGetInputInY()[1]-bbGetInputInY()[0];
+ a[2] = bbGetInputInZ()[1]-bbGetInputInZ()[0];
+
+ double b[3];
+
+ /*Second Vector*/
+ b[0] = bbGetInputInX()[2]-bbGetInputInX()[0];
+ b[1] = bbGetInputInY()[2]-bbGetInputInY()[0];
+ b[2] = bbGetInputInZ()[2]-bbGetInputInZ()[0];
+
+ double normal[3];
+
+ normal[0]=(a[1]*b[2])-(a[2]*b[1]);
+ normal[1]=(a[2]*b[0])-(a[0]*b[2]);
+ normal[2]=(a[0]*b[1])-(a[1]*b[0]);
+
+ vtkMath::Normalize(normal);
+
+ std::cout << "Normal axis : " << "X: " << normal[0] << " Y: " << normal[1] << " Z: " << normal[2] << std::endl;
+