Найти радиус и центр окружности, проходящей через N (N>=3) точек первого множества

using System; 
class point_in_plane 
 { 
 /* int main()*/
static void Main() 
   {class point_in_plane {   
    public:
    float x; 
    float y;
    void input (int i);
    float lenght(point_in_plane T);  
};
 
// ввод координат точек
void point_in_plane::input (int i){
  string buf;
  Console.WriteLine("\n Точка  "<<i<<" \n");
  Console.WriteLine("Введите значение координаты х  ");
  buf = Console.ReadLine();
  float x=Convert.ToInt32(buf);  
  Console.WriteLine("Введите значение координаты y  ");
  buf = Console.ReadLine();
  float y=Convert.ToInt32(buf);
}
 
// расстояние до другой точки
float point_in_plane::lenght(point_in_plane T){
    return sqrt(pow(x-T.x,2)+pow(y-T.y,2));}
 
//определение класса circle_in_plane
class circle_in_plane:public point_in_plane{
    public:
    float r;
    void solve(point_in_plane,point_in_plane,point_in_plane);
};
 
// нахождение окружности по трём точкам
void circle_in_plane::solve(point_in_plane a,point_in_plane b,point_in_plane c){
    const float eps=1e-5;
    point_in_plane centr;
    float Det,t1,t2,a11,a12,a21,a22;
    a11=2*(b.x-a.x);
    a12=2*(b.y-a.y);
    a21=2*(c.x-b.x);
    a22=2*(c.y-b.y);
    Det=a11*a22-a21*a12;
    if (fabs(Det)<eps)
        {x=0;y=0;r=0;}
    else{
        t1=b.x*b.x-a.x*a.x+b.y*b.y-a.y*a.y;
        t2=c.x*c.x-b.x*b.x+c.y*c.y-b.y*b.y;
        x=(t1*a22-t2*a12)/Det;
        y=(a11*t2-a21*t1)/Det;
        centr.x=x;centr.y=y;
        r=centr.lenght(a);
        };
};
 
//число точек внутри окружности
int point_in_circle(circle_in_plane CR,point_in_plane * A,int Len){
    int n=0,i;
    for(i=0;i<Len;i++) if (CR.lenght(A[i])<CR.r) n++;
    return n;
};
 
 
 
int main()
{
    
    setlocale(LC_CTYPE, "Russian");  
    //ввод кол-ва точек первого и второго множеств
    int i,N,M;
    Console.WriteLine("Введите количество элементов  1 множества "); 
    buf = Console.ReadLine();
    int N=Convert.ToInt32(buf);                             
    /*while(cin.fail()|N<3){
        Console.WriteLine("ошибка ввода, количество элементов первого множества не может быть меньше 3,\n повторите ввод  ");
        cin.clear(); 
        _flushall();
        cin >> N; } 
    */  
    Console.WriteLine("Введите количество элементов  2 множества  ");
    buf = Console.ReadLine();
    int M=Convert.ToInt32(buf);        
    /*while(cin.fail()|M<0){
        Console.WriteLine("ошибка ввода, повторите   ");
        cin.clear();
        _flushall();
        cin >> M; 
    }
    */  
    ;
    
    point_in_plane *B =new point_in_plane[M];   
    point_in_plane *A =new point_in_plane[N];
    
    Console.WriteLine("\nВведите координаты точек первого множества \n");
    for(i=0;i<N;i++)
    A[i].input(i+1);
 
    Console.WriteLine("\nВведите координаты точек второго множества \n");
    for(i=0;i<M;i++)
    B[i].input(i+1);
    
    circle_in_plane c;
    int i1,i2,i3;
    
    for(i1=0    ; i1<N-2 ;i1++)
    for(i2=i1+1 ; i2<N-1 ;i2++)
    for(i3=i2+1 ; i3<N   ;i3++)
       {
         c.solve(A[i1],A[i2],A[i3]);
         if (c.r)
         {
           int n1,n2;
           n1=point_in_circle(c,A,N);
           n2=point_in_circle(c,B,M);
           if( n1 == n2) 
           {
              Console.WriteLine("\nX="<<c.x");
              Console.WriteLine("\nY="<<c.y"); 
              Console.WriteLine("\nR="<<c.r");
              Console.WriteLine("\nn="<<n1<<"\n\n");
           };  
         };
       }; 
 
    delete [] A;
    delete [] B;
    _getch(); 
    return 0;
}
 } 
 }