Метод хорд для решения системы нелинейных уравнений

#include <iostream>
#include <cmath>
 
using namespace std;
class SystemOfEquations{
public:
    double EpsX, EpsY;//точность для икс и для игрэк
    double x1,x2,y1,y2,x3,y3;
    double x4,y4; //насколько я понимаю пока, это будет транзитная вещь, ПОТОМУ ЧТО ЕЕ НИГДЕ НЕ ДОЛЖНО БЫТЬ!!!!
    double f4,g4;
    double a,b,q; //для того, чтобы хранить в них значения определителей
    double f1,f2,f3,g1,g2,g3;
    int i,j;
    //зададим функции
    double f(double x, double y)
       {
          // return  x +y -3; // Здесь вводим ту функцию, которую будем эксплуатировать
        return x*x+y*y-4;
       }
    double g(double X, double Y)//разный регистр, так как "не может быть перегружена"
       {
         //  return  X*X+Y*Y-9; // Здесь вводим ту функцию, которую будем эксплуатировать
        return X*X-Y*Y-1;
       }
    //задали функции
 
    void EnterrData() {
         cout<<"Enter EpsX"<<endl;cin>>EpsX;   //лучше извне, наверно
        cout<<"Enter EpsY"<<endl;cin>>EpsY;   //и назвать ф-ию "точность"
        cout<<"Enter x3"<<endl;cin>>x3;       //1
        cout<<"Enter y3"<<endl;cin>>y3;       //5
    }
    void Body(){
       //в этой функции вычисляем на базе введенных x3 и y3 все остальное x2,x1,y2,y1,f1,f2...etc
                    //тестовый вывод на экран
       x1=1.01*x3; /* cout<<endl<<"x1= "<<x1<<endl;
                      cout<<endl<<"y1= "<<y1<<endl;
                      cout<<endl<<"x2= "<<x2<<endl;cout<<endl<<"y2= "<<y2<<endl; */
       y1=y3;
       x2=x3;
       y2=1.01*y3;
}
    void CalcFunc(){
       //вычисляем значение функций
       f3=f(x3,y3); g3=g(x3,y3);
       f2=f(x2,y2); g2=g(x2,y2);
       f1=f(x1,y1); g1=g(x1,y1);
       //тестовый вывод на экран
       /* cout<<"The value of function: "<<endl;
       cout<<"f3= "<<f3<<endl; cout<<"f2= "<<f2<<endl;cout<<"f1= "<<f1<<endl;
       cout<<"g3= "<<g3<<endl; cout<<"g2= "<<g2<<endl;cout<<"g1= "<<g1<<endl;*/
}
 
 
 
    //эти массивы все время будут переопределятся (во всяком случаи я надеюсь), если они будут идти после Body
    void Determinants(){
      double  A[3][3]={ //задали массив, он 3*3, так мы решаем систему из двух уравнений
              {x1,x2,x3}, //массив дет(икс)
              {f1,f2,f2},
              {g1,g2,g3}
                       };
      // вывод определителя на экран
      //тестовый вывод на экран
    /* cout<<endl<<"DetIks "<<endl;
      for(i=0;i<3;i++){
          for(j=0;j<3;j++){
              cout<<A[i][j]<<"  ";
          }
          cout<<endl;
      }*/
      double  B[3][3]={
              {y1,y2,y3}, //массив дет(игрэк)
              {f1,f2,f2},
              {g1,g2,g3}
                       };
      // вывод определителя на экран
      //тестовый вывод на экран
     /*  cout<<endl<<"DetY "<<endl;
      for(i=0;i<3;i++){
          for(j=0;j<3;j++){
              cout<<B[i][j]<<"  ";
          }
          cout<<endl;
      }*/
      double  Q[3][3]={
              {1,1,1}, //массив дет
              {f1,f2,f2},
              {g1,g2,g3}
                       };
      // вывод определителя на экран
      //тестовый вывод на экран
     /* cout<<endl<<"Det  "<<endl;
      for(i=0;i<3;i++){
          for(j=0;j<3;j++){
              cout<<Q[i][j]<<"  ";
          }
          cout<<endl;
      }*/
 
      //Ладно, хоть и кривовато, но посчитаю через формулу
a=A[0][0]*A[1][1]*A[2][2]+A[0][1]*A[1][2]*A[2][0]+A[0][2]*A[1][0]*A[2][1]-A[0][0]*A[1][2]*A[2][1]-A[0][1]*A[1][0]*A[2][2]-A[0][2]*A[1][1]*A[2][0];
b=B[0][0]*B[1][1]*B[2][2]+B[0][1]*B[1][2]*B[2][0]+B[0][2]*B[1][0]*B[2][1]-B[0][0]*B[1][2]*B[2][1]-B[0][1]*B[1][0]*B[2][2]-B[0][2]*B[1][1]*B[2][0];
q=Q[0][0]*Q[1][1]*Q[2][2]+Q[0][1]*Q[1][2]*Q[2][0]+Q[0][2]*Q[1][0]*Q[2][1]-Q[0][0]*Q[1][2]*Q[2][1]-Q[0][1]*Q[1][0]*Q[2][2]-Q[0][2]*Q[1][1]*Q[2][0];
//вычислили значения определителей //можно и красиво сделать через триангуляцию
//тестовый вывод на экран
/*cout<<endl<<endl;
cout<<"The value of determinants: "<<endl;
cout<< "detIks "<<a<<endl;cout<< "detY "<<b<<endl;cout<< "det  "<<q<<endl;*/
}
    void XfourYfour(){
        x4=a/q;y4=b/q;
    }
 
