Найти проекцию точки на прямую (С++)

/* файл struct.h: */
#pragma once
struct point
{
   double x,
   y;
};

/* файл struct.h: */
#pragma once
 
struct Point2D
{
   double x, y;
};
 
struct Line2D
{
   double x, y;
   double A, B, C;
};
 
 
/* файл Unit1.cpp: */
#include <vcl.h>
#include <iostream.h>
#include <conio.h>
#include "struct.h"
 
 
int main()
{
Point2D point;
Line2D line;
 
cout << "Zadayte koordinaty tochki (x, y)\n";
cout << "Vvedite x koordinatu tochki ";
cin >> point.x;
cout << "Vvedite y koordinatu tochki ";
cin >> point.y;
cout << "Zadayte koefficienty v uravnenii pryamoy A*x + B*y + C = 0\n";
cout << "Vvedite A koefficient pryamoy ";
cin >> line.A;
cout << "Vvedite B koefficient pryamoy ";
cin >> line.B;
cout << "Vvedite C koefficient pryamoy ";
cin >> line.C;
 
//ïðîâåðÿåì ÷òî îáà êîýô. â óð-èè ïðÿìîé íå=0
if ((line.A == 0) && (line.B == 0))
  {
  cout << "Nekorrektniy vvod";
  getch();
  return 0;
  }
 
double M, N, L;
double b;
 
//ñëó÷àé åñëè êîýôôèöèåíò ïðè x=0 èëè îáà êîýôôèöèåíòà íå=0 (ò.å. íîðìàëüíûé ñëó÷àé)
if( line.A == 0 || ((line.A != 0) && (line.B != 0)) )
{
 
//line.y = -M*line.x - L; - Âûðàçèì ó â óðàâíåíèè ïðÿìîé, ãäå M = (line.A / line.B);
//Òåîðåìà. Ïóñòü äàíû äâå ïðÿìûå  a è b è èçâåñòíû èõ óðàâíåíèÿ ó=kx+b è y2=k2x+b2.
//         È åñëè ïðÿìûå a è b ïåðïåíäèêóëÿðíû, òî k = -1 / k2
 
// â ñëó÷àå åñëè êîýô. ïðè õ=0, òî N=0
if (line.A == 0)
  N = 0; //ò.ê. N = (line.B / line.A) = 0
else
  N = (line.B / line.A);
 
  M = (line.A / line.B);
  L = (line.C / line.B);
 
//ïîäñòàâëÿåì êîîðäèíàòû òî÷êè âî 2-îå óð-èå è íàõîäèì êîýô. b
//line.y2 =  N*line.x + b
b = point.y - N*point.x;
 
//èç 2-óõ ïîëó÷åííûõ óð-èé íàõîäèì x è y - ýòî êîîðäèíàòû ïðîåêöèè
//ïðèðàâíèâàåì ïðàâûå ÷àñòè óð-èé: line.y = -M*line.x - L  è  line.y =  N*line.x + b
//N*line.x + b + M*line.x + L = 0
//(N + M)*line.x = - L - b
//òîãäà line.x = (- L - b) / (N + M)
if((N + M) == 0)
  line.x = point.x;
else
  line.x = (- L - b) / (N + M);
 
line.y = -M*line.x - L;
}
 
 
//ñëó÷àé åñëè êîýôôèöèåíò ïðè y = 0
if( line.B == 0 )
{
 
//line.x = -M*line.y - L; - Âûðàçèì x â óðàâíåíèè ïðÿìîé, ãäå M = (line.B / line.A)
//Òåîðåìà. Ïóñòü äàíû äâå ïðÿìûå  a è b è èçâåñòíû èõ óðàâíåíèÿ x=ky+b è x2=k2y+b2.
//         È åñëè ïðÿìûå a è b ïåðïåíäèêóëÿðíû, òî k = -1 / k2
 
  N = 0; //ò.ê. N = (line.A / line.B) = 0
  M = 0; //ò.ê. M = (line.B / line.A) = 0
  L = (line.C / line.A);
 
//ïîäñòàâëÿåì êîîðäèíàòû òî÷êè âî 2-îå óð-èå è íàõîäèì êîýô. b
//line.x2 =  N*line.x + b
b = point.y - N*point.x;
 
//èç 2-óõ ïîëó÷åííûõ óð-èé íàõîäèì x è y - ýòî êîîðäèíàòû ïðîåêöèè
//ïðèðàâíèâàåì ïðàâûå ÷àñòè óð-èé: line.x = -M*line.y - L  è  line.x =  N*line.y + b
//N*line.y + b + M*line.y + L = 0
//(N + M)*line.y = - L - b
//òîãäà line.y = (- L - b) / (N + M)
if((N + M) == 0)
  line.y = point.y;
else
  line.y = (- L - b) / (N + M);
 
line.x = -M*line.y - L;
}
 
cout << "\nProekciya tochki (" << point.x << ", " << point.y << ") " << "na pryamuyu " << line.A << "x";
if(line.B > 0) cout << "+";
cout << line.B << "y";
if(line.C > 0) cout << "+";
cout << line.C << "=0 - (" << line.x << ", " << line.y << ")";
 
getch();
return 0;
}