Построить множество точек на плоскости, симметричных заданном множестве относительно заданной прямой ax + by + c = 0

 #include <iostream>
 #include <cmath>
 #include <cassert>
 
 struct point
 {
   double  x;
   double  y;
   point ( double x, double y )
       : x( x ), y( y )
   {
   }
   point()
   {
     x = 0.0;
     y = 0.0;
   }
 };
 
 //-----------------------------------------------------------
 void reflect( point *original, point *reflection, 
               int _size, double a, double b, double c )
 {
   assert( ( a != 0.0 ) && ( b != 0.0 ) );
   const double sinus = sin( atan( a / b ) );
   const double cosin = cos( atan( a / b ) );
       for ( int i = 0; i < _size; ++i )
       {
          double X  = original[i].x;
          double Y  = original[i].y;
          double x  = ( X + c / a ) * cosin - Y * sinus;
          double y  = ( X + c / a ) * sinus + Y * cosin;
          double xx =  x * cosin - y * sinus - c / a;
          double yy = -x * sinus - y * cosin;
          reflection[i] = point( xx, yy );
       }
 }
 //-----------------------------------------------------------
 int main()
 {
    const double a =  1.0;
    const double b =  3.0;
    const double c =  -5.0;
    point original[] = { point( 0, -3 ), point( 1, 2 ), point( -2, 5 ) };
    const int _size = sizeof( original ) / sizeof( original[ 0 ] );
    std::cout << "Original points: \n";
        for ( int i = 0; i < _size; ++i )
        {
          std::cout << original[i].x << ", " << original[i].y << std::endl;
        }
    point reflection[ _size ];
    reflect( original, reflection, _size, a, b, c );
    std::cout << "\nReflection points relative to the line \n"
              << a << "x + " << b << "y + " << c << std::endl << std::endl;
        for ( int i = 0; i < _size; ++i )
        {
          std::cout << reflection[i].x << ", " << reflection[i].y << std::endl;
        }
    //std::cin.get();
    return 0;
 }