как занести корни уравнения в массив?

template <class S>
double signum (S &X)
{
if (X > 0)
X = 1;
if (a < 0)
X = -1;
}
struct point
{
    double X;
    double y;
    double T;
};
 
int main(int argc, char *argv[]) 
{ setlocale(LC_ALL, "rus");
int pt;
cout<<"Введите количество точек:  ";
cin>>pt; 
point* ArrayPoint = new [pt];  //создадим дин. массив.
int N=3, i, j;
double  A, omega, Th0, R, Tk,r, h;
cout<< "введите шаг:\n";
cin >> h; 
cout << "введите амплитуду:\n"; 
cin >> A; 
cout << "введите omega:\n"; 
cin >> omega; 
cout << "введите Th0:\n"; 
cin >> Th0; 
cout<< "введите Tk:\n";
cin >> Tk; 
cout<< "введите T0:\n";
cin >> ArrayPoint[i].T; 
cout<< "введите y0:\n";
cin >> ArrayPoint[i].y; 
cout<< "введите X0:\n";
cin >> ArrayPoint[i].X; 
m1: 
cout << "введите R:\n"; 
cin >> R;
if (0<R<1){} 
else { 
cout<< "не правильный R\n"; 
goto m1; 
}
//заполнение массива 
for( int i=0; i<pt; ++i )
   {
       cout<<"введите X0: ";
       cin >> ArrayPoint[i].X;
       cout<<"введите yo: ";
       cin>>ArrayPoint[i].y;
       cout<<"введите T0: ";
       cin>>ArrayPoint[i].T;
   }
{ 
 
ArrayPoint[i].T=2*M_PI/omega*N;
m2:
r=Th0+ArrayPoint[i].X*signum(ArrayPoint[i].X);
 
{
    ArrayPoint[i+1].X=ArrayPoint[i].X+h*ArrayPoint[i].y;
ArrayPoint[i+1].y=ArrayPoint[i].y+h*(sin(r)-A*sin(omega)*ArrayPoint[i].T*cos(r));
 
}
    for( int i=0; i< pt; ++i)
   {
        cout << ArrayPoint[i].y;
       cout << ArrayPoint[i].X;
       cout << ArrayPoint[i].T << endl;
    }
     return 0;
 
}
delete [] ArrayPoint; ArrayPoint=NULL; //не забываем освобождать память!!
getch(); 
return 0;
}

point* ArrayPoint = new point[pt];