Какое условие поставить на выходе из цикла в методе Рунге Кутта 4 порядка для системы ДУ 1 порядка?

#include <stdio.h>
#include <iostream>
#include <math.h>
using namespace std;
 
double f1(double x1, double x2, double t)// первая функция 
{
    return(-52*x1-100*x2+exp(-t));               
}
 
double f2(double x1, double t)// вторая функция
{
    return(x1+sin(t));               
} 
 
int main()
{
    setlocale(LC_ALL,"rus");
    
    double h, eps;
    //cout << " Введите  h = ";//шаг
    //cin >> h;
    int n;
    n=1000;
    //cout << "Введите n = ";//число точек
    //cin >> n;
    double *x1 = new double[n];
    double *x2 = new double[n];
    double *t1 = new double[n];
    double *t2 = new double[n];
    //cout << "enter x[0] = ";
    //cin >> x[0];
    //cout << "enter y[0] = ";
    //cin >> y[0];
    t1[0] = 0; // начальне значения
    t2[0] = 0; // начальне значения
    x1[0] = 1; // начальне значения
    x2[0] = 0; // начальне значения
    eps = 0.0001; // погрешность
    h = 0.1;   // шаг
    
    for (int i = 1; i < n; i++)
    {
        x1[i] = x1[0] + i*h;
        x2[i] = x2[0] + i*h;
        t1[i] = t1[0] + i*h;
        t2[i] = t2[0] + i*h;
        
        double k1 =  f1(x1[i-1], x2[i-1], t1[i]);
        double k2 =  f1(x1[i-1] + h / 2., x2[i-1] + (h*k1) / 2., t1[i]);
        double k3 =  f1(x1[i-1] + h / 2., x2[i-1] + (h*k2) / 2., t1[i]);
        double k4 =  f1(x1[i-1] + h, x2[i-1] + h*k3, t1[i]);
 
        double m1 =  f2(x1[i-1], t2[i]);
        double m2 =  f2(x1[i-1] + h / 2., t2[i]);
        double m3 =  f2(x1[i-1] + h / 2., t2[i]);
        double m4 =  f2(x1[i-1] + h, t2[i]);
        
        x2[i] = x2[i-1] + (h / 6.)*(k1 + 2 * k2 + 2. * k3 + k4);
        x1[i] = x1[i-1] + (h / 6.)*(k1 + 2 * k2 + 2. * k3 + k4);
 
        cout<< t1[i]<<' '<< x1[i] << ' ' << x2[i] << endl;
    }
    
    delete[] x1;
    delete[] x2;
    delete[] t1;
    delete[] t2;
    system("pause");
    return 0;
}