    //если все будет нормально, то эту функцию мы будем вызывать, в main() в конструкции if
    bool Precise(){
     //   cout<<endl<<pow(10,-9)*abs(a);cout<<endl<<pow(10, 9)*q;
     //   cout<<endl<<pow(10,-9)*abs(b);cout<<endl<<pow(10, 9)*q;
        if((pow(10,-9)*abs(a)<pow(10, 9)*q)&&(pow(10,-9)*abs(b)<pow(10, 9)*q))//проверить, возможно тут Лажа!!
            return true;
                    else{
                    //тут надо что-то добавить обязательно!
                    return false;}
    }
    void Rederict(){
        f4=f(x4,y4);g4=g(x4,y4);               //посчитали новую функцию
        x1=x2;x2=x3;x3=x4;y1=y2;y2=y3;y3=y4;   //переопределили иксы и игреки
        f1=f2;f2=f3;f3=f4;g1=g2;g2=g3;g3=g4;   //переопределили функции
    }
 
    bool PreciseEps(){
     //   cout<<endl<<abs((x4-x3)/x4);cout<<endl<<EpsX;
     //   cout<<endl<<abs((y4-y3)/y4);cout<<endl<<EpsY;
        if((abs((x4-x3)/x4)<EpsX)&&(abs((y4-y3)/y4)<EpsY))//проверить, возможно тут тоже Лажа!!//ааа, понятно, нельзя так рано
                                                          //переприсваивать
            return true;
                    else{
         /*  x1=x2;x2=x3;x3=x4;y1=y2;y2=y3;y3=y4;
            f4=f(x4,y4);g4=g(x4,y4);
            f1=f2;f2=f3;f3=f4;g1=g2;g2=g3;g3=g4;*/
                    return false;}
    }
 
 
};
 
int main( )
{
    int y=1;
    SystemOfEquations W;
    W.EnterrData();//ввели данные
 
    for(;;){
        W.Body();          //посчитали x1,x2,y1,y2
        W.CalcFunc();      //посчитали знач. функций f1,f2,f3,g1,g2,g3
        W.Determinants();  //задали и посчитали дерминанты
        W.XfourYfour();    //посчитали x4, y4
 
     //теперь заходим в цикл, который откорректирует нач. приближение
while(W.Precise()){   //если выполняется точность, то выйти, если нет:
        W.Rederict(); //сначала снова переопределяем
        W.Body();     //потом опять считаем x1,x2,y1,y2
        W.Determinants(); //детерминанты снова считаем
        W.XfourYfour();   //и вычисляем x4 и y4
} //и так до тех пор, пока нач. приближение не будет нормальным
 
if(W.PreciseEps()){break;}
W.Rederict();      //переопределили
    }
 
    cout<<endl;//вывели на экран
    cout<<"x= "<<W.x4<<endl;
    cout<<"y= "<<W.y4<<endl;
//while(W.PreciseEps());
}//конец процедуры main()

if((1E-9*abs(a)<1E+9*q)&&(1E-9*abs(b)<1E+9*q))

double  A[3][3]={ //задали массив, он 3*3, так мы решаем систему из двух уравнений
              {x1,x2,x3}, //массив дет(икс)
              {f1,f2,f2},
              {g1,g2,g3}

a=A[0][0]*A[1][1]*A[2][2]+A[0][1]*A[1][2]*A[2][0]+A[0][2]*A[1][0]*A[2][1]-A[0][0]*A[1][2]*A[2][1]-A[0][1]*A[1][0]*A[2][2]-A[0][2]*A[1][1]*A[2][0];

 /*тестовый вывод на экран*/