Определение взаимного положения плоскостей

  #include <iostream> 
#include <cmath> 
 
using namespace std; 
 
int main() 
{ 
 setlocale(LC_ALL, "Russian"); 
 
 float a1, 
  a2, 
  b1, 
  b2, 
  c1, 
  c2, 
  d1, 
  d2, 
  l, 
  a3, 
  b3, 
  c3, 
  d3; 
 float deta, 
  detb, 
  detc, 
  det1, 
  det2; 
 
 cout « "Enter a1: " « endl; 
 cin » a1; 
 cout « "Enter b1: " « endl; 
 cin » b1; 
 cout « "Enter c1: " « endl; 
 cin » c1; 
 cout « "Enter d1: " « endl; 
 cin » d1; 
 
 cout « "Enter a2: " « endl; 
 cin » a2; 
 cout « "Enter b2: " « endl; 
 cin » b2; 
 cout « "Enter c2: " « endl; 
 cin » c2; 
 cout « "Enter d2: " « endl; 
 cin » d2; 
 
 // Нормализуем 
 l = sqrt(pow(a1, 2) + pow(b1, 2) + pow(c1, 2)); 
 if (l < 0.0001) 
 { 
  cout « "1 - не плоскость" « endl; 
   
 } 
 a2 /= l; 
 b2 /= l; 
 c2 /= l; 
 d2 /= l; 
 
 l = sqrt(pow(a2, 2) + pow(b2, 2) + pow(c2, 2)); 
 if (l < 0.0001) 
 { 
  cout « "2 - не плоскость" « endl; 
   
 } 
 a2 /= l; 
 b2 /= l; 
 c2 /= l; 
 d2 /= l; 
 
 // Совпадение нормалей? 
 if (abs(a1 - a2) + abs(b1 - b2) + abs(c1 - c2) < 0.001) 
 { 
  if (abs(d1 - d2) < 0.001) 
   cout « "Плоскости совпадают" « endl; 
  else 
   cout « "Плоскости паралельны" « endl; 
  
 } 
 
 // Развернули один вектор наоборот 
 a3 = -a2; 
 b3 = -b2; 
 c3 = -c2; 
 d3 = -d2; 
 
 // Опять проверяем на совпадение нормалей 
 if (abs(a1 - a3) + abs(b1 - b3) + abs(c1 - c3) < 0.001) 
 { 
  if (abs(d1 - d3) < 0.001) 
   cout « "Плоскости совпадают (нормали в разные стороны)" « endl; 
  else 
   cout « "Плоскости паралельны (нормали в разные стороны)" « endl; 
   
 } 
 // Считаем направляющую для линии 
 deta = b2 * c2 - c1 * b2; 
 detb = c1 * a1 - a1 * c2; 
 detc = a1 * b2 - b1 * a2; 
 
 // Дальше везде будут d отрицательным знаком 
 d1 = -d1; 
 d2 = -d2; 
 
 // Решаем относительно z 
 // Пусть z = 0 
 // a1*x + b1*y = -d1 
 // a2*x + b2*y = -d2 
 det1 = d1*b2 - b1*d2; 
 det2 = a1*d2 - d1*a2; 
 
 a3 = det1 / detc; 
 b3 = det2 / detc; 
 
 cout « "Точка " « a3 « ", " « b3 « ", " « ", 0" « endl; 
 cout « "Плоскости совпадают" « endl; 
 cout « " x = " « a3 « " * a*" « deta « endl; 
 cout « " y = " « b3 « " * a*" « detb « endl; 
 cout « " z = 0 * a*" « detc « endl; 
  
 
  
 system("pause"); 
 return 0;